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line  17, 

for  dv/n  read  hdv 
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line  18, 
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Figure  2  (after  pane  451:  next  to  P^  designation  should  be  He  ,instead  of  He 


CONTENTS 


TABLE  OF  CONTENTS .  i 

LIST  OF  FIGURES .  iii 

TABLE  .  iv 

ABSTRACT .  v 

CHAPTER  1  THE  THEORY  OF  PARAMAGNETIC  RESONANCE  LINE- 
SHAPES  IN  SOLIDS  AT  VERY  LOW  TEMPERATURES 

1 .  1  Introduction .  1 

1 . 2  Outline  of  moment  calculation  when  hv  >  kT  .  4 

1.  3  Application  of  the  moment  formulas  to  interaction  between 

like  neighbors  .  12 

1.4  Line  moments  from  interaction  between  unlike  spins .  18 

1.5  Discussion  of  demagnetizing  effects .  24 

1 . 6  Comment  on  the  fourth  moment  and  the  effect  of  exchange  .  28 

1. 7  Comment  on  McMillan  and  Opechowski's  second  approximation  32 

1. 8  Comment  on  the  possibility  of  magnetostatic  modes  in 

paramagnetic  crystals  at  low  temperature  .  35 

CHAPTER  2  DISCUSSION  OF  EXPERIMENTAL  APPARATUS  AND 

TECHNIQUE 

41 

2.  1  Description  of  microwave  and  low  temperature  systems  .  .  41 

2.  2  Discussion  of  possible  microwave  saturation  . .  49 

2.  3  Discussion  of  heat  contact  between  sample  and  ^He  system  53 
2.4  Heat  leak  through  capillary  .  61 

2.  5  Discussion  of  spectrometer  linearity  and  large-signal 

corrections  .  64 

CHAPTER  3  MEASUREMENTS  ON  NEODYMIUM  ETHYL  SULFATE 

3.  1  Description  of  neodymium  ethyl  sulfate .  71 

3.  2  Measurements  of  line  areas .  74 

3.  3  First  moment  changes .  77 

3.  4  Second  moment  changes .  82 

3.  5  Discussion  of  possible  antiferromagnetic  coupling  in 

neodymium  ethyl  sulfate .  86 

CHAPTER  4  MEASUREMENTS  ON  NICKEL  FLUOSILICATE 

4.  1  Description  of  nickel  fluosilicate  .  93 

4.  2  First  moment  changes .  97 

4.3  Second  moment  changes  .  101 

4.4  Discussion  of  nickel  fluosilicate  results .  105 

-i- 


CHAPTER  5  MEASUREMENTS  ON  COPPER  POTASSIUM  SULFATE 
AND  COPPER  SULFATE 

5.  1  Description  of  copper  potassium  sulfate .  108 

5.  2  First  moment  changes .  Ill 

5.3  Discussion  of  exchange  in  K2Cu(SO^)2*  . 

5.4  Description  of  copper  sulfate^ .  120 

5.  5  Measurements  on  copper  sulfate  and  discussion  of 

observed  magnetostatic  modes  .  123 

CHAPTER  6  RESONANCE  LINE  BROADENING  FROM  VARIATION 
IN  CRYSTALLINE  FIELD 

6.  1  Introduction .  129 

6.2  Measurements  .  132 

6.  3  Discussion . 136 

CONCLUSION .  140 

ACKNOWLEDGMENTS .  141 

APPENDIX  I  .  A1 

APPENDIX  II  .  A7 

REFERENCES  .  Rl 


-ii“ 


I 


LIST  OF  FIGURES 


Figure  Number 

1  Block  diagram  of  waveguide  spectrometer .  44a 

3 

2  Schematic  drawing  of  He  cryostat  .  45a 

3  Detail  showing  sample  mounted  inside  microwave  cavity 

and  in  heat  contact  with  ^He  reservoir .  46a 

4  Reflected  power  from  cavity  as  function  of  absorption  in  sample  .  .  68a 

Sa  Unit  cell  of  Nd  (  C2  H^SO^ ).  •  9H2O  showing  string  of 

nearest  neighbors .  71a 

Sb  Probable  arrangement  of  water  molecules  around  the 

neodymium  ion .  71a 

6  Resonance  lineshape  in  Nd(C2HcSO^),  *  9H2O  as  function 

of  temperature .  74a 

« 

7  Measured  areas  of  the  three  resonance  lines  in 

Nd (  C2HgSO^),  *  9H2O  compared  to  the  theoretical 

absorption  rAios .  76a 

8  Peak  position  of  lines  in  Nd  (  C2  H^SO^  )-  *  9H2O  compared 

to  calculated  first  moment  shifts .  77a 

9  Schematic  representation  of  various  neighbor  arrangement  in 

Nd(C2HgSO.  ),  '  9H2O  showing  how  the  possibility  for  spin^flips 
will  influence  the  energy  levels  of  a  spin  .  78a 

10  Calculated  square  roots  of  second  moments  as  function  of 

temperature  in  Nd(  C2HgSO.  ).  •  9H2O  compared  to  some 
measured  values .  84a 

11  Rhombohedral  unit  cell  of  NiSiF^’  6H2O  .  93a 

12  Measured  resonance  lineshape  in  NiSiF^  *  6H2O  as 

function  of  temperature .  97a 

13  Measured  peak  positions  of  resonance  lines  in  NiSiF^  •  6H2O 

compared  to  calculated  first  moments .  98a 

Figure  13a  Lower,  strong  line  .  98a 

Figure  13b  Upper,  weak  line  .  98a 

-iii- 


Figure  Number 

14  Measured  and  calculated  square  root  of  second  moment  of  lower 

line  in  Ni  Si  *  6  H2  O  when  Hq  /  z  .  1 02a 

1 5  Measured  and  calculated  square  root  of  second  moment  in 
NiSi  F/  ■  6H2O  when  the  magnetic  field  is  in  a  direction  where 

the  two  resonances  overlap .  104a 

16  Shape  of  resonance  line  I  in  K2Cu(SO^)2'  ^H20  as 

function  of  temperature .  112a 

17  First  moment  shift  of  resonance  in  K2Cu(SO^)2‘  6H2O  .  112b 

Figure  17a  Field  Angle  45^  with  respect  to  crystal  axes 

a  and  b  .  112b 

Figure  17b  Field  parallel  to  b-axis .  112b 

18  Shape  of  resonance  line  in  CuSO^  *  SH2O  as  function 

of  temperature .  125a 

19a  Direction  of  magnetic  axes  in  Ni(NH^)2  (SO^)2  ’  6H2O .  130a 

19b  Energy  levels  in  Ni  (  NH^  >2  ( SO^  >2  *  ^H2  0  when  Hq/z  .  130a 


TABLE 

Table  6-1  Measured  linewidths  and  spin -Hamiltonian  parameters 

in  nickel  Tutton  salt .  135 


-iv- 

\ 


TR378 


ABSTRACT 

The  main  part  of  this  report  is  concerned  with  the  shape  of  paramnagnetic 
resonance  lines  under  conditions  where  hv  >  kT.  In  such  a  case  the  spin  levels 
are  not  equally  populated  and  the  line  moments  resulting  from  spin- spin  inter¬ 
actions  change.  From  the  temperature  dependence  of  the  lineshape  it  is  possible 
to  determine  the  strength  and  nature  of  the  coupling  between  the  ions. 

The  condition  hv  >  kT  was  obtained  experimentally  by  cooling  the 
sample  in  a  He cryostat  to  0.  3°K  . 

The  extension  of  the  Van  Vleck  line  moment  theory  to  low  temperature 
by  weighting  the  population  of  each  spin  level  by  the  proper  Boltzmann  factor  has 
been  considered  by  a  number  of  authors,  and  this  theory  is  reviewed  in  Chapter  1. 
The  moment  changes  depend  upon  whether  the  interacting  spins  are  identical  or 
not.  In  the  interpretation  of  the  experiments,  it  is  found  that  two  spins  may  be 
alike  or  unlike,  depending  upon  the  local  field  due  to  their  neighbors  or  to  their 
magnetic  nuclei.  Particular  consideration  is  also  given  to  the  consequences  of 
the  large  magnetization  of  paramagnetic  salts  at  low  temperatures.  Demagneti¬ 
zing  effects  and  the  excitation  of  magnetostatic  modes  are  discussed. 

In  Chapter  2  the  microwave  spectrometer  and  low  temperature  system 

3 

are  described.  The  problem  of  heat-contact  between  the  sample  and  the  He 
reservoir  is  discussed  in  detail.  Approximate  correction  formulas  for  the 
spectrometer  response  are  given  in  the  case  of  strong  absorption  in  the  sample. 

Measurements  on  neodymium  ethyl  sulfate  are  discussed  in  Chapter  3. 
The  temperature  dependence  of  the  absorption  lineshape  of  this  salt  can  be 
almost  completely  explained  on  the  basis  of  dipolar  interactions  alone.  However, 
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there  seems  to  be  evidence  for  some  additional,  weak,  antiferromagnetic  ex¬ 
change  coupling.  Chapter  4  is  concerned  with  measurements  of  nickel  fluosilicate. 

The  exchange  interaction  has  been  found  to  be  isotropic,  ferromagnetic,  and  of 

- 1 8 

magnitude  3.  9  '  10  ergs. 

Because  of  the  hyperfine  interaction  in  CuK^CSO^)^  *  ,it  has  been 

possible  to  obtain  rough  estimates  of  two  different  exchange  constants  in  this  salt. 
This  problem  and  a  discussion  of  the  observation  of  magnetostatic  modes  are 
presented  in  Chapter  5. 

Chapter  6  is  not  directly  related  to  the  rest  of  this  report.  In  it  are 
reported  measurements  on  diluted  nickel  double  sulfates  undertaken  in  an  effort 
to  obtain  a  suitable  salt  for  a  zero  field  maser.  However,  the  linewidths  observ¬ 
ed  were  much  too  large  for  such  a  device,  the  broadening  being  the  result  of 
crystalline  field  variations. 


CHAPTER  1 


THE  THEORY  OF  PARAMAGNETIC  RESONANCE  LINESHAPES 
IN  SOLIDS  AT  VERY  LOW  TEMPERATURE 
1.1.  Introduction. 

Electron  and  nuclear  spins  are  quantized  in  a  magnetic  field. 

In  resonance  experiments  we  induce  transitions  between  the  energy 
levels  of  the  spins  by  a  suitable  microwave  or  radio-frequency  field 
and  observe  the  power  absorbed.  In  a  constant  field  a  single,  non¬ 
interacting  spin  would  have  infinitely  sharp  energy  levels  (except  for  the 
broadening  due  to  the  very  small,  but  finite  probability  of  spontaneous 
emission),  and  if  we  could  observe  this  spin  alone,  the  resonance  line 
would  look  like  a  6-function.  The  actual  lines  observed  in  solids  are  al¬ 
ways  somewhat  broadened,  and  there  are  several  reasons  for  this  : 

a)  The  N»10^^  or  so  spins  in  the  sample  may  not  be  identical. 
The  electric  crystalline  field  each  spin  sees  from  its  neighboring  ions 
may  vary  somewhat  since  the  interionic  distances  can  change  slightly 
because  of  distortions,  impurity  stresses  and  microscopic  cracks  in 
the  solid.  Through  the  spin-orbit  coupling,  this  will  give  a  spread  of 
the  g -value  and  of  any  crystalline  field  splitting  of  the  spin  ground 
state.  Some  measurements  on  the  variation  of  the  zero  field  splitting 
will  be  reported  in  Chapter  6  . 

b)  The  spins  will  interact  among  themselves  through  dipolar 
and  exchange  forces.  Each  spin  sees  an  addition  I , varying  field  from 
the  randomorientationof  all  its  magnetic  dipole  neighbors,  and  this  will 


broaden  the  resonance  lines.  It  is  then  immediately  recognized  that  if  we 
do  the  resonance  experiment  at  a  temperature  where  all  spins  will  be 
frozen  in  their  lowest  state,  all  spins  will  see  the  same  ordering  of  their 
neighbors,  and  there  will  be  no  line  broadening  from  this  cause.  In  the 
next  sections  we  will  review  the  theory  of  the  temperature  dependence  of 
resonance  lineshapes  caused  by  spin- spin  interaction..  In  Chapters  3  -  5  the 
results  of  some  experiments  made  at  30  kMc/s  and  down  to  0.4°K  will 
be  compared  to  the  theory. 

c)  The  spins  may  interact  so  strongly  with  the  lattice  that  this 
causes  broadening  of  the  resonance.  The  spin-phonon  interaction  is  a 
strong  function  of  temperature,  and  it  is  usually  not  difficult  by  cooling 
to  make  the  spin-lattice  relaxation  time  Tj^  so  long  that  this  type  of 
broadening  is  not  important. 

The  shape  of  a  resonance  line  can  be  given  by  a  normalized  ab- 

oo 

sorption  amplitude  f(v)  at  each  frequency  v  so  that  J  f(v)dv=  1  .  To 

0  N 

calculate  f(v)  we  would  have  to  find  the  eigenvalues  of  the  (ZS-t-l)  -  fold 
matrix  for  the  spin  system  and  sum  up  for  each  energy  interval  dv/n 
around  v/n  .  This  is  impossible.  But  it  was  recognized  by  Van  VIeck 
in  a  now  classical  paper  ^  that  we  can  get  some  information  about  the 
line  shape  through  the  moments  defined  as 

oo 

<  v”  >  =  f  v”f(v)dv  (1-1) 

b 

The  moments  can  be  expressed  as  diagonal  sums  of  matrices  in¬ 
cluding  all  spins,  but  since  traces  are  invariant  under  similarity  trans- 


formations,  we  can  calculate  the  trace  in  the  most  convenient  represen¬ 
tation,  namely  when  all  spins  are  individually  diagonalized. 

ri  n  ^  zi 

From  <  V  >  ,  the  moments  <Av>  =  <(v-v)  >  around  the 

unperturbed  resonance  frequency  v  oan  be  foxind.  In  the  Van  Vleck 

treatment,  which  is  a  high  temperature  approximation,  the  odd  moments 

vanish,  <  Av^  >  =  0  for  n  =  odd.  The  second  and  sometimes  the  fourth 

moments  have  been  calculated  and  compared  to  experimental  values  in  a 

number  of  cases. 

The  moment  calculations  can  be  extended  to  low  temperatures 

where  hv  >  kT  by  weighting  each  term  in  the  traces  by  the  proper 

Boltzmann  factor  exp(-hv/kT).  This  has  been  done  by  Bloembergen,^ 

3  4 

Pryce  and  Stevens,  Kambe  and  Usui,  and  more  explicitly  by  McMillan 
5 

and  Opechowski.  It  turns  out  that  the  odd  moments  are  of  ten  finite,  and 
the  first  moment, in  particular, may  show  interesting  temperature  shifts. 
Except  for  the  early  demonstration  by  Bloembergen^  that  the  proton 
resonance  in  CuSO^  *  ^^2^  splits  up  into  a  number  of  components  at 
helium  temperatures  because  of  the  alignment  with  the  field  of  the  Cu^^ 
spins,  there  seems  to  have  been  no  measurements  of  low  temperature 
changes  in  resonance  line  moments. 

The  method  of|Calculating  line  moments  is  exact,  and  in  principle 
it  should  give  all  information  needed  about  the  line  shape  .  But  in 
practice,  it  is  hard  to  calculate  the  higher  moments,  and  since  the 
measured  moments  depend  strongly  upon  the  absorption  far  out  in  the 
weak  wings  of  the  line,  it  is  oftenhard  to  compare  theory  with  experiment. 


4 


Another,  more  direct  approach  to  the  problem  of  resonance  line- 

6  7 

shapes  has  been  attempted  by  Anderson  and  Weiss  ’  and  by  Kubo  and 
8  9  10 

Tomita.  ’  ’  They  are  able  to  derive  the  actual  lineshape  from  the 
Fourier  transform  of  the  relaxation  function  of  the  magnetic  moment. 

By  proper  use  of  Boltzmann  factors  it  can  also  be  shown  by  this  method 
that  the  lines  narrow  at  low  temperature,  as  has  been  done  for  a  special 
case  by  Moriya  and  Obata.^^  But  again  the  calculations  are  complicated 
and  so  many  approximations  necessary,  that  for  quantitative  comparison 
with  experiments  this  approach  is  no  better  than  the  moment  method  in 
most  cases.  Therefore,  only  the  qualitative  conclusions  from  this 
theory  will  be  discussed  briefly  in  Section  1.6. 

1.2.  Outline  of  moment  calculation  when  hv  ^  kT. 

In  this  section  we  review  the  moment  theory  as  given  by  Kambe 

4  5 

and  Usui,  but  we  will  use  the  notation  of  McMillan  and  Opechowski. 

The  theory  can  be  applied  to  nuclear  resonance  as  well  as  to  electron 

resonance,  but  only  in  the  latter  case  can  we  experimentally  realize  the 

condition  hv  m  kT  . 

The  Hamiltonian  for  a  system  of  N  interacting  spins  is 

i  =  1  i, j 

where  is  to  be  taken  over  all  values  i  and  j  so  that  i  i  sind 

(0) 

3C>  is  the  energy  operator  for  an  individual  spin  in  a  constant  exter¬ 
nal  magnetic  field  ifg  .  It  will  now  be  assumed  that  all  spins  are 


I 


5 


identical  in  the  sense  that  their  spectra  are  identical  and  overlapping. 
In  Section  1.4  the  formulas  will  be  extended  to  the  case  when  there  are 


two  or  more  different  species  of  interacting  spins. 

The  perturbation  is  the  sum  of  the  dipolar  and  possible 

exchange  interactions  between  spins  i  and  j  . 

The  classical  interaction  potential  }C>j  for  two  dipoles 


j?.  and  jt.  joined  by  a  vector  .  is^^ 
^  J 


ill.  •  3(ir.  •?.)  (p.  •  ?. .) 

\  IJ  /  IJ 


(1-3) 


In  general,  the  g -value  of  a  spin  is  a  tensor,  so  p.  =  pg  S.  where  p 

is  the  Bohr  magneton  and  is  the  spin  operator.  The  matrix  elements 

of  K--  J-  are  found  by  expanding  it  in  terms  of  S_.  ,  S  +  .  *  S  .  +  S,  .  . 

ij  dip  '  *  "  zi  —  1  XI  —  yi 

As  an  example,  when  spins  i  and  j  have  g^^  along  the  direction 
of  iIq  (z-axis)  and  g^  in  the  x-y  plane,  we  have 


V..,.  =-fir(a  +  b+c  +  d+e  +  f) 
ij  dip  3 

A  •  • 

IJ 


(1-4) 


where 


a  =  g//S„,  S,.  (1  -  3  cos  0,.) 


zi  zj 


ij 


c  * 


•  '““ij 


d  »  +S_j  S^.  )  singj.  co.«jj  exp(i<j.) 


6 


e  =  S^.]  sin^  O.j  exp(-2i 

f  =  -|g^  [S_.  S_j]  sin^  0.j  exp(2i 

0..  is  the  polar  angle  that  makes  with  ,  and  .  is  the 
ij  ^  ®  ij  0  ij 

azimuthal  angle  of  . 

The  ordinary  exchange  interaction  is  caused  by  a  slight  overlap 
of  the  wavefunctions  from  different  spins  together  with  the  exclusion 
principle  that  prevents  two  electrons  from  being  in  the  same  state  at 
the  same  time.  Superexchange  is  interaction  via  other  ions,  usually 
oxygen,  lying  between  the  spins.  This  is  probably  the  source  of  the 
weak  coupling  in  most  of  the  hydrated  salts  we  have  investigated.  Since 
the  wavefunctions  in  the  crystals  are  at  best  only  crudely  known,  it  is 
not  possible  to  calculate  the  strength  of  the  exchange  from  first  prin¬ 
ciples.  It  has  to  be  found  from  experiments,  and  we  usually  take  it  to  be 
isotropic  , 

3C..  =  A.,  i!  =  A..  (S  .  S  .  +1  (S..  S  .  +  S  .  S.,)).  (1-5) 

^••ex  ^  J  ij  '  zi  zj  2  '  +1  -j  -I  i-y' 


Since  the  wavefunctions  fall  off  exponentially  vdth  distance,  it  is 
usually  assumed  that  this  coupling  is  important  only  between  nearest 
neighbor  spins.  A  positive  A.^  in  (1  -  5)  corresponds  to  anti- 
ferromagnetic  coupling,  a  negative  A.j  to  ferromagnetic  coupling. 

The  approximation  in  the  moment  calculation  is  to  treat  = 


(0) 


ij 


JCij  ^.p  +  3Cjj  gjj  a®  small  compared  to  the  Zeeman  energy  3q'  '  and 
keep  only  first-order  perturbation  terms,  that  is,  the  terms 
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that  commute  with  and  that  can  be  diagonalized  in  a  representa¬ 
tion  where  is  diagonal.  These  terms  will  slightly  shift  the 

energy  levels  a^  of  a  single  spin.  If  we  look  upon  the  N  spins,  each 
with  (2S  -f  1)  =  R  levels, as  a  coupled  system,  the  states  are  . .  .  n  , 

I 

.  .  . .  n  ,  ...  with  energies  ...£  ....Ei.  In  the  absence  of  the  per- 

“  n  n 

turbation,  a  large  number  N  of  states  has  the  same  energy  E  s 

Or  a 

N,a,  +  N,a-  +  ..N^a-  where  Ns  N  The  first-order  inter- 

11  2  2  R  R  r 

action  terms  couple  states  within  this  manifold  of  degenerate  states,  the 
degeneracy  is  lifted,and  the  energy  states  are  spread  out  slightly  around 

E  . 

a 

The  second-order  terms  in  would  weakly  couple  together 

different  levels  a^  .  .a^t  of  a  single  spin  and  permit  transitions  in  the 
microwave  field  of  Am  ^  f  1  .  For  a  system  of  spins  without  crystal¬ 
line  field  splitting,  the  terms  c  -  f  in  (1-4)  would  thus  cause 

weak  satellite  lines  at  hv  =  0  ,  2gpHQ  ,  Sg^Hg  in  addition  to  the  main 

line  at  hv  =  .  Since  their  separation  from  the  main  line  is  large, 

these  weak  lines  would  give  considerable  contribution  to  the  higher  mo¬ 
ments.  The  calculation  would,  therefore,  not  be  improved,  but  made  in¬ 
correct  by  including  the  second-order  terms  in  the  approximation,  and 
the  noncommuting  part  of  the  interaction  operator  has  to  be  thrown  away. 

The  properly  truncated  we  will  designate  and  we 

will  also  use  JC  =  .  It  should  be  remembered  that  if  a  crystal¬ 

line  field  of  low  symmetry  mixes  the  states  of  the  unperturbed  spin,  then 
may  also  contain  contributions  from  the  terms  c-f  in  (1-4). 


I 
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When  a  crystalline  field  is  present  and  the  spectrum  consists 
of  more  than  one  line,  the  microwave  field  operator  M  also  has  to  be 
truncated  to  M  to  include  only  the  matrix  elements  for  the  resonance 
transition  in  which  we  are  interested. 

The  power  P(v)dv  absorbed  by  the  spin  system  from  the 
microwave  field  of  frequency  v  is  proportional  to 


A  E  ^  . 

(b  "  -  b  "  )|  (n  |M|n')|^ 

I 

nn 


(1-6) 


where  ^  means  summation  over  all  states  |n)  and  |n  )  such 
nn' 

that 


hv  <  E  ,  -  E  <  h(v  +  dV) 
n 

and  the  Boltzmann  factors  are  written 


-E  /kT 
n 
e 


(1-7) 


The  line  shape  function  is  then 

P(v) 


^  6  E  E  l  /\  ,  y 

2,,  (b  "-b  ")|(n|M|n)|^ 


f(v)  = 


nn 


^  £  £  I 

f  P(v)dv  2  "')|(n|Mln')|^ 


(1-8) 


nn 


where  ^  means  summation  over  all  states  |n)  and  |n  )  such 


nn 

that  E  '  >  E  . 
n  —  n 


th 


From  the  definition  (1^1)  we  can  write  the  <t  moment  as 
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E  E 


X!"(E  ,  -  E  -  b“n')|(„|iJ|„’‘'2 


h^<v^>  =  22 


E  E  , 


(1-9) 


A  ,  I  ,2 

»  »  I  _  »  I  “ 


,b  ”  -  b  “  )  I  (n|M|n  ) 


To  avoid  the  summation 


operators 


Z.  • 


Kambe  and  Usui  introduce  the 


t,  A  I 

(n  Min)  for  E  i  <  E 

n  " 

0  otherwise 

r(n  |M|n  )  for  E  ,  >  E 

(n|M_|n)  =  <  ” 

0  otherwise 


.  (I- 10) 


Then  the  moments  can  be  written  in  invariant  trace  form 


h<  v>  = 


V  ^  ^ 

r  \.Ji  T  ftn  \M 


Trace  (b^[M_,  ]  ) 


(1-11) 


h^<V^>  = 


Trace  (b^[[M_,3C]  »  [Sf  ,M^  ]] 
Trace  (b^[  M_  ,  ]  ) 


(1-12) 


In  these  expressions  M.  and  3C  are  supposed  to  act  on  the 

▼ 

system  eigenfunctions  |n,k)  where  k  stands  for  the  de¬ 

generacy  of  the  state  with  energy  E^  .  But  the  states  |n,k)  are  made 
up  by  all  linear  combinations  of  the  eigenfunctions  of  a  single  spin,  and 
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we  can,  therefore,  evaluate  the  traces  with  respect  to  the  functions 
|a)j  |b)2  Ic)^  ...  Id)^  =  •••  d)  using 


Here 


|g)^  is  the  g 


ij 


eigenfunction  of  the 


.  (1-13) 

spin,  so  that 


if 


a 


a 


I  a,  b,  ...d)  =  I  a,  b . d) 

+  a.  +  a  +  . . .  +  a  ,  =  . 

DC  an 

The  density  matrix  can  be  approximated  by 


b 


K 


-  If  b 


a. 


(1-14) 


In  this  first  approximation  the  interaction  terms  are  neglected 
in  the  density  matrix  expansion.  This  is  only  valid  far  above  the  Curie 
temperature. 

With  the  substitutions  of  (l-i3)and(l-14),  the  traces  in  (1-11)  and 
(1-12)  become 


Trace  (b^  [M_  ,  ] 


2 

r  =1 


a  +  a  +  .  .  .  +  a  + 

r,  r,  r. 

12  1 


] 


=  (Trace  (b®))^"^  ^ 

i 


m+il  ) 


(1-15) 


I 
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Trace  (b^[M 

=  ss 

i  r 

[(3c/°^  +  I 

=  (Trace  (b^)  ) 
+  (Trace  (b^)  ) 


Trace  b^  [[M 


=  (Trace  (b*)) 
+  (Trace  (b*))^ 


.  [3C.  M]  ]  ) 


r.  •  “r.  ^  •  •  • 

1  J  >  r  /  ^  » 

Zj  l  (ni.  +  m.)  , 


1  y 


l  - 

r  3C  ij  )  ,  (m.  +  )]  ] 


N- 


^  ^  [*ni.[3c/°^  m.  ]  ]  ) 

i 

'*'^2  “i  ij*'*’  “i  +  “j]  1  > 


(1-16) 


.  5cl  .  [jc  ,  1] 


a  +  . .  a  + 
'i  ^• 


k  I 


Trace  (b*‘t  (Xjl®).{S^]  .  [x,<‘’>,  ft.j)  ] 


'■^E'  +  m*j)  I  . 

ij 


,  m_.]]  +  ([3Cj^°\ro^j  ]  ,  [3Cij^^^(m_.  +  TO  j)]] 
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+  (Trace  (b^)  ^  Trace  (b  ^ 

ijk 


a.  +  a.  +  a, 

1  j  k 


[  •  (^+i  +^+j^  ]•  •  <^-i  ■‘■“-k^  1  1  ) 


(1-17) 


These  formulas  were  derived  by  Kambe  and  Usui.  They  can 
be  applied  to  special  cases  by  substituting  the  proper  values  of 
and  . 


1.3.  Applications  of  tiie  moment  formulas  to  interaction  between  like 


neighbors  . 


A  simple  example  to  which  Kambe  and  Usui  applied  the  form¬ 
ulas  (1-11)  to  (1-17)  is  the  case  of  no  crystalline  field  splitting  of  the 
lowest  state.  Their  results  will  be  repeated  here  since  it  will  be  of 
use  later.  The  unperturbed  Hamiltonian  for  isotropic  g  is 


3Ci"”  =  gPHo®i 


(1-18) 


where  the  field  is  taken  to  be  along  the  z-direction.  The  trun¬ 
cated  interaction  will  contain  the  exchange  (1-5)  plus  the  terms 

a  and  b  of  the  dipolar  energy  (1-4)  .  This  we  write  as 


=  A..  S.  •  ^.  +  B..  S  .  •  S,. 
ij  ij  »  J  iJ 


(1-19) 


where 
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A..  + 
ij 


SLL. 


(3  COB 


0.. 

ij 


1) 


(1-20) 


B..  =  -  (3co8^e..  -1)  (1-21) 

2r?. 

ij 

The  oscillating  microwave  field  is  along  the  x-axis]  hence, 

m.  =  gps^. 

When  (1-18)  ,  (1-19)  and  (1-22)  are  inserted  in  (1-15),  (1-16),  and 
(1-17),  the  commutators  can  be  simplified  by  using  the  relations 


[S^.S.j]  .2«yS, 


(1-23) 


In  this  way,  Kambe  and  Usui  find  the  mean  deviation  from  the 

* 


unperturbed  transition  frequency  v  =  gP  ^Q/h 


h<A  v>  =  -  S  (y)2^  B.. 


(1-24) 


and  the  mean  square  deviation  about  v  =  V  4<Av> 


^  B.(y)^  By' 


(1-25) 


B_(y)  is  the  Brillouin  function 

S 


(1-26) 


14 


where  y  =  gpPi^S/kT  .  At  low  temperatures  >  1  for  T  — ^  0; 

so, depending  upon  the  lattice  sum,  we  may  have  a  finite  first  moment 

I 

shift.  The  sum  Jj  B*-  gives  the  effect  of  sample  magnetization 

j 

which  will  be  discussed  in  more  detail  in  Section  1.5.  Here  we  only  note 
that  for  a  cubic  lattice  and  spherically- shaped  sample  it  vanishes.  On 


the  other  hand,  dB  (y)/dy  — ^  0  for  T  — >  0  and  the  second  moment 

8 

vanishes  at  low  temperatures  as  it  should.  Isotropic  exchange  be¬ 


tween  like  spins  does  not  influence  the  first  or  the  second  moment. 

For  spin  S  =  -^  ,  the  g -tensor  in  many  cases  has  axial  sym  - 
metry.  If  the  magnetic  field  then  is  along  the  symmetry  axis,  the 
moments  are  already  given  by  (1-24)  and  (1-25)  if  we  use 


~  g,  P 

A..  + 

‘J  2r^ 

£r  .j 


(3  cos*'  0..  -  1) 
ij 


(l-20a) 


(2g/  +g^)P^  2 

B.j  - - S-5 - - -  (3  cos^  0  -  1)  (l-21a) 

^'ij 


Another  case  that  can  be  calculated  directly  and  which  applies 
to  the  case  of  Ni  SiF^  '  discussed  later  is  when 

S=1  =  gpS.  Sq  +D(sf.  -2/3)  .  (1-27) 

For  Hq  along  the  z-axis,  we  will  then  have  two  resonances,  and 

^nd  m^.  have  to  be  truncated  more  than  in  the  previous  case. 
Then  we  must  be  very  careful  in  using  the  commutator  relations  (1-23) 
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so  that  only  the  right  terms  are  included.  But  if  everything  is  written 
out  in  matrices,  the  multiplications  are  straightforward.  We  have 


e  +  •*•  D 


(0)  . 


0 

-Id 

0 


0 

0 


'C  +^D  / 


a3  0  0  \ 


0  a^  0 


0  0  aj  / 


(1-28) 


where  gpHf.=  e 


m 


+  i 


where  x  a  g  p  /T  H 


rf 


If  we  observe  the  transition  (a^  -  a^)  =  (e  -  D)  ,  we  have  to 


truncate  m  , .  to 
+1 


™+i  = 


0 

0 

0 


°\ 

X 

0 


(1-29) 


Similarly , 


A 

m 


-1 


0 

0 

X 


o\ 

0 


(1-30) 


By  nnatrix  multiplication,  we  find 
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/(A,,+B..)  0 


i 


ij  \y 
0 

0 

0 

0 

0 

0 

0 

0 


0 

0 


A.. 

IJ 


0  -(i\j+R,j)0 

0  0 


A.. 

ij 

0 

0 

0 

0 

0 


(0) 

0 

0 

0 

0 


0 

0 

0 

0 

0 


0 

0 

(0) 

0 

0 

0 

(0) 

0 

0 


0 

0 

0 

0 

0 

0 

0 


A.. 

IJ 


0 

0 

0 

0 

(0) 

0 


0 

0 

0 

0 

0 


A.. 

IJ 


-(A.j+B.j)  0 


0 

0 


0 

0 

0 

0 

0 

0 

0 

0 


A 


0  A,.+B^/(l-34) 


where  the  positions  (0)  show  Hie  positions  of  truncated  spin-flip 
terms. 

(m^  -I*  mj)  can  be  written  in  the  same  way,  and  we  find  for  the 
diagonal  terms  of  [m.,  ,  (m.  +  “^j)]  1 


(1-35) 
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Combining  (1-31),  (1-32)  and  (1-35)  with  (1-11),  (1-15)  and  (1-16), 
we  find  the  first  moment  for  the  transition  (a^  a^) 


h  <  Av  >  =  / 


vn  '  [(A..  +  B..)  b  ^  -  A.,  b  ^  -  B..  b  ^  ] 
\  '■'ij  ”  ■' 


ij 


^1  ^2  *3 

(b  ^  +  b  ^  +  b  "*  ) 


(1-36) 


The  trigonal  lattice  of  NiSiF^*  6H2O  is  sufficiently  close  to 
cubic  so  that  for  a  spherical  sample  the  dipolar  part  of  the  interaction 
in  (1-36)  cancels.  The  exchange  part  of  A.j  ,  however,  will  give  a 
characteristic  temperature  dependent  shift  in  the  first  moment,  and 
from  the  observed  shift  we  can  find  the  magnitude  and  sign  of  A.j  . 

The  first  moment  for  the  transition  (a^  >  a^)  can  be  calcu¬ 

lated  the  same  way.  But  to  find  the  second  moment,  we  would  have 

3 

to  work  with  (2S  <f  1)  =  27  -  fold  matrices,  and  while  this  should  not 

be  impossible,  it  is  an  unnecessarily  hard  way  of  doing  the  calculation. 
McMillan  and  Opechowski  have  calculated  the  commutators  in  (1-11) 
and  (1-12)  in  a  general  way  using  the  technique  of  projection  operators. 
Their  formulas  for  the  first  and  second  moment  are  given  in  Appendix 
1  .  These  formulas  only  involve  summations  over  the  truncated  ma¬ 
trix  elements  connecting  the  unperturbed  energy  states. 

1.4.  Line  moments  from  interaction  between  \mlike  spins. 

The  previous  moment  calculation  can  be  extended  to  include 
the  effect  of  unlike  spins.  The  importance  of  differentiating  between 
"like"  and  "unlike"  neighbors  comes  from  the  fact  that  for 


1 
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interaction  between  unlike  spins  a  further  trtincation  of  is  neces¬ 

sary,  because  then  terms  like  S^.  S  .  do  not  connect  states  of  the 
same  energy. 

The  criterion  for  "unlikeness"  is  that  the  spins  have  iinper- 
turbed  transitions  sufficiently  different  in  energy  so  that  the  resonance 
lines  have  little  or  no  overlap.  It  is  obvious  that  different  chemical 
species  of  ions  may  have  different  magnetic  properties.  But  also 
chemically  identical  ions  may  be  unlike  in  this  sense,  for  instance,  if 
the  crystalline  field  is  different  in  strength  or  direction  on  different 
crystalline  sites,  ox  if  identical  spins  interact  so  strongly  with  a  par¬ 
ticular  neighbor  or  the  nuclear  spin  that  separate  lines  are  created. 
Whether  the  interaction  is  like  or  unlike  depends  to  some  extent  on 
the  strength  of  spin  interaction;  strong  interaction  will  tend  to  make 
separate  lines  overlap. 

The  Hamiltonian  for  two  kinds  of  interacting  spins,  primed  and 
unprimed,  is 


The  density  matrix  now  has  to  be  approximated  by 
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(1-38) 


We  may  suppose  that  we  are  looking  at  a  resonance  of  the  un¬ 
primed  system. 

When  now  the  traces  in  (1-11)  and  (1-12)  are  expanded  the  same 
way  as  before,  we  find  that  the  additional  terms  in  the  moments  due  to 
interaction  with  unlike  neighbors  are  : 


^  ‘  ’•unlike 

Trace  {b  ‘  ''  [in.  ,  (jc  ,***  ,  m.  ]  ]  ) 


1  k 

I 


ik 


k*  A  ^ 

Trace  b  Trace  (b  ^  t  ]  ) 


^  ,  2 
h  <  A  V  > 


unlike 


a.  +  a 


=  2  Trace  (b  ^  »  ^+il  '  ™-il  1 


(1-39) 


+  [  .  m-i  ]  .  [X  ,  m  .  ]  ] 

^  ik  ‘ 

+  [  [3C  .  m.j  1  A  ) 

ik  ”  ik  ' 

r  *^1.'  ^  ^ 

sTrace  b  *  Trace  (b  ^  ]  )  r 


-1 


I 
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+  2  Trace  (b 


‘  ■'  *  (  (5c  . 

ik  ” 


[5c  .  A  j]  ]  ) 

ie 


a  ,  +a  , 

Trace  (b  ^  ®  ) 


rA  A 


Trace  (b  [m_.,  m^.  ]  ) 


a.  +  a.  +  a  , 
Trace  (b  ^  J  ^ 


[  [K  .  m  .  ]  , 


[3C  ,  (m^.  +  m  .)  ]  ] 

ij  J 


+  [  [3C  .  (m  .  +  m  .)  ]  ,  [jc  ,  rn  .  ]  ]  ►  ) 

ij  ik  ‘ 


Trace  b 


a  ,  a 
k 


Trace  b  ^ 


(1-40) 


We  see  that  for  more  than  one  kind  of  different  neighbors  the 
contribution  to  h  <  ^  v  >  can  simply  be  added  for  each  specie.  The 
last  terms  in  h  <  V  >„„nvp(  cross  terms  between  like  and  \m- 
like  neighbors.  But  such  cross  terms  can  be  expected  to  cancel  out 
when  we  calculate  the  more  useful  central  second  moment  of  a  line. 

Equations  (1-39)  and  (1-40)  will  now  be  applied  to  a  special 
case  :  Two  magnetic  species,  both  of  which  have  spin  S  -  and 
g -Censors  widi  symmetry  axis  along  the  direction  of  the  magnetic 
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field.  The  results  will  be  used  later  to  explain  the  measurements  on 
neodymium  ethylsulfate. 

Using 


(1-41) 


=  A.j  1  <s.i  s^j  "  S-i'  +  <Ay  +  ®ij>  Sxj  '  *-«> 


jc  ,<*>-  (A  ,  +  B  ,)S  jS  • 
ik  ik  Ik 


where 


\  =  *ic  *  2  V'"  ^  “ic  ■ 

'ic 


B.,  *  -  -y  - g -  (3  cos^  0._  -1) 

^  3 

'iC 

I 

and  (  stands  for  j  or  k  ,  the  first  moment  shift  becomes 


h  <  iv  >  =  h  <  iv  ♦  h  <  4V 


<  s  .  > 

zi 


2 


By  +  <S 

Ik 


I 

,>y  (A  ,  +B  ,) 

ik  ik 


k 


(1-43) 
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Here 


<S 


H 


1  (1  -  b 

>  =  -I  ^ - 


®//C 


) 


8//’^ 
{l+b''c 


(1-44) 


The  second  central  moment  becomes,  in  this  case, 


(<  s  <s .  >^)  y  B./ 

'  Z1  Zl  '  IJ 


+  (<S 

zk 


<S  ,>^)y  (A  ,+B  ,)^  (1-45) 

zk  ik  Ik 


where 


<S^^>  .  4  S^(Sj  +  l)  =  I 


(1-46) 


Here  we  should  note  that  B. ^  and  the  dipolar  part  of  are 
defined  with  opposite  signs  so  that  they  partially  cancel  when  they  ap¬ 
pear  together  in  a  sum.  Then  we  see,  as  is  well  known, ^  the  unlike 
neighbors  are  less  efficient  in  dipolar  broadeningof  a  resonance  line. 
But  exchange  between  unlike  neighbors  will  contribute  to  the  second 
moment.  We  should  also  note  that  at  low  temperature  where 

<  S  ,  >  is  finite,  we  will  have  an  exchange  induced  shift  of  the  first 
zk 

moment. 


I 
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1.  5.  Discussion  of  demasnetizing  effects. 


We  have  seen  that  in  the  calculation  of  line  monnents  we  have 


to  evaluate  sums  over  the  lattice  of  the  form 


j 

Such  a  sum  is  only  slowly  converging.  But  following  the  classical 

13 

theory  of  demagnetization,  the  summation  can  be  performed  by  di¬ 
viding  it  into  two  parts.  Close  to  the  spin  i  expression  (1-47)  has  to 
be  summed  over  the  lattice  points.  But  outside  a  radius  Q  from  i 
we  can  treat  the  magnetization  as  being  continuous  rather  than  asr* 
sociated  with  discrete  lattice  points,  and  this  part  of  the  sum  can  be 
converted  to  an  integral.  This  integral  will  depend  on  the  shape  of  the 
crystal;  and  it  will,  in  general,  be  different  for  different  lattice  points. 
Only  for  an  ellipsoid  will  the  integral  be  the  same  for  all  points  i  ,  and 
only  for  a  sphere  will  it  vanish.  For  a  cubic  lattice  the  sum  inside  Q 
also  vanishes. 

When  the  field  is  along  one  of  the  ellipsoid  axes,  the  shape 
dependent  part  of  the  sum  (1-47)  can  be  described  by  a  demagnetizing 
factor  N 

z 

E 

j" 

■'ij  >C 

where  is  the  density  of  spins. 


3(cos  Oj.  -  1) 


=■  N,  (N.-  ^  ) 


(1-48) 
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The  factor  equals  4ir/3  for  a  sphere.  For  ellipsoids 


14 


of  revolution  around  the  z-axis,  is  given  by 


a)  Prolate  spheroid 


r  2 

(m'  -  X)  l2(m^  - 1)  '  IT)  -  (n>  -  1)  -* 

where  tha  axis  ratio  m  =  ’^z  /  rx 
b)  Oblate  spheroid  : 

1/2 

N  •  ft  -  —i—  .rc  .in  1 


(l-«) 


(1-50) 


where  m  =  ^x  ^r^ 

From  (1-24)  we  now  get  the  change  in  first  moment  for  free 
spins  in  a  cubic  lattice  : 


<  Av  >  «  -  I  ^ 

where  the  magnetisation  is  given  by 

<  >  =  Nj  gp  <  S^>  =  -  Nj  gp  S  »  -  M  .  (1-52) 

Here  Nj  is  the  number  of  spins  per  unit  volume  and  is  flie  Bril- 
louin  function  (1-26)  . 

It  is  interesting  to  compare  the  first  moment  shift  (1-51)  to  the 

15  ■ 

Kittel  formula  for  ferromagnetic  resonance 
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vh  =  gp  ^[HQ  +  (N^-N^)M][HQ  +  (Ny-N^)M]|  (1-53) 

Here  N  and  N  are  the  demagnetization  factors  in  the  x-  and  y- 
X  y 

directions.  They  are  related  to  N  by 

2 

N  +  N  +  N  =  4ff  (1-54) 

X  y  z 

For  axial  symmetry,  so  that  =  •^(4ir  -  N^),  the  resonance 

condition  (1-53)  becomes 

vh  =  gP(HQ  +  (N^  -  N^)  M) 

Comparing  (1-55)  to  (1-51),  we  see  that  the  shape  dependent  lineshifts 
depend  in  the  same  way  on  the  magnetization  in  these  two  cases  as  we 
expected.  For  a  general  ellipsoid  with  N  i  N  ^  N  ,  we  see  that  the 
first  moment  shift  (1-51)  differs  in  order  M^/Hq  from  the  resonance  con¬ 
dition  (1-53).  This  is  because  the  component  of  the  magnetization  in  the 
direction  of  the  applied  field  is  not  a  constant  of  the  motion  in  the  case  of 
a  general  ellipsoid,  as  has  been  pointed  out  by  Van  Vleck.  The  moment 
method  is  not  strictly  valid  unless  the  sample  has  cylindrical  symmetry 
around  Hq  .  Only  in  this  case  can  we  truncate  the  interaction  Hamiltonian 
in  the  representation  where  is  constant. 

The  appearance  of  the  factor  that  signifies  a  dynamic  de¬ 
magnetizing  effect  in  the  formula  for  the  first  moment  change,  may  be 
thought  of  as  being  caused  by  the  spin-flip  terms  S.j  in  the 


interaction.  Suppose  these  terms  were  trxincated  out  (in  the  case  of 
no  exchange)  and  we  would  get  only  the  static  demagnetization 

h  <  Av>  =  -  gp  (^  -  N  )  <  M,  >  (1-56) 

It  is,  therefore,  not  surprising  that  in  the  interaction  between 
nonequal  neighbors  or  when  crystalline  field  splitting  is  present  and 
further  tmncation  is  necessary,  the  Kittel  resonance  condition  is  not 
immediately  valid  in  the  sense  that  in  (1-55)  M  cannot  simply  be  re¬ 
placed  by  the  static  magnetization  -  <  >  .  In  the  case  of 

NiSiF^'  6H2O  ,  for  example,  the  magnetization  relevant  to  the  shape 
dependent  first  moment  shift  of  the  lower  line  is  found  from  (1-36)  to  be 


(1-57) 


(1-58) 


The  formulas  (1-57)  and  (1-58)  should  replace  <  M.  >  in 

z 

(1-51).  This  correction  for  elliptically- shaped  samples  is  later  ap¬ 
plied  to  our  measurements  on  NiSiF^  *  . 

If  the  crystal  lattice  is  not  cubic,  we  also  have  to  consider 
the  part  of  the  sum  (1-47)  taken  over  r.j  <  Q  which  may  give  an 
additional  line  shift. 
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The  second  central  moment  formula  involves  only  the 
rapidly -converging  sum 

2  ^  (3  cos^  -  1)  . 

j  'ij 

and  the  second  central  momentwill  not  depend  on  the  sample  shape  as 
long  as  it  is  spheroidal. 

When  a  sample  deviates  from  exact  ellipsoidal  shape,  the  de¬ 
magnetizing  field  will  vary  from  spin  to  spin.  We  may  expect  this  to 
have  a  relatively  larger  effect  on  the  second  moment  than  on  the  first 
moment,  and  in  particular,  we  may  not  see  the  narrowing  of  the  reso¬ 
nance  lines  at  low  temperature. 

1.6.  Comment  on  the  fourth  moment  and  the  effect  of  exchange. 

The  method  of  calculating  moments  can,  in  principle,  be  ex¬ 
tended  to  the  n^  moment,  and  more  information  about  the  line  shape 
may  be  gained.  Van  Vleck^  has  calculated  the  fourth  moment  at  high 
temperature  for  the  case  of  no  crystalline  field  splitting.  He  found: 

a)  Isotropic  exchange  with  like  spins  will  not  change  the  second 
moment,  but  it  will  increase  the  fourth  moment.  This  can  only  happen  if 
the  resonance  line  is  being  so-called  exchange -narrowed  --  made  more 
narrow  and  sharp  in  the  center  and  with  longer  wings  on  the  sides. 

b)  If  two  unlike  systems  of  spins  are  interacting,  isotropic 
exchange  between  the  spins  in  one  system  will  exchange -narrow  the 
resonance  line  of  the  other  system. 
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The  fourth  moment  at  low  temperature  can  be  written  as 

.  .  Trace  (b^  [[3C  .  [5C,  M  ]],[5c,[K.M  ]]]) 

h^<v^>  =  - ;; - : - - i -  .  (1-59) 

Trace  (b^  [  M_  ,  ]  ) 


The  traces  can  again  be  calculated  as  in  Section  1.  2  by  in¬ 


serting 


and  taking 


ij 


b*'  -TTb** 

i 


as  the  approximation  for  the  density  matrix.  This  will  lead  to  even 
more  complicated  formulas  than  those  given  by  Van  Vleck.  At  low 
temperature  where  Trace  (b  S  .)  and  Trace  (b  S.,  )  are  no  longer 

SI  Z1 

xero,  many  terms  involving  these  factors  will  appear. 

These  terms  require  that  we  perform  awkward  summations 

t 

over  the  lattice  like  A  A..  A..  B..  B..  even  before  we  can  tell 

^  ij  jk  ij  ik 

whether  they  are  positive  or  negative  or  in  other  words  whether  they 
add  to  or  subtract  from  the  fourth  moment  at  low  temperatures.  Be¬ 
cause  of  the  labor  necessary  and  the  limited  information  to  be  gained, 
we  will  not  here  try  to  calculate  the  fourth  moment. 

A  more  direct  approach  to  the  problem  of  resonance  lineshapes 

6  7 

and  exchange-narrowing  has  been  attempted  by  Anderson  and  Weiss.  ’ 


t 
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They  consider  the|  exchange  interaction  to  give  a  random  frequency 

modulation  to  the  dipolar -broadened  resonance,  similar  to  the  motional 

modulation  of  nuclear  spin- spin  interaction  in  liquids  treated  by 

12 

Bloembergen,  Purcell  and  Pound.  The  line  shape  can  then  be  found 
from  the  Fourier  transform  of  an  autocorrelation  function 


^(t)=<p..  (t  +  T)  |i..  (t)  > 

ij  ij  average  over  t 


(1-60) 


^(T)  is  a  measure  of  the  "memory"  of  the  state  p..  at  the 
time  t  T  for  what  it  was  at  a  time  t  earlier.  Strong  exchange  will 
cause  rapid  spin  flips  that  reduce  the  correlation  time  and  average 
out  the  dipolar  field  fluctuations,  and  the  line  will  appear  exchange - 
narrowed. 

Assuming  Gaussian  random  modulation,  they  find  that  the  cen¬ 
tral  part  of  an  exchange -narrowed  line  must  have  Liorentzian  shape. 
But  the  intensity  of  the  line  wings  will  decrease  more  rapidly,  as 
expected,  since  a  strictly  Lorentzian  line  would  have  an  infinite  se¬ 
cond  moment.  Qualitatively,  in  the  case  of  strong  exchange,  the  line- 
width  is 

(1-61) 

where  is  the  effective  dipolar  field  and  is  the  effective  field 
from  exchange  interactions.  For  more  quantitative  results,  the  line 
parameters  would  have  to  be  fitted  so  that  they  give  the  same  second 
and  fourth  moment  as  given  by  Van  Vleck's  theory. 
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In  this  model,  Anderson  also  treats  the  case  of  two  lines  close 
together  that  merge  into  one  as  the  exchange  is  increased.  This  is 
beyond  the  power  of  the  moment  method.  From  Fig.  2  in  his  article, 
it  is  clearly  seen  that  increasing  exchange  between  the  unlike  spins 
first  broadens  the  two  lines  so  that  they  overlap  more  and  more,  and  at 
the  same  time  the  peaks  move  together.  This  behavior  is  observed  in 
the  case  of  CuSO^  ■  5H2O  to  be  described  later. 

For  many  lines  close  together,  Anderson's  qualitative  result 
is  that  the  exchange  will  broaden  them  and  wash  out  their  identity,  but 
while  still  separated ,  the  shift  of  each  peak  is  negligible.  If  there  are 
p  lines  separated  by  hAv  ,  we  expect  qualitatively  that  an  exchange 
3Cg'  »hAv  will  make  exchange  spin  flips  possible  between  neighboring 
lines  and,  therefore,  broaden  the  separate  lines  into  one.  But  a  con¬ 
siderably  stronger  exchange  'I  s(p-l)hAV  is  required  to  exchange 
narrow  appreciably  the  composite  line  .  .  This  seems  to  be  the  ob¬ 
served  behavior  of  the  hyperfine  lines  of  Mn  'wdien  the  exchange  is 
changed  by  varying  the  concentration.  For  an  exchange  3c  in  the 
approximate  range  hAv  <  }C-<  (p-1)  h  AV  ,  we  have  to  treat  the  inter- 
action  as  partly  between  like  and  partly  between  imlike  neijghbors.  Ex¬ 
change  between  unlike  neighbors,  however,  will  lead  to  a  shift  of  the 
resonance  frequency  at  low  temperatures,  and  we  have  found  it  neces¬ 
sary  to  use  this  model  to  explain  the  line  shift  in  K2Cu(SO^)2*  6H2O 
where  the  Cu  resonance  has  p  =  4  hyperfine  lines. 
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8  9  10 

Kubo  and  Tomita  '  ’  have  given  the  ideas  of  Anderson  and 
Weiss  a  more  rigorous  foundation  by  starting  from  the  theory  of  ir¬ 
reversible  processes  in  quantum-mechanical  systems.  Their  treat¬ 
ment  includes  the  Boltzmann  factors  necessary  if  we  go  to  temperal- 
tures  where  hv^kT  .  The  theory  is  essentially  equivalent  to  the 
moment  method  :  each  higher  order  in  the  expansion  of  the  relaxation 
function  corresponding  to  ^(T)  in  (1-60)  gives  only  the  same  informa¬ 
tion  as  the  corresponding  line  moment,  and  requires  about  the  same 
amount  of  work.  Only  for  partly -over lapping  lines  or  ^en  the  inter¬ 
action  fields  are  comparable  in  strength  to  the  external  field,  does 
this  approach  seem  to  be  superior. 

It  might  be  pointed  out  that  in  the  latter  case  when  the  exchange 
is  strong,  all  terms  in  the  dipolar  interaction  (1~4)  contribute  to  the 
linewidth  and  should  not  be  truncated  out.  ^  Physically,  we  can  see 
that  dipolar  rearrangements  that  do  not  conserve  energy  become  pos¬ 
sible  when  the  fluctuations  in  exchange  energy  are  very  great.  Thus, 
the  complete  second  moment  for  a  highly  exchange -narrowed  line  will 
be  10/3  times  its  value  calculated  from  the  truncated  dipolar  inter¬ 
action. 

1.7.  Comment  on  McMillan  and  Opechowski*s  second  approximation. 

The  approximation  for  the  density  matrix  (1-14)  can  be  made 
more  accurate  by  using  the  fact  that  the  interactions  are  weak,  and  the 
part  of  the  Boltzmann  factor  containing  these  interactions  can  then 
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be  expanded  in  powers  of  (kT)"^  .  Such  higher  approximations  in  the 

derivation  of  the  line  moments  have  been  considered  by  McMillan  and 
5  18 

Opechowski.  ’  Using  the  projection  operator  formalism  explained 
in  Appendix  I  ,  they  write 

P 

bol  (P^  3C  P^)  =  bol(E^  P^)[  1  -  - i  +  . . .  ]  .  ( 1-62) 

In  the  first  approximation,  only  the  first  term  in  the  bracket  is  kept, 

and  this  is  the  approximation  we  have  worked  with  so  far.  We  write 

the  resulting  moments  as  h*'<Av  .  If  the  term  P  X^^^P  /kT 

1  a  a 

is  also  kept  in  (1-62),  this  leads  to  the  second  approximation  moments 
h^  <  A  >2  •  The  labor  involved  in  calculating  higher  approximations 
increases  rapidly,  and  McMillan  has  given  a  formula  only  for  h<^v>2* 
His  general  equation  is  copied  in  Appendix  II  . 

Because  of  the  extra  term  in  (1-62),  we  expect  the  first 
moment  in  the  second  approximation  to  contain  terms  in  (X^^^)^  . 

Even  for  a  cubic  lattice,  this  will  give  a  finite  shift,  which  is  related 
to  the  magnitude  of  the  second  moment  in  the  first  approximation. 
McMillan  and  Opechowski  show  quite  generally  that 

V*  <  Av  >j  +  <  A  >1 

<Av>2  - >  - i - - - - 

T-^£X>  V  +  <  A  V  >1 

<  Av^>j 

for  <  AV  >,  - >  0 

T-400 


8 


V 


.  (1-63) 
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Here  V  is  the  unperturbed  frequency  of  the  resonance.  The 
odd  moments  of  the  line  generally  vanish  at  high  temperature  in  the 
first  approximation. 

Equation  (1-63)  can  easily  be  derived.  The  spin-spin  inter  - 
action  causes  a  variation  in  local  field,  and  the  possible  transition 
frequencies  will  have  a  distribution  S(v)  .  The  intensity  of  the  ab¬ 
sorption  at  each  frequency  V  is  proportional  to  the  population  dif- 
ference  (b"^  -  b  '  between  the  initial  and  final  energy  level. 

When  the  Boltzmann  factors  are  expanded,  the  first  moment  in  the 
second  approximation  is  by  definition 


<4V>^. 


(1-64) 


We  assume  that  the  distribution  S(V)  can  be  calculated  from 
the  spin-spin  interaction,  and  [S(V)]^  calculated  will,  in  general,  be 
different  for  each  approximation  n  in  3C..^^VkT  .  For  T  — >  oo  , 
however,  the  population  difference  between  energy  levels  goes  to  zero 
and  in  this  limit  [S(v)^,j  ^  S(v)  •  From  (1-64),  we  then  get 


[  x[S(v*  +  x)]j(v*  +  x)  dx 

<  AV  >2  - >  ^ - 

T-»oo  y*  jg  (V*  +x)]j(  V*  +x)  dx 


(1-65) 
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This  directly  gives  (1-63)  ,  since  we  also  have,  by  definition 


00 

j  X*”  [S  (V*  +x)]j  dx 

<  Av*”  s  — — —  —  -  '  (1-66) 

f  [S(V*+x)3i  dx 
-00 

The  correction  (1-63)  is  small  and  we  have  never  seen  it  ap¬ 
plied  to  measurements  of  g -values.  It  may  be  of  importance  when 
broad  lines  are  measured  at  low  frequency.  If  a  resonance  with 
second  moment  of  (500  gauss)^  ,  say,  is  measured  in  a  field  of  2500 
gauss,  the  relative  error  in  g  ,  if  not  corrected  by  (1-63)  will  be 


M  -  -L 

g  ■  H. 


<  A 


*  4  i 


(1-67) 


McMillan  and  Opechowski  have  applied  their  general  eqiiation 
for  <Av  >2  to  the  transition  in  NiSiF^  *  6H2O  .  The  difference 
from  the  first  approximation  is  small  at  all  temperatures,  and  we 
will  neglect  higher  approximation  corrections  in  our  discussion  of 
experimental  data, 

1.  8.  Comment  on  the  possibility  of  magnetostatic  modes  in  para  - 
magnetic  crystals  at  low  temperature. 

It  is  well  known  that  under  certain  circumstances  ferromagnetic 
resonance  experiments  can  be  complicated  by  the  appearance  of  mag¬ 
netostatic  modes  in  the  sample,  propagation  effects  in  large  samples 


I 


36 


and  spin  wave  instabilities  at  high  microwave  power  levels.  Para¬ 
magnetic  salts  become  magnetized  and  similar  to  ferromagnets  at  low 
temperature  when  the  spins  are  aligned  by  an  external  field.  We  may 
then  ask  whether  we  should  see  some  of  these  effects  in  our  line  shape 

experiments.  In  particular,  magnetostatic  modes  have  recently  been 

19 

observed  in  disks  of  K2CuCl^  ■  ZH^O,  and  we  have  also  seen  some 
extra  lines  in  the  spectra  of  CuSO^  *  and  K2Cu(SO^)2*  6H2O 

below  l^K  that  probably  are  caused  by  such  modes.  Here  we  will 
briefly  review  the  theory  for  the  various  effects  observed  in  ferromag¬ 
netic  resonance.  The  experimental  evidence  for  nnagnetostatic  modes 
will  be  discussed  in  Section  5.  5  . 

If  the  magnetizations  in  different  parts  of  the  sample  precess 
with  the  same  frequency  but  in  different  phases,  the  regions  will  inter¬ 
act  via  the  dipolar  fields,  and  the  resonance  condition  may  be  changed 
relative  to  that  of  the  uniform  mode.  Walker^^'^^  has  given  the  reso¬ 
nance  frequencies  for  a  number  of  these  magnetostatic  modes  by 
solving  the  equations 

dil?/dt  =  Y  (M  X  H  )  (1-68) 


curl  H, 


in 


=  0  div(Hj  +  4Tr  Mj  )  =  0  (1-69) 

the  limit  of  small  deviations  from  the  equilibrium  magnetization 
So  and  with  the  proper  boundary  conditions  for  a  spheroid.  The 
energy  differences  between  the  modes  are  proportional  to  Mq  and  for 
a  sphere,  they  fall  over  a  frequency  range  of  Y  ■  4ir  Mg  . 
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Fletchev  et  al.  and  others  have  obtained  good  agreennent  between 
Walker's  theory  and  measured  mode  positions. 

Magnetostatic  modes  are  a  result  of  dipolar  interaction  only. 

The  exchange  forces  are  of  short  range  compared  to  the  mode  wave¬ 
length  and  they  are  neglected  in  Walker's  theory.  His  results  should, 
therefore,  be  directly  applicable  to  paramagnetic  resonance  at  low 
temperature. 

Whether  a  particular  mode  is  excited  or  not  depends  on  the 
symmetry  of  the  microwave  field  at  the  sample.  In  a  perfectly  homo¬ 
geneous  driving  field,  only  the  fundamental  mode  would  be  excited.  It 
23 

is  found  that  \mless  special  care  is  taken  in  mounting  the  sample 
away  from  field-disturbing  dielectric  materials,  one  is  likely  to  excite 
a  whole  series  of  modes,  and  they  may  partly  be  so  close  together  in 
frequency  that  they  are  not  observed  as  separate  lines. 

23 

Consider  now  the  modes  in  a  non -spherical  sample.  Dillon 
found  that  small  deviation  from  sphericity,  such  as  results  from 
tumble -grinding  an  anisotropic  material,  would  slightly  change  the  sepa¬ 
ration  of  the  modes.  Large  deviations,  however,  will  cause  the  de¬ 
magnetizing  field  to  vary  through  the  sample,  the  magnetization  in  dif¬ 
ferent  regions  will  process  with  different  frequencies,  and  there  will 
be  no  discrete  modes  for  the  crystal  as  a  ediole.  It  is  hard  to  give  a 
good  criterion  for  how  close  to  spheroidal  shape  the  sample  has  to  be 
before  separate  magnetostatic  modes  can  be  observed,  but  it  may  well 
be  possible  that  a  sample  is  good  enougl^  for  measurements  of  first 
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moment  lineahifts  and  still  have  sufficiently  crude  shape  to  prevent 
observation  of  discrete  modes.  Since  we  have  observed  these  modes 
only  in  rather  good  spheres,  it  may  even  be  advantageous  to  use 
crudely- shaped  samples. 

It  should  be  noted  that  homogeneous  broadening  of  the  line  by 
dipolar  and  hyperfine  interactions  will  not  prevent  the  possible  exci¬ 
tation  of  magnetostatic  modes.  The  macroscopic  magnetization  of 
one  part  of  the  sample  will  still  precess  on  the  average  with  the  same 
frequency  as  that  of  another  part.  However,  if  the  paramagnetic  salt 
has  two  or  more  separate  resonances,  the  microwave  field  will  usually 
excite  only  one  of  them  in  a  given  external  field.  The  energy  dif¬ 
ference  between  magnetostatic  modes  will  then  be  given  by  the  mag¬ 
netization  of  that  sub-lattice  alone.  The  spectrum  is  likely  to  be 
complicated  when  the  possible  mode  patterns  of  two  separate  resonances 
overlap. 

We  expect  any  single  mode  to  have  the  same  linewidth  as  the 
uniform  mode.  If  higher  modes  are  excited,  line  areas  can  no  longer 
be  used  to  measure  resonance  intensity,  since  each  spin  may  resonate 
in  a  succession  of  modes  as  we  sweep  through  the  line  . 

If  the  sample  is  made  so  big  that  it  becomes  comparable  to  the 

driving-field  wavelength,  the  finite  propagation  speed  of  the  dipolar 

fields  has  to  be  accounted  for  by  using  the  complete  Maxwell's  equa- 

24 

tions  instead  of  the  approximations  (1-69).  This  will  change  the 
resonance  frequencies  of  the  first  few  modes  of  long  wavel^Bngth  and. 
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in  particular,  the  resonance  condition  for  the  uniform  mode  in  a 
sphere  becomes^® 


Sir  y 

0)  =  Y  Hq  ^  1  -  (I  +  ^  )  (kEg)  + 


where 

c  =  dielectric  constant 
k  3  2  ir/X  propagation  constant 

Rq  =  sample  radius 

The  magnetization  at  T  =  0^  for  NiSiF^  *  6H2O,  for 
example,  is  from  (1-52) 


(1-70) 


Mq  =  NjgpS  =  86  gauss 

for  N|^  =4.16  *  10^^  ions/cm^  ,  g  s  2.  24  and  S  3  i  .  If  we  take  the 
dielectric  constant  for  this  salt  to  be  c  w  2.  5  ,  work  at  30  kMc/s  so 
that  k  =  2ir  cm'^  and  use  the  sample  size  Rq  «  0.075  cm  (wave¬ 
guide  spectrometer),  we  find  the  size  induced  resonance  shift  for 
constant  w  at  T  3  0°K  to  be 

A  H  »  —5—^  (1  +  V^qI  =  80  gauss. 

Our  cavity  spectrometer  in  which  most  of  the  escperiments  on 
NiSiF^  *  6H2O  were  done  requires  a  sample  size  of  about  Rq  » 

0.015  cm.  Then  the  shift  becomes  6H^4  gauss,  udiich  is  negligible. 
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In  Nd{C2H5S04)3  •  9H2O  and  K2Cu(SO^)2  *  6H2O  the 
saturation  magnetization  is  smaller,  and  we  have  ignored  the  size  ef¬ 
fect  even  for  measurements  in  the  waveguide. 

The  higher  order  magnetostatic  modes  can  be  thought  of  as 
going  gradually  over  into  traveling  spin  waves  obeying  the  dispersion 
relation^^ 

1 

uKk)  =  y[  (Hq  -  N^Mq  +  (Hq  -  +  4it  sin^O]^  . 


(1-71) 

Here,  is  the  demagnetization  factor  in  the  direction  , 

d  is  the  lattice  spacing  and  0  is  the  angle  between  *1^  and  the  spin- 

wave  propagation  vector  k  .  The  term  is  important  in  fer  - 

romagnets  where  the  exchange  field  is  strong,  but  it  is  negligible  in 

paramagnets  magnetized  in  strong  fields  at  low  temperatures. 

The  various  magnetostatic  modes  and  spin  waves  are  non- 

interacting  in  the  first  approximation,  but  they  are  weakly  coupled  by 

the  dipolar  fields,  surface  irregularities  and  lattice  defects.  The  way 

the  fundamental  mode  is  coupled  to  high-order  spin  waves  is  very  im- 

27  28 

portant  for  the  spin-lattice  relaxation  process  in  fer  romagnets,  ’ 
and  the  coupling  may  also  cause  instabilities  of  the  resonance  at  high 
power  levels.  ’  However,  we  will  not  here  speculate  on  the  pos¬ 
sibility  of  such  spin  wave  effects  in  paramagnets  at  T  w  0  since  we 
observed  no  instabilities  and  did  not  try  to  measure  relaxation. 
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CHAPTER  2 


DISCUSSION  OF  EXPERIMENTAL  APPARATUS  AND  TECHNIQUE 

2. 1.  Description  of  microwave  and  low-temperature  systems. 

In  the  line  moment  theory  reviewed  in  the  previous  chapter,  all 

-hv/kT 

shifts  are  functions  of  the  Boltzmann  factor  e  .  Resonance  ex¬ 

periments  have  usually  been  done  in  magnetically  diluted  salts  under  con¬ 
ditions  where  hv«kT  ,  and  no  observations  of  the  temperature  effect 
seem  to  have  been  reported  before.  The  line  moment  changes  can  be 
clearly  demonstrated  only  if  the  measurement^  are  made  down  to  a  tem¬ 
perature  where 

T<hv/k  (2-1) 

and  on  concentrated  salts  where  the  magnetic  ions  are  close  enough  to 
interact  moderately. 

We  satisfied  the  temperature -frequency  condition  (2-1)  by 

3 

working  at  30  kMc/s  and  using  a  He-cryostat  for  cooling  the  sample 
down  to  0.3°  -  0.4°K.  This  temperature  is  three  times  lower 
than  the  approximately  1°  K  one  can  reach  by  pumping  on  ordinary 

4 

liquid  He  .  Still  lower  temperatures  can  be  obtained  by 
adiabatic  demagnetization,  but  that  way  of  cooling  is  not  very  con¬ 
venient  for  resonance  experiments,  and  the  problem  of  securing  heat 
contact  between  the  sample  and  the  cold  reservoir  is  likely  to  be 

s 

serious.  Spin  resonance  work  at  a  frequency  higher  than  30  kMc/s 
nnay  require  a  stronger  magnetic  field  than  the  available  magnets  can 
give.  A  field  of  10,700  gauss  is  needed  for  resonance  in  a  salt  with 


42 

g  M  2  at  30  kMc/s  .  This  is  just  below  the  approximately  11,500 
gauss  field  limit  of  our  12"  Varian  magnet  with  2"  gap. 

For  h  V  =  30  kMc/s  ,  the  Boltzmann  factors  are 

-hv/kT 
e 

0.  71 
0.027 

At  T  =  0.  4°  K  .  all  but  a  few  percent  of  the  spins  are  in  the 
lowest  state,  while  at  4.  2^  K  .  the  levels  of  a  Kramers'  doublet  are 
nearly  equally  populated. 

For  a  Kramers'  doublet,  the  intensity  of  die  resonance  goes  as 
(I  -  e"^'’^^^)/(l  +  ,  and  the  line  area,  therefore,  changes  by 

a  factor  of  5.  5  as  we  lower  the  temperature  from  4.  2^  to  0. 4°  K  . 
If  the  line  narrows  as  expected,  the  peak  intensity  will  increase  even 
more.  This  makes  the  measurements  difficult,  since  we  may  easily 
get  too  little  resonance  absorption  for  good  signal- to  •noise  ratio  at  the 
high  temperatures  or  too  much  absorption  at  low  temperatures  for 
linear  response  of  the  system.  Spectrometer  linearity  will  be  discus¬ 
sed  later  in  Section  2.  5  . 

In  most  microwave  spectrometers,  magnet  modulation  and 
lock-in  detection  are  used  to  obtain  a  good  signal-to-noise  ratio. 
Working  with  our  concentrated  salt  we  expect  linewidths  of  at  least  a 
few  hundred  gauss,  and  then  we  could  not  use  a  modulation  scheme 
because  : 


T 

4.  2°  K 
0.4°  K 
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a)  The  field  modulation  inside  the  shielding  brass  cavity  and 
container  would  have  to  be  a  reasonable  fraction  of  the  linewidth  for 
sensitivity,  and  we  would  have  had  trouble  trying  to  modulate  this  big 
magnet  with,  say,  100  gauss  at  30  c/s. 

b)  Even  if  we  could  modulate  the  magnet,  the  eddy  currents  set 

3 

up  in  the  He  system  would  represent  a  serious  heat  input  and  prevent 
us  from  reaching  the  lowest  temperatures.  In  practice,  at  0.4°  K  we 
began  to  see  the  heating  effect  if  we  increased  the  field  modulation  to 
more  than  the  minimum  of  a  few  gauss  necessary  to  see  the  proton 
resonance  for  field  determination. 

Since  we  could  not  use  a  lock-in  system,  we  had  to  use  larger 
samples  and  record  the  chaises  in  reflected  power  as  we  swept  the 
magnetic  field  slowly  through  resonance.  This  requires  a  stable  mi¬ 
crowave  system  which  is  hard  to  obtain  at  this  high  frequency,  and  we 
found  it  necessary  to  have  the  klystron  frequency  locked  either  to  the 
sample  cavity  or  to  an  external  cavity.  Even  so,  noise  and  klystron 
amplitude  fluctuations  usually  were  of  the  order  of  one  percent  of  the 
power  level,  and  a  useful  signal  had  to  be  well  above  that. 

For  work  around  0.  5°  K  ,  problems  with  sample  saturation 
and  poor  heat  contact  require  that  the  microwave  power  be  kept  very 
low  as  will  be  discussed  in  Sections  2.  2  and  2.  3  .  At  the  lowest  tem¬ 
peratures  we  are  forced  to  use  the  more  sensitive,  but  more  compli¬ 
cated  superheterodyne  detection  scheme  instead  of  simple  crystal  de¬ 
tection. 
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The  microwave  circuits  used  were  modified  occasionally;  some¬ 
times  the  sample  was  sitting  in  the  maximum  H-field  of  a  resonant 
cavity,  and  sometimes  it  was  sitting  at  the  reflecting  bottom  of  a  wave¬ 
guide.  From  the  calculations  in  the  next  section,  it  appears  that  in  the 
latter  arrangement  the  danger  of  saturating  the  resonance  is  less,  and 
while  a  larger  sample  is  then  required,  it  is  usually  no  problem  to  get 
enough  absorption  at  this  frequency  and  temperature  range.  ' 

The  main  features  of  the  waveguide  spectrometer  are  shown  in 
Fig.  1  .  The  measuring  power  comes  from  an  EMI  9518  klystron  which 
is  frequency-locked  to  a  fixed  cavity  for  stability.  The  reflected  power 
from  the  sample  arm  is  detected  in  the  balanced  mixer,  amplified  and 
displayed  with  a  recorder  as  the  magnetic  field  is  swept  slowly  throxigh 
resonance.  The  local  oscillator  is  a  Varian  X-12  klystron  whose  fre¬ 
quency  is  doubled  in  a  crystal  harmonic  generator  made  by  Mr.  E. 
Royce.  The  X-12  frequency  is  locked  to  a  wavemeter  adjusted  so  that 
the  beat  frequency  from  the  mixer  is  30  Mc/s  .  The  frequency  stabi¬ 
lizing  circuits  were  designed  by  Mr.  S.  Prevot.  The  klystrons  are 
frequency-modulated  slightly  at  26  kc/s  or  455  kc/s  .  If  the  mean 
frequency  drifts  slightly  off  the  frequency  of  nnaximum  absorption  in 
the  cavity,  the  microwaves  will  be  amplitude -modulated  at  the  modu¬ 
lation  frequency.  A  phase -discriminating  network  then  gives  a  nega¬ 
tive  feedback  to  the  klystron  repeller,  returning  the  mean  frequency  to 
the  maximum  cavity  absorption. 
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Fig.  1  Block  Diagram  of  Waveguide  Spectrometer  . 
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A  30  db  ferrite  isolator  is  necessary  to  prevent  local  oscilla¬ 
tor  power  from  reaching  the  sample.  The  waveguide  going  into  the 
cryostat  is  the  somewhat  oversized  K  dimension;  a  tapered  transition 
connects  this  with  the  rest  of  the  system  in  Ka-band. 

The  low  temperature  system  is  shown  schematically  in  Fig.  2. 

4  o 

The  He-bath  can  be  pumped  down  to  about  1. 25  K  with  the  big  rotary 

p\imp,  ;  the  oil  diffusion  pump,  ,  could  take  it  somewhat  lower 

in  temperature,  but  has  seldom  been  used.  The  bath  temperature  can 

be  found  from  its  vapor  pressure  measured  in  the  mercury  and  oil 

4 

manometers  .  The  Cartesian  manostat  CM  can  stabilize  the  He 

pressure  at  any  level  above  a  few  nun  Hg  ,  and  this  will  keep  the  temper- 

3  X 

ature  constant  to  within  0.  Ol^K. 

3 

The  vapor  pressure  of  He  is  considerably  higher  than  that  of 

4  3 

He  at  the  same  temperature.  Also,  He  is  not  a  superfluid  with  film 

flow  so  heat  leaks  can  be  kept  low.  It  is,  therefore,  relatively  easy  to 

3 

pump  a  He  bath  down  to  about  IpmHg  pressure  which  corresponds 

to  0. 3®K  . 

3 

A  He  cryostat  similar  to  this  one  has  been  described  by  Seidel 
32  3 

and  Keesom.  The  expensive  He  gas  is  normally  kept  in  the  storage 
containers  >  ^2  ’  ^3  '  thin- 

A  Q 

walled  stainless  steel  tube  T  going  through  the  He  bath  at  1.25  K  , 
it  condenses  and  drips  down  into  a  little  container  which  is  inside  an 
evacuated  brass  can  that  also  encloses  the  sample  S  sitting  in  the 


VACUUM  UNE 


Fig.  2  Schematic  Drawing  of  He  Cryostat- 
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3 

microwave  waveguide  W  or  cavity  C  .  From  here  the  He  can  be 
pumped  back  to  storage  with  the  oil  diffusion  pump  and  the 
vacuum-sealed  rotary  pump  (Welch  1400  B),  and  in  the  process,  the 
sample  is  cooled.  The  heat  leak  together  with  the  pumping  speed  de- 

3 

termines  the  equilibrium  temperature  of  the  He -bath  which  can  be 

measured  from  the  gas  pressure  with  the  oil  manometer  and  McLeod 

gauge  M^  .  The  effective  pumpii^  speed  can  be  reduced  by  a  con  - 

trolled  leak  back  through  the  needle  valve  N  and  the  heat  input  can  be 

increased  by  dissipation  in  an  attached  electric  heater  element.  M2 

are  thermocouple  and  ionization  gauges  that  are  used  for  checking  the 
3 

vacuum  of  the  He  system  between  operations . 

Seidel  and  Keesom  measured  a  heat  removal  rate  from  their 

cryostat  by  ^He  evaporation  of  about  50  erg/s  at  T  =  0.  3**K  and 

about  1000  erg  /s  at  T  =  0.  5^  K  .  We  are  using  about  the  same 

geometry  and  pump  sizes,  and  this  cooling  capacity  should  be  sufficient, 

since  for  sample  saturation  and  heat  contact  reasons  we  will  have  to 

limit  the  microwave  dissipation  in  the  sample  to  less  than  10  erg  !  sec  . 

3  3 

In  practice,  about  1  cm  of  condensed  He  will  last  for  several  hours 

before  it  is  pumped  away  if  the  heat  leak  is  kept  small. 

Figure  3  shows  the  means  for  heat  contact  between  sample  and 
3 

He  bath.  The  sample  S  is  glued  onto  a  quartz  post  inside  a  quartz 
tube  going  through  a  hole  in  the  cavity  (or  waveguide)  bottom.  The 
quartz  tube  is  sealed  to  a  platinum  tube  which  again  is  soldered  to  a 
copper  tube  going  to  the  wall  of  the  He  container.  This  tube  system 
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SCALE  INCHES 

Fig.  3  Detail  Showing  Semple  Mounted  Ineide  Microwave  Cavit/  and  in 

3 

Heat  Contact  with  He  Reeervoir. 
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4 

can  be  filled  with  He  gas  through  a  stainless  steel  capillary  about 

4 

10  cm  long,  with  0.003"  inside  diameter.  The  He  will  condense  on 
the  tube  walls,  form  a  superfluid  film  that  flows  to  the  sample  and 
evaporates  there  to  flow  as  gas  carrying  latent  heat  back  to  the  wall 

3 

of  the  He  bath  where  it  will  recondense.  The  efficiency  of  this  heat 
contact  will  be  estimated  in  Section  2.  3  . 

It  should  be  noted  that  the  waveguide  and  cavity  are  not  cooled 
by  ^He  ,  but  remain  at  1.  25^  K  .  This  requires  care  in  centering  the 
quartz  tube  in  the  hole  so  that  it  does  not  touch  the  walls.  The  hole 
diameter  d  should  be  small  enough  to  prevent  any  microwave  power 
from  coupling  out.  The  size  of  the  quartz  tube  and  the  necessary 
clearance  forced  us  to  use  d  =  0.  43  cm.  This  hole,  when  empty,  will 
not  propagate  waves  of  frequency  30  kMc/s  ,  but  it  is  not  certain  that  the 
frequency  is  below  cut-off  when  the  tube  is  partly  filled  with  quartz  of 
dielectric  constant  3.8  .  For  frequencies  below  cut-off,  the  micro- 
wave  fields  extend  exponentially  into  the  hole,  and  the  hole  should, 
therefore,  be  reasonably  long  to  prevent  couplixig  out. 

The  temperature  was  measured  by  the  resistance  of  two 

3 

standard  10  Q  ,  1/10  W  resistors  R  calibrated  against  He  vapor 
pressure.  One  was  glued  onto  the  ^e  can  and  one  was  attached  to 
the  quartz  tube  as  close  as  possible  to  the  cavity  hole.  Whenthe  micro- 
wave  power  was  kept  sufficiently  small, we  saw  no  temperature  dif  ^ 
ference  between  the  two  thermometers.  From  the  discussion  in 
Section  2.  3  ,  it  will  appear  that  the  most  likely  bottleneck  for  heat 
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transfer  from  the  sample  is  the  He  gas  flow  in  the  tube.  If  the  two 
resistors  read  the  same  during  measurements,  it  is  then  reasonable 
to  believe  that  the  sample  is  not  significantly  heated. 

In  Fig.  3  the  seal  E  is  for  electrical  leads  to  thermometers 
and  heater.  The  waveguide  also  serves  as  a  pumping  tube  for  evacua¬ 
tion  of  the  brass  can,  and  it  is  vacuum-sealed  on  top  with  a  seal  of 
mica  onO-ring.  A  piece  of  polyfoam  in  the  waveguide  prevents  room 
temperature  radiation  from  coming  down  and  heating  the  quartz  and 
sample . 

The  samples  were  grown  from  water  solutions  at  room  tem¬ 
perature.  Some  of  the  crystals  could  fairly  easily  be  ground  into 
spheres  by  blowing  them  around  in.a  cylinder  lined  with  emery  paper. 
However,  crystals  like  NiSiF^  •  6H2O  that  would  easily  split  along  the 
axial  direction  usually  turned  into  ellipsoids,  and  in  others  like  CuK2Cl^  * 
2H2O  ,  the  crystal  lattice  seemed  to  distort  under  this  treatment.  Crys  ¬ 
tals  of  the  latter  type  could  be  made  into  approximate  spheres  by 
carefully  dissolving  away  corners  with  wet  filter  paper.  The  sample 
size  needed  in  a  waveguide  was  about  1  mm^  ,  in  a  cavity  100  times 
less,  and  such  small  sizes  of  soft  crystal  were  rather  hard  to  handle. 
The  samples  were  x-ray  oriented  and  glued  on  the  quartz  rod  with 
G£  1202  varnish  or  Shell  Epon  828  epoxy.  The  crystal  orientation  in 
the  horizontal  plane  could  be  checked  by  rotating  the  magnet  and  ad¬ 
justing  its  angle  for  observation  of  extremal  g- value  or  narrowest  or 
broadest  resonance  line  as  required  by  the  particular  salt. 
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It  it  believed  that  crystal  orientation  usually  was  better 

•I  -O 

taan  +  2 

2.  2.  Discussion  of  possible  microwave  saturation. 

The  spectrometer  must  operate  at  a  low  enough  microwave 
level  so  that  the  spin  system  is  kept  essentially  at  the  temperature  of 
the  ^He  bath.  This  will  be  most  critical  at  the  lowest  temperatures 
where  the  spin -lattice  relaxation  time  is  long  and  the  heat  contact  be¬ 
tween  crystal  lattice  and  bath  is  poor.  For  the  design  of  the  spectrom 
eter,  it  is  helpful  to  estimate  the  order  of  magnitude  of  the  heat  trans¬ 
fer  bottlenecks  to  see  what  the  permissible  microwave  power  levels 

are  and  whether  a  superheterodyne  detection  method  is  necessary. 

12 

To  avoid  spin  saturation,  we  require 

Y^TjT2Hj^<<1  ,  (2-2) 

where  Y  =  gp/1l  =  gyromagnetic  ratio  , 

Tj  =  spin-lattice  relaxation  time, 

T^  =  spin-spin  relaxation  time, 

=  amplitude  of  circularly  polarized  microwave  magnetic 
field  at  the  sample  . 

T^  is  related  to  the  width  of  the  line  by 

GO 

0 


(2-3) 


For  a  Lorentzian  line  shape,  ZT^  equals  the  width  between  half-power 
points. 

The  microwave  field  amplitude  for  a  given  input  power  will 

now  be  calculated  and  compared  in  the  cases  of  cavity  spectrometer 

and  nonre sonant  waveguide  spectrometer. 

a)  We  used  a  cylindrical  cavity  having  the  dimensions:  radius 

a  =  8.  25  mm  ,  length  Zq  =  14.  8  mm  .  It  was  operated  in  the 

mode,  and  the  sample  was  sitting  on  the  vertical  aucis,  one-qtiarter 

wavelength  from  the  bottom.  The  field  distribution  in  the  cavity  in 

33 

terms  of  the  microwave  field  at  the  sample  is 

Hz  =  Jq  (u'f  )  cos  p  z  (2-4) 

I 

”r  “  ”s  -^l  I  ^  P  *  * 

I 

where  u  s  3.  832  =  first  zero-point  argument  of  Bessel  function 

Jj  (r) 

,  i 

c  3  speed  of  light 
(d  s  2ir  *  30  kMc/s 

The  energy  content  of  the  cavity  is 
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1 

i 

I 

\ 

\ 

I 

=  Zq  (^-r)^]  (  Jq('1  ))^ 

=  18.6  •  10"^  H  ^  erg  (2-6) 

8 

The  Q  of  the  cavity  was  never  accurately  measured,  but  in  order  of 
magnitude  it  was  about  10,000  at  helium  temperatures. 

The  wall  losses  in  terms  of  the  microwave  field  at  the  sample 
would  then  be 

Wg  =  ^  =  3.5-  10^  erg/s  .  (2-7) 

For  a  critically  coupled  cavity,  the  wall  losses  would  equal  the 
input  power  W.  .  In  this  case  , 

W.  i 

H  =  ( - i - 1^  )  ^  gauss  (2-8) 

®  3:  5  •  10^ 

To  avoid  saturation,  we  demand  that 

-  ,2-2  H  , 

Y^TjT^Hj  =  S-|-  TjT2(-2^)  f 

Tl 

i 

s 

where  the  factor  2  under  H  comes  in  because  the  is  linearly  j 

9  8 

polarized.  Substituting  for  H  from  (2-8),  we  then  find  the  maximum 

8 

input  power  as 

w.  <  erg/.  .  (2-10) 

‘  “  g2  Tj  t^ 
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-1  -8 

For  the  reasonable  values  Ti  *  10  sec  ,  T ^  *^10  sec  , 


1 


g  »  2  ,  this  becomes 


W.  <  0.45erg/s  i^O.OSjiW 


(2-11) 


b)  The  power  flow  in  the  dominant  mode  of  a  waveguide  in 


terms  of  the  maximum  transverse  magnetic  field  is 

c  H  .  ^  a^  brt 

SI  0  0 


.  34 


(2-12) 


For  a  standard  K-band  waveguide,  the  dimensions  are  : 
a^  =  1.  07  cm  ,  bQ  =  0.432  cm  .  If  the  sample  is  sitting  in  the  maxi¬ 
mum  field  on  a  perfectly  reflecting  waveguide  short,  we  have 


=  2H. 


SI 


In  this  case  the  saturation  condition  (2-9)  becomes 


W,  < 


i  T?r' 


‘o^o' 


16»  T, 


s  100  erg/s  =:  10  |iW  (2-13) 

-1  -8 

when  we  again  use  ^  10  sec  ,  ’1*2^  >8^  ^  * 

The  sensitivity  of  a  spectrometer  is  proportional  to  the  sqxiare 
of  the  field  amplitude  at  the  sample.  Here  we  have  estimated  the 
cavity  spectrometer  to  be  about  200  times  more  sensitive  than  the 
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waveguide  spectrometer.  This  is  in  fair  agreement  with  our  experience 
that  the  sample  for  the  waveguide  would  have  to  be  about  100  times 
bigger  than  the  sample  for  the  cavity,  and  it  does  not  necessarily 
contradict  our  assumption  of  a  cavity  Q  10,000  ,  since  the  field 
configurations  in  the  two  spectrometers  are  completely  different. 

It  is  doubtful  from  this  order  of  magnitude  calculation  whether 
a  cavity  spectrometer  can  be  safely  used  with  only  simple  crystal  de¬ 
tection  The  signal-to-noise  ratio  from  the  detector  will  probably  be 
too  poor  if  the  microwave  power  level  has  to  be  kept  below  0.1  . 

But  straight  detection  should  be  possible  without  danger  of  saturation 
if  a  waveguide  spectrometer  is  used  with  its  smaller  microwave  field. 
This,  however,  means  that  we  have  to  use  a  bigger  sample,  more 
power  is  dissipated  in  it,  and  lattice-bath  thermal  bottlenecks  may 
become  important.  This  will  be  discussed  in  the  next  section. 

3 

2.3.  Discussion  of  heat  contact  between  sample  and  He  system. 

We  would  like  to  choose  the  sample  size  so  that  about  10  ^  of 
the  incident  microwave  power  is  absorbed  at  resonance.  According  to 
our  guess  in  the  last  section,  this  absorbed  power  may  perhaps  amount 
to  1  before  we  start  saturating  the  transition  in  the  waveguide 
spectrometer  case.  Now  we  want  to  calculate  whether  we  can  safely 
dissipate  this  power  without  heating  the  sample.  We  should  then  be 
forewarned  that  even  rough  order  of  magnitude  calculations  may  be 
misleading  since  heat  transfer  below  l^K  is  complicated  by  little- 
known  material  constants  and  surface  effects.  Also,  the  tube  connecting 
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the  sample  and  the  He  reservoir  has  a  complicated  geometry,  and 
to  simplify  the  calculation  we  have  to  assume  that  this  is  equivalent  to 
the  same  length  of  uniform  tubing. 

We  may  visualize  the  heat  transfer  as  being  divided  in  these 
separate  processes  : 

a)  Heat  conduction  from  the  interior  of  the  sample  to  the 
surface. 

b)  Heat  resistance  on  the  surface  of  the  sample  due  to  phonon 
mismatch. 

4 

c)  Evaporation  of  the  He  film  on  sample,  which  has  to  be 
replenished  by  film  flow. 

d)  Flow  of  the  gas  containing  the  heat  of  vaporization  through  a 

3 

narrow  tube  that  offers  resistance  back  to  the  wall  of  the  He  reser> 
voir  for  recondensation.  Where  the  tube  has  thick  copper  walls,  heat 
conduction  in  the  copper  may  reduce  the  necessary  gas  flow  consider¬ 
ably. 

It  will  turn  out  from  the  following  estimates  that  effect  d)  is 
the  bottleneck  that  prevents  us  from  reaching  really  low  temperatures 
during  measurements. 

f 

a)  Since  we  expect  the  resonating  spins  to  give  up  energy  in 
e 

quanta  of  hv  ,  we  might  look  for  the  energy  carrying  properties  of  a 

e 

band  of  "hot  phonons"  around  V  which  would  take  the  dissipated 
energy  to  the  sample  surface.  However,  there  are  conflicting  opinions 
about  the  reality  of  such  "hot  phonons"  in  most  cases  of  spin 
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relaxation.  It  appears  likely  that  the  Zeeman  energy  is  mostly 
transferred  to  phonons  of  higher  frequency  via  impurities,  spins  at  ir¬ 
regular  lattice  points  and  strongly  exchange-coupled  pairs.  Presum¬ 
ably,  the  phonon  spectrum  comes  rapidly  into  thermal  equilibrium  in 

4 

most  relaxation  experiments  at  liquid  He  temperatures,  and  the 

energy  transfer  through  the  lattice  can  then  probably  be  treated  as 

simple  heat  conduction.  We  will  here  assume  that  this  is  valid  even 

for  temperatures  below  0.  5°K  where  the  phonon  mean  free  path  is  of 

the  same  order  as  the  sample  dimensions. 

36 

Casimir  ”  has  derived  a  formula  for  the  heat  conductivity  K 
in  pure  crystals  when  only  surface  scattering  of  phonons  is  important. 
Measured  values  of  conductivities  are  few  in  this  temperature  range, 
but  they  are  generally  somewhat  smaller  than  predicted  by  Casimir' s 
theory.  At  T  =  0.  3°K  Garrett^^  found  K  *  lO^ergsec  ^cm  ^deg^for 
K  Cr  {SO^)^  '  *  which  corresponds  to  a  temperature -indepen¬ 

dent  phonon  mean  free  path  of  0.5  mm  .  We  will  assume  that  the  heat 

conductivity  is  approximately  the  same  for  our  samples  of  hydrated 

4  3  -1 

crystals.  A  uniform  heat  dissipation  of  P  w  10/(y  ir  a  )  erg  sec 

_3 

cm  in  a  sphere  of  radius  a  will,  in  the  stable  state,  give  a  tem- 

38 

perature  difference  between  the  center  and  the  surface  of 

AT  =  Pa^/6K  SO.  0005°  for  a  =  0.  075  cm"^  .  (2-14) 

This  is  insignificant, and  unless  this  estimate  for  the  heat  con¬ 
ductivity  is  too  small  by  more  than  an  order  of  magnitude,  heat 
trapping  inside  the  sample  does  not  seem  to  be  important. 
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b)  Phonon  mismatch  at  the  boundary  between  crystal  and  He 

film  will  appear  as  a  thermal  resistance.  The  mismatch  will  de¬ 

pend  on  the  vibration  spectra  of  the  two  media,  but  the  resistance  even 

4 

between  polished  copper  and  bulk  He  does  not  agree  too  well  with  the 
existing  theories.  In  our  case  we  have  additional  complications  :  First, 
the  sample  surface  is  rather  rough  from  grinding,  and  part  of  it  is 
covered  with  glue  of  unknown  properties.  Second,  the  film  thickness  is 
of  the  order  of  the  wavelength  of  the  excited  phonons  at  this  tempera¬ 
ture.  The  phonon  spectrum  is,  therefore,  likely  to  be  changed  relative 
to  that  of  bulk  liqmd. 

Instead  of  relying  on  the  theory,  we  may  then  as  well  use  the 

4 

heat  transfer  rate  between  copper  and  bulk  liquid  He  as  measured  by 
Fairbanks  and  Wilks^^  in  the  range  0.  3  to  l^K 

cm^  deg/ watt  (2-15) 

For  W  =  10**^  watt,  T  m  0.4°K  and  a  sample  surface  of 
«  4ir  a^  where  a  =  0.  075  cm  ,  we  get 

dT  0.004°K 

The  boundary  resistance  seems  to  be  unimportant. 

c)  The  critical  flow  in  a  saturated  ^He  film  is  proportional  to 

the  perimeter  of  the  narrowest  constriction  of  the  flow.  At  T  w  0.  5^ 

-4  3  42 

the  liquid  flow  rate  is  about  10  cm  /(sec  cm).  The  specific  weight 
4  3 

of  He  is  0. 1455  g /cm  and  its  molecular  weight  4*,  so  with  a 
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smallest  diameter  taken  to  be  m  0.  05  cm  ,  the  maximum  flow  rate 
should  be 

F  «  3.14  •  10~®  cm^/sec  =  1.15  •  10'^  mole/sec. 

The  heat  of  vaporization  of  He  is  about  70  joules/mole  at 
o  43 

T  a  0.  5  K  .  Therefore,  to  carry  away  an  energy  of  1  fiW  a  flow 
of  0. 14  ■  10~^  mole/sec  is  needed,  and  the  film  flow  should  easily  take 
care  of  this. 

4  44 

d)  The  saturation  pressure  of  He  is  given  by 

'“Sio  ■■mm  Hg  =  + 1  ‘“*10  •  ( *-  “•> 

Ok  k 

Other  terms  appearing  in  more  accurate  expressions  for  p(T)  are 
negligibly  small  below  1°K  .  In  this  empirically-fitted,  theoretical 
equation,  p  and  T  have  to  be  taken  as  dimensionless  quantities  in  the 
units  indicated. 

The  mode  of  flow  of  the  evaporated  gas  back  to  the  cold  reser- 

4 

voir  will  depend  upon  whether  the  mean  free  path  of  the  He  atoms  in 
the  vapor  is  smaller  or  larger  than  the  diameter  of  the  tube.  However, 
it  is  here  hardly  appropriate  to  speak  of  a  mean  free  path  as  an  ana¬ 
log  to  the  classical  case  of  hard,  elastic  spheres  colliding  with  each 
other.  For  helium  molecule  collisions  at  low  temperature  wave  me¬ 
chanical  resonance  effects  appear  and  the  cross'^sections  vary  widely 
with  the  energies  of  the  colliding  particles.  To  see  the  order  of  mag¬ 
nitude  of  the  mean  free  path,  we  may  then  as  well  use  the  following 
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equation  given  by  Keesom,  although  the  empirical  formula  for  the 
viscosity  used  in  deriving  it  is  not  very  good  below  3°K  . 

=  2.159  •  10  ^  i -  cm  .  (2-17) 

”mm  Hg 

For  T  =0.5°  we  find  from  (2-16)  ,  p  =  16. 10"^  mmHg,  and 
the  mean  free  path  is  approximately 

0  5°  K  *  ^  ■ 

The  vapor  pressure  is  varying  very  rapidly,  and  at  T  =  0.  4°  K 
and  T  =  0.6°K  we  find 


4  ^  =  30  cm 

0.4°K 


*0.6°K 


0.  04  cm 


We  see  that  below  T  m  0.  5°K  the  gas  flow  has  to  be  treated  as 
free  molecular  diffusion  or  Knudsen  flow.  The  formula  for  the  flow 
through  a  tube  of  radius  a  and  length  L  under  a  pressure  of  tP 
-  (Pi  “  P2) 

^freemol^-^  ^  ^  mole/sec  (2-18) 

where 

7 

R  =  8,317' 10  ergs /(deg  mole)  =  gas  constant 

4 

M  s  4  g/mole  for  He 


Equation  (2-18)  is  derived  under  the  condition  that  the  mole¬ 
cules  are  re -emitted  in  completely  random  directions  after  qpch  col¬ 
lision  with  the  film-covered  tube  wall.  If  the  molecules  were  reflected 
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from  the  film,  keeping  part  of  their  forward  momentum,  the  gas  flow 
could  be  considerably  larger.  In  view  of  the  molecule -molecule  col- 

4 

lilion  anomalies  in  the  He  at  this  very  low  temperature,  it  might 
perhaps  be  interesting  to  try  to  check  experimentally  by  heat  conduc¬ 
tion  in  narrow  tubes  whether  molecule -film  collisions  are  specular 
reflections  or  essentially  condensations  and  re -evaporations. 

From  (2-16) 


P 


19260^  ^2 

3.1077T 


10 


2.06 


dyn/cm^ 


(2-19) 


By  differentiating 

I 

c  (|  T  +  3.107  fn  10)  AT  , 

Ap  =  2.06  •  10^  - "  q^.107/T -  dyn/cm^  .  (2-20) 


Combining  (2-20)  with  (2-18)  and  multiplying  with  the  heat  of  vaporiza 
tion  C  =70  joule/mole,  we  get  the  flow  of  latent  heat  : 


^free  mol  "  ^  ^ 


0.  265  •  10^  a^  (|t  +  7.17)  AT 


L  10 


Twrrr 


watt  . 


(2-21) 


The  total  length  of  the  tube  connecting  the  sample  to  the  cold 

3 

wall  of  the  He  can  is  about  10  cm  .  Most  of  it  is  made  of  copper. 
The  heat  conductivity  of  the  copper  walls  will  vary  greatly  with  metal 
purity,  heat  treatment,  and  machining,  and  soldered  joints  will  reduce 
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it.  However,  by  extrapolating  from  the  data  summarized  by  White, 

we  may  take  the  copper  heat  conductivity  to  be  of  the  order  0. 1  W 

cm'*'  deg'*  at  0.  5°K  .  A  tube  8  cm  long  and  with  w  10  mm*  aver- 

-5 

age  wall  cross  section  will  then  conduct  about  10  watt  for  a  tempera¬ 
ture  difference  of  0.01°  between  ends.  We  can,  therefore,  assume  that 
the  gas  has  to  flow  only  through  the  quartz  tube  to  the  beginning  of  the 
copper  tube.  In  (2-21)  we,  therefore,  take  L  2.  5  cm  ,  a  »  0.1  cm  , 
and  if  we  specify  that  the  temperature  difference  6T  is  to  be  kept  be¬ 
low  0.  01°  ,  the  permissible  heat  input  is 

T  =  0.30  0.40  0.50  (0.60)°K 

Qfree  ^ol  -  ‘  ‘  ’ 


At  the  temperatures  above  0.  5  K  it  is  more  correct  to  treat 
the  flow  of  gas  as  viscous  or  Poiseuille  flow.  The  formula  for  this 

u  « 


ira  pAp  , 

'  StiRTL  ^  free  mol 


(2-22) 


The  last  term  in  this  formula  is  a  theoretical  correction  term 
for  "slip"  along  the  tube  walls.  t]  is  the  viscosity  coefficient  which 
we  take  for  ^He  gas  around  0.  6°K  as  n  3. 10'^  poise.  Plugging 
(2-19)  and  (2-20)  into  (2-22)  and  using  the  previous  values  of  a  ,  L  , 
and  C  ,  we  get  for  the  flow  of  latent  heat  in  this  case  : 
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^visc  "  ^visc  ■ 

rl8.6  •  10^  (|  T  +  7.  17)AT  °  63  {|  T  +  7.  17)  AT 

=  L - JJ57ZRTT -  ^  — [^nUTTT - J  watt.  (2-23) 

If  we  again  specify  AT  ^  0.  01°  ,  we  get  the  values 

T  =  0.6  0.7  0.8  0.9  1.0°K 

Q  .  <  6.3  •  10“^  1.3  •  10'^  1.8  •  10“*  1.7  •  10"^  1.  MO'^watt 

vise  — 

It  should  again  be  stressed  that  the  calculated  numbers  repre¬ 
sent  at  best  only  the  order  of  magnitude  heat  flow  in  the  actual  case, 
since  we  can  only  guess  the  magnitude  of  the  material  constants  in¬ 
volved.  and  perhaps  the  formulas  used  are  not  even  approximately 
valid.  But  in  spite  of  these  uncertainties,  it  seems  clear  that  the  gas 
flow  presents  the  most  serious  bottleneck  to  the  heat  transfer  from  the 
sample,  and  that  a  |iissipation  of  10'^  W  may  appreciably  heat  up  the 
sample  below  about  0.  6^  Only  with  a  heat  input  of  about  10  ^  watt 
or  less  can  we  hope  to  take  measurements  at  0.4°K  ,  and  this  certainly 
requires  a  superheterodyne  system. 

2.4.  Heat  leak  through  capillary. 

4 

The  He  film  used  to  cool  the  sample  is  supplied  through  a 
thin-walled,  stainless  steel  capillary.  Flow  of  gas,  film  and  perhaps 

4 

bulk  liquid  He  through  this  capillary  will  represent  a  heat  leak  from 
the  ^He  bath  at  approximately  1. 25°K  to  the  ^He  system.  We  can 
reduce  this  heat  leak  by  pumping  the  gas  away  from  the  warm  end  of 
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the  capillary.  When  the  warm  end  is  not  pumped,  we  estimate  the  heat 

-4  3 

leak  to  be  about  10  watt  from  the  observed  rate  of  He  evaporation, 

and  we  get  only  down  to  about  0.  in  temperature. 

The  capillary  is  about  10  cm  long  and  has  an  internal  diameter 

of  0.003  "  (a  =  3.  8  '  10  cm).  If  we  apply  equations  (2-18)  and  (2-22) 

to  the  problem  of  heat  conduction  in  this  capillary,  we  find  that  the 

temperature  gradient  is  concentrated  very  close  to  the  cold  end  which 

3 

is  kept  at  the  He  temperature.  The  heat  leak  from  film  and  gas 
flow  calculated  this  way  is  of  order  10~^  watt  when  the  warm  end  is 
not  pumped.  This  may  seem  surprisingly  large  compared  to  the 
rather  poor  heat  contact  we  calculated  for  the  sample.  It  is  a  result  of 
the  good  heat  conduction  in  the  larger  part  of  the  capillary  which  has 

4 

almost  the  temperature  of  the  He  bath.  However,  it  does  not  ex¬ 
plain  the  observed  heat  leak  D^ich  is  larger  by  two  orders  of  magni¬ 
tude. 

We  must  then  assume  that  the  capillary  and  the  sample  tube 

4 

will  get  filled  with  liquid  He  if  we  do  not  pump  the  gas  away  after 

condensation.  The  flow  of  heat  in  -filled  capillaries  can  be  con- 

50 

siderable,  but  it  is  a  nonlinear  fxinction  of  the  temperature  gradient, 
and  no  estimate  of  such  a  heat  leak  will  be  attempted.  We  should  note, 
however,  that  if  the  sample  tube  is  filled  with  liquid  helium,  there 
should  be  no  problem  with  heat  conduction  away  from  the  sample. 

If  we  pump  on  the  warm  end  of  the  capillary,  we  may  be  able  to 
bring  the  gas  pressure  down  to  the  saturation  pressure  of  ^He  at 
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T  Si  0.  8  K  .  In  this  case,  the  heat  leak  will  be  given  by  the  gas  flow 
formula  (2-18)  which  we  now  write 

0.265  •  10“*  a^  (|t  +  7.  17)  dT 


Q  = 


10 


3.io7rr 


dL 


12.4  •  10 


dT 


dL 


(2-24) 


where  we  have  used  (y  T  +  7 . 17)  »  8.  5  . 

We  take  the  cold  end  of  the  capillary  to  be  at 

integrate  along  the  tube 

L  0.  8 

J  Q  dL  =  12.4  •  10"“* 


«0.4°K  and 


0.4 


10 


dT 

rwr 


(2-25) 


The  integral  on  the  right  side  cannot  easily  be  solved.  But 

since  the  overwhelming  temperature  dependence  is  in  the  exponent, 

2  2 

we  multiply  under  the  integral  by  0.  6  /T  w  1  .  Then  we  get  ap¬ 
proximately 

JQ-3.107/T 
3.107  -tn  10 

=  0.8*10"®  .  (2-26) 

For  L  s  10  cm  we,  therefore,  calculate  the  approximate  heat  leak 
when  the  gas  is  pumped  away 

•  _Q 

Q  10  watt. 


0.  8 

QL  «  12.4  •  10"“*  -0.36  J 

0.4 
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Now,  helium  is  removed  from  the  sample  tube  at  the  critical 

-4 

rate  of  film  flow  in  the  capillary.  The  critical  flow  rate  is  a  ^10 
3  42 

cm  /(cm  sec),  and  for  the  capillary,  this  corresponds  to 

-3  3 

F  =  a.Zira  =  8.6*  10  cm  /hour 

Filling  the  sample  tube  goes  somewhat  faster,  since  it  is  done 
by  gravity  flow  of  bulk  helium.  Once  the  tube  is  filled,  we  should  be 
able  to  work  for  a  few  hours  before  re-filling,  and  this  is  also  our 
practical  experience. 

2.  5.  Discussion  of  spectrometer  linearity  and  large-signal  corrections. 

As  shown  in  Section  2.1  ,  we  require  that  the  spectrometer  have 
a  fairly  large  dynanciic  range  udiere  the  output  preferably  is  a  linear 
function  of  the  imaginary  part  of  the  susceptibility  X  ",  or  is  at  least 
a  known  function  of  X  "  .  The  large-signal  behavior  of  a  cavity  and 
waveguide  spectrometer  will  now  be  discussed. 

First,  we  note  that  it  does  not  matter  for  the  linearity  (but 
certainly  for  the  sensitivity)  whether  we  use  direct  detection  or  super¬ 
heterodyne  amplification  followed  by  a  low-level  square-law  detector. 

In  both  cases  the  final  recording  will  be  proportional  to  tiie  nriicrowave 
power  of  the  detector.  The  linearity  of  the  detecting  system  can  be 
simply  checked  by  inserting  and  varying  a  precision  attenuator,  which 
has  been  done. 

Secondly,  any  attempt  to  correct  for  large-signal  nonlinearity 
resulting  Drom  too  large  a  sample  will  be  very  unreliable  because  we 
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then  violate  the  requirement  that  the  presence  of  the  sample  be  only  a 
small  perturbation  on  the  microwave  field.  A  big  field  perturbation 
will  have  the  following  effects  : 

a)  The  absorption  in  the  sample  is  proportional  to  X. 

where  is  the  microwave  field  at  the  sample.  For  a  large  perturba¬ 
tion,  the  average  may  change. 

b)  With  a  different  field  distribution  the  wall  losses  may  change. 

c)  If  the  sample  in  a  cavity  perturbs  the  field  at  the  coupling 
iris,  the  cavity  coupling  may  change. 

It  is  very  hard  to  estimate  the  importance  of  the  effects  a)  -  c) 
in  a  given  experiment.  If  for  part  of  the  measurements  large  signals 
cannot  be  avoided,  the  following  formulas  for  the  response  of  the 
spectrometer  should  only  be  taken  as  an  indication  of  how  the  output  is 
likely  to  change. 

A.  Cavity  spectrometer. 

The  magnetic  resonance  is  accompanied  by  dispersion  which 

will  change  the  frequency  of  the  loaded  cavity.  But  when  the  klystron 

frequency  is  coupled  to  the  cavity,  this  will  have  no  first-order  effect. 
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The  power  absorbed  in  the  sample  is 

wX”  J  dV  (2-27) 

where  w/ 2v  is  the  frequency  and  V  is  the  volume  of  the  sample. 
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-  '> 

+  *  Tj2 - 5  •*.  „  Q  X. 


4(2  - 

4  -a - 2-^  X")^  +  ...  > 

(^4.  J 


(2-31) 


where  is  the  incident  power. 

We  see  that  for  linear  response  we  want 


(2  -  ^  ) 


(2-32) 


and  in  particular,  we  must  avoid  critical  coupling  B  =  Q  /Q  -  I 

^  *^0  ex  o 

for  the  unloaded  cavity.  The  second -order  term  is  seen  to  vanish  for 

~  ^  which  corresponds  to  an  unloaded  cavity  reflection 

of  P  =  P  (1  -  B  )^/(l  +  B  K  P  we  should  try  to  use  an 

r  o  *^0  *^0  7  0 

iris  corresponding  to  this  coupling. 

The  ratio  can  be  easily  measured  on  an  oscilloscope 

when  the  klystron  frequency  is  swept  through  its  mode.  Doing  this  for 
no  magnetic  resonance  and  maximum  magnetic  resonance,  say,  we  can 

II 

then  calibrate  the  chart  records  and  find  the  relative  values  of  X 
from  (2-31).  The  reflected  power  according  to  formula  (2-31)  is 
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plotted  in  Fig.  4  for  varying  coupling  and  magnetic  absorp  - 
tion. 

B.  Waveguide  spectrometer. 

For  a  crystal  mounted  close  to  the  bottom  of  a  perfectly-re¬ 
flecting  waveguide  short,  dispersion  is  not  expected  to  have  any  first- 
order  effect  on  the  reflected  power. 

A  wave  traveling  through  an  absorbing  medium  will  be  attenu¬ 
ated  exponentially,  and  we  have 

II 

P  =  P  e  (2-33) 

Since  the  crystal  is  small  compared  to  the  wavelength  and  is  sitting  in 
a  node  of  the  magnetic  field,  it  may  be  more  appropriate  to  use  the 
following  formulas  : 

The  field  amplitude  at  the  crystal  is 

^  +  /PJ  .  (2-34) 

The  absorbed  power  P^^  is 

**a  =  ^ 


We  also  have 


P  *  P  -  P.  . 
r  o  a 

Solving  these  equations  we  get 

2  -  b^  X  -  b  X  y  4  -  3  b^X 


(2-36) 


)  =p  [ 
r  o  *• 


2  (1  +b  X”  )  ^ 
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Fig.  4  Reflected  Power  from  Cavity  as  Function  of  Absorption  in  Sample  . 
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II  Q  ?  Ats  'K 

a  P^[l-3bX  +  j  h  t  +^b^)(  +...  ].(2-37) 

Including  terms  to  second  order,  this  expansion  is  identical  to 
the  expansion  of  (2-33)  with  a  =  3b  ,  and  it  does  not  matter  much  then 
whether  we  use  correction  formula  (2-33)  or  (2-37)  . 

C.  Waveguide  spectrometer  with  microwave  bucking. 

In  order  to  reduce  the  dc  bias  on  the  signal  in  a  waveguide 
spectrometer  and  reduce  the  klystron  noise,  one  sometimes  balances 
or  "bucks  out"  part  of  the  microwave  power  in  a  bridge.  This  may 
take  the  form  of  a  T  where  the  sample  sits  in  one  of  the  equivalent 
arms  and  the  other  arm  contains  an  adjustable  attenuator  and  a  sliding 
short.  Then  we  can  add  a  signal  of  arbitrary  phase  and  amplitude. 

Such  a  system  will  be  very  similar  to  an  overcoupled  cavity. 

As  for  a  cavity,  the  magnetic  dispersion  in  the  sample  will  change  the 
tuning  and  the  klystron  frequency  should  be  locked  for  right  balance. 

From  the  same  reasoning  as  in  Part  B  ,  the  output  of  a  balanced 
system  is 

P  =  P  {J  +  e’^^TF"  )^  (2-38) 

r  o 

The  balancing  signal  of  relative  power  b  should  be  added  in 

•  1 

opposite  phase,  so  that  e~  *'  =  -1  .  Expanding  the  exponential  in  this 
case,  we  get 

2 

P^  =  P^  [  (1  +  b  -  2/B" )  -  (I  -  )  a^  X  +  . .  .  ] 

(2-39) 
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II  1 

For  complete  balance  when  X  =  0  ,  we  need  to  have  b  s 

II 

and  as  for  a  cavity  we  have  then  no  linear  term  in  X  in  the  response. 
The  second-order  term  will  only  vanish  when  b  =  4  ,  and  this  requires 
some  attenuation  in  the  sample  arm. 

In  practice,  it  is  hard  to  determine  the  amount  of  bucking; 
proper  large -signal  correction  is,  therefore,  very  difficult  and  should 
not  be  attempted. 
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CHAPTER  3 


measurements  on  neodymium  ethyl  sulfate 

3.1.  Description  of  neodymium  ethyl  sulfate. 

The  rare  earths  form  a  series  of  isomorphic. hexagonal  ethyl 

sulfates,  M(C2HgSO^)^'  9H2O.  and  these  salts  have  been  widely 
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used  in  paramagnetic  experiments.  Ketelaar  determined  the 
structure  from  x-ray  studies,  and  gave  the  space  group  as  The 

unit  cell  is  shown  in  Fig3a,  and  the  two  rare  earth  ions  are  at  the  posi¬ 
tions  (y  y  ^  dimensions  of  the  unit  cell  vary  slightly  with 

the  size  of  the  rare  earth  ion;  for  Nd(C2HgSO^)^  *  9H2O  it  is  a^ 

O  O 

=  13. 992  A  and  c  =  7.07  A.  Each  Nd  ion  will,  therefore, 
o 

o 

have  two  nearest  Nd  neighbors  7.07  A  away  along  the  hexagonal  axis 

o 

as  indicated  in  the  drawing.  The  other  ion  in  the  same  cell  is  8.85  A 
away  in  a  direction  making  an  angle  0  =  66.4°  with  the  hexagonal  axis. 

Each  Nd^^^  ion  is  surrounded  by  9  water  molecules,  and  the 
arrangement  is  probably  as  shown  in  Fig.  5b  .  The  water  dipoles,  to¬ 
gether  with  some  contribution  from  the  more  distant  surroundings, give 
an  electric  field  of  axial  symmetry  at  the  position  of  the  rare  earth 
ion.^  The  free  ion  Nd^  has  the  ground  state  4f  ^9/2*  flte  field  will 
split  this  into  five  Kramers'  doublets,  and  the  field  symmetry  is  reflected 

w 

in  the  axial  symmetry  of  the  g-tensors.  The  first  excited  doublet  is  about 
130  cm~^  above  the  lowest  doublet,  and  only  the  latter  will  be  populated  at 


I 


5b  Probable  Arrangement  of  Water  Molecules  Around  die  Neodymium 
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liquid  helium  temperatures .  The  magnetic  behavior  is  the  same  for  both 
ions  in  the  unit  cell,  and  it  is  given  by  the  approximate  spin-Hamilton- 
ian  for  spin  S  =  . 

K  =  6g„HS,  +Bg  (HS  +  HS)+AIS 
z  z  '  X  X  y  y  a  a 

+  B  (S^  +  Sy  ly)  +  P  -  31  (I  +  1)  )  .  (3-1) 

The  g -values  for  neodymium  ethyl  sulfate  measured  in  salts 
where  the  have  been  diluted  1:200  with  unmagnetic  La'^^^^are^^ 


%ll  =  3.535  +  0.001 


2.072  +  0.001 


The  even  isotopes  of  Nd  of  mass  142,  144,  146,  148,  and  150, 
which  have  a  total  natural  abundance  of  79.  5  ^  ,  have  no  nuclear  spin 
and  ,  for  these  ,  the  three  last  terms  in  (3-1)  vanish.  The  isotopes 
143  (abundance  12.  2  ^)  and  145  (abundance  8.  3  ^)  both  have  spin 
1  =  7/2  which  will  give  a  number  of  weak  hyperfine  lines.  Their  po¬ 
sitions,  when  the  field  is  parallel  to  the  crystal  axis,  can  be  given  more 
54 

accurately  by 


2  2 

g//PH  =  hv  .Am(l+^.j]  [I(I  +  l)-m^] 

4(hv) 


[  2m^  -  m  (21^  +  21-1)] 


2(hv) 


(3-2) 


where  m  is  the  nuclear  magnetic  quantum  number,  and  the  interaction 
constants  are  : 
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Isotope  10^A(cm"^)  10^B(cm  10^P(cm 

143  380. 3  +  0  1  198. 9  +  0.5  <2 

145  236.4  +  0.1  123.7  +  0.5  <2 

In  concentrated  salts  we  expect  the  spin-Hamiltonian  constants 

to  have  about  the  same  values  as  in  the  diluted  salts.  We  find  a 

little  larger,  however;  g^^  =  3.60  ,  and  this  value  will  be  used  in  the 

following  calculations  But  a  new  effect  also  appears  here.  The 

main  resonance  line  splits  up  into  three  almost  separated  peaks  when 

measured  with  the  external  field  along  the  hexagonal  axis.  This  was 

first  seen  and  explained  on  the  basis  of  spin- spin  interaction  by 
55 

Bleaney  et  al. 

The  dominant  interaction  in  this  case  is  the  dipolar  field  from 
the  two  nearest  neighbors.  Since  each  spin  can  be  either  parallel  or 
anti-parallel  to  the  external  field,  four  arrangements  of  the  nearest 
neighbors  are  possible,  and  this  will  be  indicated  by  arrows  : 

I  Ho  Ul>  <1T>  'ti)  <TT) 

a)  b)  c)  d) 

At  high  temperatures,  each  of  these  arrangements  is  equally, 
probable.  The  static  dipolar  field  in  case  a)  and  case  d)  will  shift 

3 

the  line  by  +  2g^p/r  =+ 188  gauss.  Case  b)  and  c)  have  no  effect 
on  line  position,  and  the  central  peak,  therefore,  should  have  twice  the 
intensity  of  the  side  peaks  s  188  gauss  away. 
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At  low  temperatures  and  in  high  fields  that  will  line  up  the  spins 
in  the  lowest  state,  we  expect  case  a)  to  become  more  and  case  d)  to 
become  less  probable.  Professor  Bloembergen  suggested  that  we  look 
at  this  lineshape  as  an  example  of  spin  alignment. 

3.2.  Measurements  of  line  areas. 

We  measured  two  approximately  spherical  samples  of  Nd 
(C^HgSO^)^  •  9H2O  in  the  waveguide  spectrometer  at  30.00  kMc/s 
using  superheterodyne  detection.  The  crystals  were  generously  given 
to  us  by  Mr.  Rieckhoff,  Vancouver.  The  magnetic  field  was  along 
the  hexagonal  axis,  and  a  typical  change  in  lineshape  and  intensity  is 
shown  in  Fig.  6.  The  strong  line  appears  on  the  low  field  side  of  the 
unperturbed  resonance,  since  in  this  geometry  the  dipolar  field  will  aid 
the  external  field.  At  T  =0.  5°K  the  absorption  was  rather  too  strong; 
about  50  ^  of  the  incident  power,  and  the  nonlinear  response  had  to  be 
corrected  as  indicated  in  Section  2.  5  . 

The  probability  that  any  spin  is  parallel  to  the  field  is 

p^  =  l/{l+b^^)  (3-3) 

♦ 

We  tried,  without  success,  to  grow  crystals  of  Nd(C2HeSO.)>  * 
9H2O  at  room  temperature.  The  ethyl  sulfates  are  tather  unstaois, 
ana  we  got  only  the  monoclinic  NdgtSO^)^  *  8H2P.  The  spin-Haxnil- 
tonian  of  this  salt  does  not  seem  to  have  Seen  published,  although  the 
resonance  has  been  investigated  at  low  frequency.  ^  We  observed  up 
to  4  separated  lines  with  g-factors  varying  from  2. 9  to  1.  2,  but  we  did 
not  seriously  attempt  to  determine  the  principal  axis  of  the  g-tensors. 
X-ray  studies  of  this  crystal  indicate  8Nd'^^^  ions  per  unit  cell.  ' 
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RELATIVE  ABSORPTION 


Fig.  6  Resonance  Lineshape  in  Nd(C2HgSO^)2  *  9H2O  as  Function  of 
Temperature. 
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where  the  Boltzmann  factor  is  written 

AE  -g«pHQ/kT 

Here  and  in  the  following  treatment,  we  neglect  any  influence 

of  the  interaction  on  the  Boltzmann  factors.  This  corresponds  to  the 

58 

first  approximation  of  McMillan  and  Opechowski,  which  Lee  has 
found  to  agree  very  well  with  the  exact  solution  for  an  Ising  chain. 

The  probability  that  any  spin  is  anti-parallel  is 
AE  /  AE 

p^  =  b  /  (1  +b  )  .  (3-4) 

The  probability  that  the  two  neighbors  of  a  resonating  spin  are 


both  parallel  is 


P  .  .  =  =  1  /  (I  + 


(3-5) 


Similarly,  the  probabilities  for  the  other  arrangements  of  two  neigh* 


bors  are  : 


P-||  =  (P.j)  +b*V 

*’iT '  ° 


(3-6) 


(3-7) 


If  the  probabilities  (3-5) ,  (3-6),  and  (3-7)  for  the  neighbor  ar¬ 
rangements  are  multiplied  by  the  probability 


p  =  (I  -  b“=)/(i  +  b“= ) 


(3-8) 


of  resonance  absorption,  we  get  expressions  for  the  relative  intensities 
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of  the  three  lines  as  a  function  of  temperature.  This  is  plotted  in  Fig. 

7  for  AE  3  h  ‘  30  kMc/s  =  1  cm"^  ,  together  with  the  areas  of  experi¬ 
mental  lines  like  those  in  Fig.  6  measured  with  planimeter  in  arbi^ 
trary  units  and  fitted  to  the  probability  curves. 

Lines  from  two  different  samples  were  measured  and  then 
separately  fitted  as  indicated  by  different  symbols.  The  indicated  un¬ 
certainties  are  partly  given  by  the  spread  in  areas  of  different  lines 
and  partly  by  an  estimated  minimum  of  ^  20  ^  uncertainty  in  mea¬ 
sured  areas  because  of  line  overlap  and  nonlinearity  from  the  large  in¬ 
tensity  range  of  2:100  the  spectrometer  had  to  cover. 

In  spite  of  the  cited  difficulties,  the  measured  line  areas  agree 
closely  with  the  calculated  intensities,  except  that  the  weakest  line 
appears  to  be  about  10  times  too  strong  at  T  =  0.  5*^  K  . 

But  the  main  part  of  this  particular  line  can  be  explained  as 
coming  from  hyperfine  components  of  the  strong  line.  Formula  (3-2) 

with  the  given  values  of  the  constants,  predicts  hyperfine  lines  of  both 

3  5 

isotope  143  (m  =  -  -^  )  at  H  =  325  gauss  and  of  isotope  145  (m  -  ~  ) 

at  H  =  350  gauss  on  the  high  field  side  of  each  resonance  line.  These 
lines  will  overlap,  and  since  each  isotope  has  8  hyperfine  lines,  the  in¬ 
tensity  of  the  combined  line  should  be  20.  5/8  =  2.  5  ^  of  the  strong  line. 
This  is  just  the  intensity  of  the  line  we  observe  s;  350  gauss  away  from 
the  strong  line. 
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3.3.  First  moment  changes. 

In  Fig.  6  a  small  shift  in  the  position  of  each  line  is  seen  as 
we  lower  the  temperature.  The  peaks  are  also  not  quite  evenly  spaced. 
Since  each  line  looks  reasonably  symmetric  and  moments  of  each 
component  are  difficult  to  measure  because  of  the  overlap,  we  will,  in 
this  section,  use  the  peak  position  of  the  lines  and  explain  their  shift 
on  the  basis  of  the  first  moment  theory  outlined  in  Chapter  1  . 

The  measured  peak  positions,  as  a  function  of  temperature,  are 
shown  in  Fig.  8  .  The  field  was  measured  with  proton  resonance 
markers,  and  each  line  was  swept  through  a  number  of  times,  going 
alternately  up  and  down  in  field  to  average  out  possible  hysteresis  or 
recorder  lag. 

In  the  moment  theory  we  stressed  the  importance  of  differen¬ 
tiating  between  like  and  unlike  neighbors,  since  they  required  different 
truncation  of  the  interaction  Hamiltonian.  Here,  since  we  see  three 
lines, although  all  spins  are  magnetically  identical,  we  have  to  treat 
neighbors  as  like  or  \mlike  depending  upon  whether  their  respective 
instantaneous  neighbor  arrangements  are  the  same  or  not.  The  proba¬ 
bility  of  "like"  or  "unlike"  neighbor  interaction  will,  therefore^  de¬ 
pend  upon  the  temperature,  since  the  probability  for  a  given  neighbor 
arrangennent  is  a  function  of  T  . 

The  interaction  between  spin  i  and  spin  j  can  be  divided  into 


four  parts  : 
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a)  Static  interaction  from  S^j  terms  with  the  two  nearest  neigh¬ 
bors.  This  gives  the  line  splitting  and  is,  thus,  already  taken  care  of. 

b)  Spin- flip  interaction  from  terms  S^^  S_.  with  the  two  nearest 
neighbors.  This  is  only  possible  when  the  flipping  spins  conserve  dipolar 
energy. 

c)  Interaction  S  .  S  .  with  all  the  more  distant  neighbors. 

d)  Interaction  S^.  S_.  with  all  the  more  distant  neighbors  that 
are  "like.  " 

We  first  treat  the  nearest  neighbor  interactions  a)  and  b).  Whether 
or  not  the  spin-flip  interaction  b)  needs  to  be  considered  depends  on  whether 
an  ion  and  its  neighbor  are  like  or  not;  i.  e.  ,  the  nearest  neighbors  of  the  two 
interacting  spins  must  be  in  the  same  state  for  b)  to  be  included.  In  other 
words,  the  change  in  resonance  frequency  of  an  ion  as  a  result  of  nearest 
neighbor  interaction  is  a  function  of  nearest  neighbor-next  nearest  neighbor 
interactions. 

Consider  the  case  of  an  ion  contributing  to  the  low  field  line.  Speaking 
semiclassically,  its  nearest  neighbors  are  by  necessity  parallel  to  the 
external  field.  Its  next  nearest  neighbors  can  then  have  the  several  differ¬ 
ent  possibilities  illustrated  in  Fig.  9.  In  case  i) ,  both  nearest  neighbors 
are  like  the  spin  denoted  by  the  cross  and  the  resonance  frequency  of  this 

spin  is  shifted  according  to  (1-43)  with  cos  6^  s  1  by  the  amount  <  h  Ay  > 

2  2  2  "3 

~  ^  ^  ^  'ij  ^  cases  ii)  and  iii)  one  neighboring  spin  is  like  and 

the  other  unlike;  so  that  the  frequency  shift  is  (^8^^  4  J  'ij  ^  ’ 

2  2-3 

whereas  for  case  iv)  it  is  H/  P  'y  • 


I 


t  t  t  t 


Probabilities  i) 

il) 

iii) 

iv) 

( 

1  \ 

'  1  +  b^®' 

(1  +b^®)^ 

( 1  +  b^® 

\T7^I 

Fig.  9  Schematic  Representation  of  Various  Neighbor  Arrangements 
for  Low  Field  Une  in  Nd  (  C2H5SO^)3  •  9H2O  . 
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Weighting  these  four  cases  by  their  probabilities  of  occurrence,  the 
change  of  the  first  moment  of  this  low  field  line  becomes 

AE 

(1 


A  mM  ^  g  ^  2.  2^  1  A  t  2,1  2%  2b 

r.j  |^(2g,  ,__^+(2g^  t^g^  ‘rrr?®? 


V'r^>']  jttX. 


(3-9) 


Hence,  the  mean  field  shift  of  this  low  field  line  resulting  from  nearest  neigh¬ 
bor  spin- flip  interactions  b)  is 


(AH) - j  - jg- 


(3-10) 


Similar  reasoning  gives,  for  the  line  shift  of  the  upper  line  due  to  inter¬ 
actions  b) 


AE 


(AH),  =  -i±-^ 


1  +b 


and  for  the  middle  line 


1-b 


AE 


(AH)  =  - - 

g,ri/7(77b^) 


(3-11) 


(3-12) 


These  contributions  to  the  shifts  of  the  lines  are  calculated  below  for  several 
temperatures. 


I 


80 


T 

0 

0.5 

1.0 

2.0 

o 

8 

(AH)j  = 

-32 

-31 

-26 

-21 

-16  gauss 

(AH)2  = 

0 

1 

6 

11 

16  gauss 

(AH)3  = 

16 

15 

10 

5 

0  gauss 

In  case  c),  the  static  field  from  all  neighbors  except  the  two  nearest 
gives  all  three  lines  the  same  shift  given  from  (1-43) 


(AH)^  = 


1  -  b 


AE 


Afi 


1  +  b 


cos  6. .  -  1 
_ U _ 


r. . 
ij 


(3-13) 


where  the  star  on  the  siunmation  means  all  neighbors  except  the  two 
nearest  ones. 

59 

Daniels  has  given  the  total  lattice  sum  over  all  ions  closer  than 
a^^To  as  18.  7  X  a^  From  this  result  it  is  easy  to  find  the  lattice 

sum  needed  here. 
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3  cos**  -  1 


ij 


18.  7 


2x2 


13.992^-  10 


7.07-'.  10 


=  -43.  3 


10^0  cm-^ 


(3-14) 


The  shift  becomes 
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(AH). 
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1.0 
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00  ®K 

0  gauss 


I 


81 


In  case  d),  the  shift  from  spin-flip  interaction  with  all  next 
nearest  neighbors  is  proportional  to  the  probability  of  finding  the 
neighbor  in  the  same  state.  So  we  get  for  the  lower  line  shift 


Pg 


(AH), 


J-  ,1  -  b 
V’ 


AE 
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*  3  cos  0. .  -  1 

TZ —  V - SE“  ^ ^  'XE  i.  /^.  - - .(3-15) 
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middle  line  shift 
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and  upper  line  shift 
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Numerical  values  for  these  shifts  are 


T 

= 

0 

0.5 

1.0 

2.0 

4.0  oo°  K 

(AH)5 

= 

11 

9 

4 

2 

1  0  gauss 

(AH)^ 

= 

0 

1 

2 

2 

1  0  gauss 

(AH).^ 

= 

0 

0 

0 

0 

0  0  gauss 

The  various 

shifts 

1 

a-* 

O 

(AH).y 

are  added  together  and 

fitted  to  the  measured  peak  positions  in  Fig.  8  The  temperature 
shifts  and  tl^e  unequal  line  spacing  are  seen  to  be  explained  reasonably 
well  by  this  theory.  Since  the  upper  weak  line  is  mostly  hyperfine 
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components  of  the  other  lines  at  T  =  0.  5  °  K  ,  we  do  not  expect  any 
good  agreement  for  this  one. 

Two  points  about  the  line  positions  should  be  noted.  First, 
the  unperturbed  position  for  the  middle  line  seems  to  be  5960  gauss. 
For  a  frequency  of  30.  0  Mc/s,  this  corresponds  to  g  s  3.60  which  is 
somewhat  higher  than  the  reported  g  ==  3.  535  for  the  diluted  salt. 
Second,  the  splitting  caused  by  the  static  part  of  the  dipolar  interaction 
would  have  to  be  only  about  170  gauss  instead  of  the  calculated  dH 

3 

=  2g^  P  /  r.j  =  188  gauss.  The  possible  explanation  of  this  dif¬ 
ference  will  be  discussed  in  Section  3.5. 

3.4.  Second  moment  changes. 

The  second  moment  for  interaction  between  two  species  of 
spins  is  given  by  (1-45).  For  all  but  the  two  nearest  neighbors  we 
have  to  multiply  the  lattice  sums  by  the  probability  that  the  neighbors 
of  any  spin  make  it  "like"  or  "unlike"  with  respect  to  the  line  con¬ 
sidered.  For  the  lowest  line,  for  example,  this  contribution  to  the 
second  central  moment  becomes 


(h^<  6V^> 


6£«2 


“  -  -u'  2  *  ‘V  * 
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+  (1  - - TET  ^  ^  ^ 

(1  +  r  '' 


(3  C08^  0  -  1)^ 
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ij 


(3-18) 


The  spin-flip  with  the  two  nearest  neighbors  will  cause  some 

broadening  too.  As  found  before,  there  is  a  probability  p^ 

AE 

=  1/(1  ■¥  b  )  that  one  of  the  next  nearest  neighbors  will  be  arranged 

2  2  3 

so  as  to  give  the  lower  line  a  shift  of  s  =  g  p  / Z  The  resulting 

contribution  to  the  second  moment  around  the  mean  position  from  two 
neighbors  will  be 


(h*<4v^>)^  =  2[(1  -  P^)  (P^.)^  +  P|(8  -Pt*)^] 

g,  P 

,,  »  2  -L  r  1 

=  2(P|  -  P|)  8  - - TT  [ 


2  r, 
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1  +b 


AE 


- ^-tjj^].(3.  19) 

(1  +  b^^)^ 


The  total  second  moment  will  be  the  sum  of  (3-18)  and  (3-19)  ■ 

The  lattice  sum  in  (3-18)  can  be  found  by  summing  over  all 
lattice  points  inside  a  sphere  of  radius  R  and  integrating  over  an 
assumed  \miform  density  of  dipoles  outside  the  sphere.  If  the  sum.- 

o 

mation  is  taken  over  the  64  next  nearest  neighbors  within  R  =  21.  5  A, 
we  find  (please  see  following  P*go) 


84 


»  (3  co»^8..  -  1)* 
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«  (640  +  150)  •  lO'*®  =  790  •  10^°  cm“^ 


(3-20) 


Here  we  have  used  the  dipole  density 


c  = 


2  •  cnO 

Sq  Sin  60  *  Cq 


1-  A-^ 


600 


The  square  root  of  the  theoretical  total  second  moment  for  the  lower 
line  is  plotted  as  a  function  of  temperature  in  Fig.  10  for  = 

1  cm  \ 

Since  interaction  between  like  neighbors  broadens  a  line  more 
than  interaction  between  unlike  spins,  and  since  the  probability  for 
finding  like  neighbors  increases  for  this  line  as  we  lower  the  tempera¬ 
ture,  we  might  expect  the  second  moment  to  have  a  maximum  at  some 
low  temperature.  This  would  indeed  happen  if  the  ratio 
been  larger. 

In  the  same  way,  the  total  second  moment  for  the  middle  line 
becomes  : 
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ae  .4  r2i>^^ug»^  +  g,¥ 
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And  for  the  upper  line  we  find 
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(3-22) 


The  square  roots  of  these  second  moments  are  shown  as  functions 
of  the  temperature  in  Fig.  10  .  The  small  differences  in  second  mo¬ 
ment  between  the  three  lines  reflects  the  different  probability  of  spin- 
flip  broadening. 

It  is  difficult  to  compare  the  calculated  moments  to  experimen¬ 
tal  data  since  the  lines  are  partially  overlapping.  But  one  can  measure 
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second  central  moments  reasonably  well  on  the  side  of  a  line  where  it 
does  not  overlap.  This  can  always  be  done  for  the  lower  line,  and 
some  measured  second  moments  for  this  line  are  given  in  Fig.  10  . 
The  upper  line  vanishes  at  sufficiently  low  temperature,  and  then 
measurements  of  the  middle  line  can  be  done  also. 

The  agreement  is  good  above  2°K  ,  but  at  lower  temperatures 
the  lines  appear  too  broad.  This  is  probably  caused  by  variations  in 
local  field  over  the  sample  since  the  spheres  used  were  rather  crude. 

In  the  previous  section  we  saw  that  the  distant  neighbors  would 
2  2 

shift  the  lower  line  g^  /2g^  more  than  the  other  two  lines  at 
T  =  0°  K  because  of  the  possibility  for  spin-flip-caused  shift.  Like¬ 
wise,  variations  in  the  local  demagnetizing  field  will  broaden  the 
lower  line  more  than  the  others  by  the  same  amount.  Together  with 
the  extra  broadening  of  the  lower  line  shown  in  Fig.  10  ,  this  should 
account  for  the  measured  difference  in  second  moment  between  the 
lower  and  middle  line.  Qualitatively  one  sees  this  difference  just  by 
looking  at  the  lineshapes  in  Fig.  6  . 


3.  5.  Discussion  of  possible  antiferromagnetic  coupling  in  neodymium 
ethyl  sulfate. 

As  we  saw  in  the  previous  sections,  the  resonance  line  changes 
with  temperature  in  neodymium  ethyl  sulfate  are  e3q>lained  fairly  well  on 
the  basis  of  dipolar  interaction.  But  to  fit  the  data,  the  line  splitting 
from  the  static  interaction  with  the  two  nearest  neighbors  had  to  be 
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taken  as  about  170  gauss  instead  of  the  calculated  188  gauss.  Some  of 
this  difference  can,  perhaps,  be  accounted  for  as  systematic  error, 
since  we  used  the  peak  positions  instead  of  the  first  moments  of  slight¬ 
ly  unsymmetric  lines.  It  does  not  seem  likely,  however,  that  this 
error  would  amount  to  18  gauss,  and  it  is  interesting  to  postulate 
some  exchange  interaction  to  explain  this  shift. 

The  exchange  would  have  to  be  only  to  the  two  nearest  neighbors 
of  each  spin  since  the  upper  and  the  lower  lines  have  to  be  shifted  in  op¬ 
posite  directions.  It  would  have  to  be  antiferromagnetic  in  character 
since  its  effect  is  to  reduce  the  dipolar  splitting  in  this  geometry.  But 
it  can  hardly  be  isotropic  since  the  resulting  texnperature -dependent 
shift  from  spin-flip  of  varying  probability  will  tend  to  destroy  the  good 
agreement  between  calculated  and  measured  temperature  change  of 
line  positions. 


We  are  then  left  with  the  possibility  of  the  anisotropic  exchange 


3C..  =  A. .  S  .  S  . 

ij  zi  zj 


-  [a 


for  near  neighbors 


otherwise 


(3-23) 


I  \M 

which  for  A  *  6  *  10  erg  would  give  the  required  shift  of 


I  A  H  I  =  2  •  A/g^P  =  18  gauss 


This  is  a  very  small  exchange  and  alone  corresponds  roughly  to 
_ _ _ ,  60 


a  Curie  temperature  of 
0  »  - 


'Xs  (S  <(•  1)  Z 


0.0018”K 


(3-24) 


when  Z  s  2  is  the  number  of  nearest  neighbors. 
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Since  the  crystal  is  not  cubic,  the  dipolar  interaction  also  can 

give  a  spontaneous  magnetization.  From  the  crystal  structure  and 

59 

the  g-values  for  dilute  salt,  Daniels  calculated  a  powder  Curie 
temperature  of  0.  007^  while  Roberts,  Sartain  and  Borie^^  measured 
Opowder  ~  difference  is  probably  not  signifi¬ 

cant,  but  if  it  were,  it  predicts  that  the  additional  coupling  should  be 
ferromagnetic  rather  than  antiferromagnetic. 

The  specific  heat  data  allows  a  somewhat  better  comparison 

between  theory  and  experiment  for  lattice  interaction.  Daniels  calcu- 

59 

lates  the  specific  heat  from  dipolar  interactions  as 


j 


+  9(g 


2  2.2 

•‘i’ 


•  8.95  •  10'® 


(3-25) 


This  is  in  excellent  agreement  with  die  esqierimental  value  of 
-5  61 

8.  85  ^  0.  25  *  10  ,  and  does  not  seem  to  leave  room  for  any  ex¬ 

change. 

Equation  (3-25)  is  derived  without  consideration  as  to  whether 
spin  neighbors  are  "like"  or  "unlike"  as  we  had  to  do  in  the 
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calculation  of  moments.  This  is  valid  when  no  external  field  is 
present  to  quantize  the  spins  in  the  z-direction.  But  (3-25)  involves 
%!!  to  the  4th  power  and  r.j  to  the  6th  power,  and  small  errors  in 
these  constants  can  easily  upset  the  apparent  agreement.  The  values 
of  r.j  are  from  room' temperature  measurements;  the  coefficient  of 
thernnal  expansion  is  unknown  for  neodymium  ethyl  sulfate,  but  it  is 
quite  possible  that  the  dimensions  could  decrease  by  a  factor  of,  say, 

3  *  10  from  room  temperature  to  helium  tempieratures  as  is  the  case 
for  many  materials.  That  would  give  a  2  increase  in  CT^/2R  .  And 
even  more  important,  we  have  measured  g^^  =  3.60  in  the  concentra¬ 
ted  salt,  and  the  use  of  this  value  in  (3-25)  instead  of  fLji  -  3.  535  in- 
creases  CT  /2R  by  about  6  The  calculated  specific  heat  is  then 
too  large,  but  when  the  small  antiferromagnetic  coupling  (3-23)  is  in¬ 
cluded,  it  is  reduced  to  agree  again  with  the  experimental  data. 

The  easiest  way  to  include  the  effect  of  the  exchange  in  die 
specific  heat  is  to  change  g^  in  (3-25)  when  interaction  between 
nearest  neighbors  is  considered.  The  coupling  (3-23)  added  to  the  di¬ 
polar  S  .  S  .  -term  in  the  interaction  between  nearest  neighbors 

Z1  ZJ 

shifts  the  line  by  «  18  gauss  out  of  a  total  of  188  gauss,  and  the  inter¬ 
action  energy  is  reduced  to 


_  ,170, 

*^ij  ■  ^188^  “5 
*  • 


'  jvns 


3.42  . 


(3-26) 
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This  apparent  g»  should  be  used  in  (3-25)  whenever  the  sunns  ^ 

"  j 

involve  one  of  the  nearest  neighbors.  We  use  the  total  lattice  sums 
as  given  by  Daniels  for  neodynniunn  ethyl  sulfate  and  split  them  up  in 
the  contribution  from  the  two  nearest  neighbors  and  all  the  rest  as 
follows; 

Daniels  Two  nearest  The  rest 


4^  TT  13.99^  •  10"“*® 
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=  (1.59  +  1.66)  10 
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r. . 
ij 
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=  (1.59  +  0.32)  10 
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*ij 


■ij 


128. 07 
13  99^  10 


=  (1.59  +  0.11)  10 
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We  then  find  CT^/2R  »  8.  35  *  10~^  ,  and  when  this  is  increased  by 
fa  2  ^  because  of  the  decrease  in  the  interionic  distances,  we  get 


I  2 

S  T^ 


»  8. 55  •  10 


-5 


(3-27) 


This  is  almost  within  the  accuracy  of  the  eiqperimental  result 
CT^/2R  =  8. 85  4  0.  25  *  10’^  .  But  this  agreennent  with  specific  heat 
data  should  not  be  taken  as  a  proof  for  the  existence  of  the  coupling 
given  by  (3-23)  .  In  view  of  the  uncertainties  in  g-value  and  ionic 
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spacings,  probably  all  that  can  be  said  is  that  the  specific  heat  data 

certainly  does  not  rule  out  the  possibility  of  such  an  interaction. 

Some  weak  interaction, in  addition  to  the  dipolar  forces,  should 

not  be  too  unexpected  in  neodymium  ethyl  sulfate.  Demagnetization 

e:q>eriments^^  on  the  isomorphic  cerium  ethyl  sulfate  give  a  specific 

heat  for  this  salt  that  is  six  times  larger  than  what  we  expect  from  the 

dipolar  interaction  alone.  The  Curie  temperature  in  this  salt  is  of  the 

order  -0.  l^K  which  corresponds  to  an  antiferromagnetic  exchange 

m  50  times  stronger  than  the  interaction  (3-23)  . 
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Finkelstein  and  Mencher  have  suggested  that  the  anomalous 

behavior  of  cerium  ethyl  sulfate  may  be  due  to  interactions  between  the 

atomic  electric  quadrupole  moments  of  the  magnetic  ions.  Such  an 

interaction  would  depend  upon  the  relative  orientation  of  the  spins  and 

give  terms  in  S  .  S  .  and  S ..  S  .  .  But  the  relative  strength  of 

these  terms  is  likely  to  depend  on  the  anisotropy  of  the  g -tensor,  and 

we  expect  the  spin-flip  term  to  be  smaller  than  the  static  term  by  a 

factor  g  l%ii  sO.33.  Our  experimental  accuracy  is  not  sufficient 

to  decide  whether  a  term  like  A..  (S^  S  .  +  S  .  S,.  )  is  needed  in 

6  ij  +4  -J  "1  +J 

(3-23)  or  not.  The  quadrupole  interaction  would  probably  be  important 
only  for  the  nearest  neighbors  since  it  falls  oU.  as  the  5th  power  of  the 

t 

interatomic  distance. 

Comparing  cerium  and  neodymium,  the  quadrupole  interaction 
would  be  proportional  to  the  8th  power  of  the  4f  electron  orbit 
radius  and  proportional  to  the  matrix  elements  connecting  the  ground 


92 


state  with  the  lowest  excited  level  of  the  ion.  Cerium  is  the  first 
member  of  the  rare  earth  series,  and  we  expect  the  mean  radius 

for  its  4f  electron  to  be  considerably  larger  than  in  neodymium. 

3+  -1 

Also,  the  first  excited  state  of  Ce  in  ethyl  sulfate  is  only  4.8  cm 

62  3  + 

above  the  ground  level,  while  in  Nd  the  splitting  is  about 
-1  53 

130  cm  The  electronic  quadrupole  interaction  in  neodymium 

ethyl  sulfate  should,  therefore,  be  much  smaller  than  in  the  corres> 
ponding  cerium  salt,  and  it  could  well  be  smaller  by  the  two  orders  of 
magnitude  that  are  required  for  the  quadrupole  interaction  to  agree  with 
our  postulated  exchange.  It  is,  therefore,  not  unreasonable  to  suggest 
that  this  may  be  the  physical  origin  of  the  weak  antiferromagnetic 
interaction  that  we  have  observed  in  neodymium  ethyl  sulfate. 

Since  this  interaction  is  mainly  static  in  character  and  between 
nearest  neighbors  only,  it  will  not  influence  the  second  moments  of  the 
resonance  lines. 


CHAPTER  4 


i 

1 

MEASUREMENTS  ON  NICKEL  FLUOSILICATE  I 

4.1.  Description  of  nickel  fluosilicate. 

Nickel  fluosilicate,  NiSiF^  •  belongs  to  a  large  class  of 

24  44  ~ 

isomorphic  crystals  of  formula  M  N  A^  where  M  may 

44 

be  Ni  ,  Cu  ,  Fe  ,  Co  ,  Mn  ,  Zn  ,  Mg  or  Cd  ,  N  may  be  Si  ,  Ti  , 

i 

Sn  ,  Zr  ,  Pd  or  Pt  ,  and  A^”  may  be  H  ,  F  or  Cl  .  The  structure 
of  Ni  SnClg  ‘  has  been  determined  from  x-ray  measurements  by 

Pauling.  It  is  trigonal  and  the  unit  cell  containing  one  molecule  is 
shown  in  Fig.  11  .  The  structure  may  be  thought  of  as  made  up  of 
(Ni  bH^O^^and  (SnCl^)^"  complexes  arranged  ona  body-centered  cubic 

! 

lattice  that  has  been  trigonally  distorted  along  the  body  diagosnal  denoted 

by  z  in  Fig.  11  .  This  body  diagonal  is  the  same  as  the  hexagonal 

growth  axis.  Each  complex  is  assumed  to  have  the  form  of  an  almost 

regular  octahedron  surrounding  the  metallic  ion. 

The  other  members  of  the  crystal  series  are  assumed  to  have 

the  same  structure  with  only  small  changes  in  the  dimensions  of  the 

unit  cell.  For  the  case  of  NiSiF^  •  6H2O  ,  the  Ni  -  Ni  distance  OA 
65  Q 

in  Fig.  11  is  given  as  6.  21  A  and  the  angles  a  between  the  rhom- 
bohedral  axes  are  96^  5'  .  This  is  close  enough  to  cubic  symmetry 
(a  =  90®  )  so  that  the  lattice  sum  ^  (3  cos^  Ojj  -  l)/rjj^  is  found 

J  CC 

to  be  essentially  the  same  as  in  the  cubic  case. 
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Ni  Si  Fg*  6H2O 

Fig.  11  Rhombohedral Unit  Cell  of  NiSiF^  •  6H2O  . 


I 


94 


67 

This  is  the  structure  at  room  temperature.  Walsh  has 
measured  the  thermal  expansion  Of  this  crystal  with  strain  gauges.  He 
found  the  expansion  coefficients  to  be  highly  anisotropic,  and  from  his 
results  we  expect  the  deviation  from  cubic  symmetry  to  be  slightly 
larger  at  low  temperatures.  But  the  lattice  sum  can  probably  still  be 
assumed  to  vanish  for  spherical  samples. 

From  the  structure  we  expect  the  Ni^^  -  ion  to  see  a  slightly 
trigonally  distorted  cubic  electric  field  from  the  surrounding  water 

3 

dipoles.  The  cubic  field  splits  the  ground  state  of  the  free  ion 

and  leaves  an  orbital  singlet  level  some  10,000  cm~^  below  the  other 
68 

states.  The  triple  spin  degeneracy  of  this  level  is  lifted  by  the  trig¬ 
onal  field  vdiich  splits  it  into  a  doublet  and  a  singlet.  The  xxiagnetic 
behavior  can  then  be  described  by  the  following  spin -Hamiltonian  for 
S  =  l 

Hj*®*  =  +D(S^*  -  )  .  (4-1) 

No  hyperfine  interaction  term  needs  to  be  included  here  since 
the  Ni-isotope  with  nuclear  spin  has  low  abundance  aiul  small  nuclear 
moment. 

If  the  external  field  makes  an  angle  Y  with  the  x-axis  , 

and  if  we  assume  g  =  g  =  g.  =  g  >  we  find  the  energy  levels  of  (4-1) 

X  y  z 

by  firat-order  perturbation  calculation  (exact  for  y  ■  0°)  . 
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*1  =  ^D(|co8^Y  -j) 

a  =  iD  (1  -  3  coB^  Y  )  ^ 

^3  =  gpHQ+yD  (|co8^Y  --j) 

The  re8onance  transitions  are 

hvi  =  ^2  -  =  gpHjj  -  -g  D  (3  cos^  Y-l) 

hv2  =  ^3  ~  ~  *P^0  ■*'1^  co8^Y-l) 


(4-2) 


(4-3) 


For  Y  4  ^  crystalline  field  mixes  the  states  and  the 

transition  a^  -  a^^  also  becomes  allowed,  but  with  small  probability 
when  gpHQ  »  D  . 

We  have  measured  at  helium  temperatures  g  -  2.  24  ^  0. 01  , 

isotropic  within  the  accuracy  of  the  measurements  This  is  in  good 

agreement  with  the  g  =  2.  25  deduced  from  optical  rotation  data  at 

low  temperature^^  and  does  not  differ  significantly  from  the  g  w  2.  25 

70  71 

-  2.  36  reported  in  resonance  experiments  at  higher  temperatures.  .  ’ 

We  have  found  the  splitting  D  =  -  0. 113  +  0.  005  cm~^  at 

helium  temperatures.  Above  20^  K  it  is  a  function  of  temperature 

-1  70  71 

and  increases  to  D  »  -  0.  50  cm  at  room  temperature.  ’  Hy- 

67 

drostatic  pressure  will  decrease  |D  |  again,  and  Walsh  has  shown 
that  at  room  temperature  D  actually  goes  through  sero  and  changes 
sign  for  p  =  6200  kg/cm^  .  Since  the  splitting  is  so  sensitive  to  lat¬ 
tice  distortions,  it  is  not  unreasonable  to  eiqpect  that  crystalline 
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defects,  impurities  and  cooling  strains  will  give  a  slight  variation  in 

the  value  of  O  from  ion  site  to  ion  site.  This  view  is  supported  by 

the  fact  that  the  narrowest  resonance  lines  measured  in  Zn-diluted 

72 

salts  are  about  30  gauss  wide,  much  more  than  one  should  expect 

from  dipolar  interactions  with  the  surrounding  protons. 

Earlier  measurements  of  pacific  heat,  magnetlBation  and 

resonance  linewidths  of  NiSiF^  *  6H2O  have  indicated  the  presence  of 

moderately  strong  exchange  in  this  salt.  But  there  has  been  some 

confusion  as  to  the  magnitude  of  the  exchange  and  whether  it  would  have 

to  be  partly  anisotropic.  As  will  be  seen,  our  measurements  indicate 

-18 

anisotropic  ferromagnetic  exchange  A  s  -  3.9  *  10  erg  in  good 
agreement  with  some  of  the  previous  measurements. 

Our  lineshape  data  was  taken  mainly  with  the  external  fieldalong 
the  trigonal  axis  (T  s  0) .  Thenwe  see  two  exchange -shifted  and 
broadened  lines  separated  at  high  temperature  by  2|D|/gp.  However, 
we  see  from  (4-3)  that  the  two  lines  overlap  ior  the  Bold  direction  Y  *^55° 
where  (3  cos^  y  -  1)«0  .  Here  the  combined  line  is  exchange -narrowed 
and  unshifted.  But  all  directions  ywSS^do  not  give  the  same  linewidth. 

In  particular,  the  lattice  sum  ^(3  cos  -  1)  /r^^  determining  the 
second  moment  is  very  different  for  parallel  to  the  direction  7  in  Fig. 


11,  udiich  happens  to  be  almost  along  one  of  the  rhombohedral  axes,  and  for 
^  I)  which  is  obtained  by  rotating  an  angle  2  X  55**  in  the 
plane  AZB  .  This  difference  in  linewidth  is  interesting  to  observe, 
but  is  no  more  surprising  then  the  variation  in  the  square  root  of  the 
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second  moment  1  :  1.  5  :  2.4  found  for  the  directions  (111)  ,  (110)  and 

73 

(100) ,  respectively,  in  a  cubic  crystal. 

4.  2.  First  moment  changes. 

We  tried  to  grind  spheres  of  NiSiF^  *  ,  but  we  did  not 

succeed,  since  the  crystals  are  soft  and  cleave  very  easily  along  the 
hexagonal  axis.  Instead,  they  came  out  as  reasonably  good  ellipsoids 
of  revolution  with  axial  ratios  ranging  from  1.  25  to  1.  50  .  The  field 
for  resonance  then  had  to  be  corrected  for  demagnetization  effect  as 
discussed  in  Section  1.  5  .  For  an  axial  ratio  of  1.'50  and  a  spin  den¬ 
sity  of  4. 16  •  10^^  ions/cm^  ,  we  find,from  (1-49)  and  (1-57)  ,  a  shape 
correction  of  164  gauss  at  T  s  0°  K  which  is  far  from  negligible. 

The  ellipticity  of  the  samples  had  one  beneficial  effect,  how¬ 
ever;  they  were  much  easier  to  x-ray  orient  since  we  could  see 
roughly  the  direction  of  the  crystalline  axis. 

We  measured  four  different  samples  in  the  cavity  version  of 
the  spectrometer  and  two  samples  in  the  waveguide  version.  The  fre¬ 
quencies  used  were  slightly  different;  f  =  29.  80  kMc/s  for  the  cavity 
and  f  s  30.06  kMc/s  for  the  waveguide  spectrometer.  For  compari¬ 
son  of  the  results  we  reduced  the  measured  field  values  for  waveguide 
resonance  by  (29.80/30.06)  .  Typical  resonance  lineshapes  for  Hq 
along  the  hexagonal  z-axis  are  shown  in  Fig.  12  ynhere  the  maximum 
amplitudes  are  normalized  to  the  same  value.  In  this  plot  the  field 
values  are  not  corrected  for  the  effect  of  a  nonspherical  shape,  but 
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the  lines  at  0.4**  K  and  0.93**  K  are  corrected  for  too  much  absorption 

by  using  formula  (2-31)  .  Completely  off  resonance  we  had 

M  0.  05 ;  at  resonance  at  T=0.4**K  we  found  P  /P  w0.70. 

r  o 

From  Fig.  12  we  see  that  at  4.  2**  K  we  have  two  resonances 
of  almost  equal  intensity  as  expected.  Since  |D|  is  small,  the  line  - 
splitting  is  not  as  large  as  we  might  have  liked  to  have  it,  and  the  lines 
overlap  somewhat.  As  we  lower  the  temperature,  the  transition  a^ 

-  a^  at  the  higher  field  decreases  in  strength  relative  to  the  other 
transition,  and  it  also  changes  position.  Below  1**  K  it  cannot  be 
separated  from  the  wing  of  the  line  a^  -  a^^  in  the  lower  field.  The 
lower  line  narrows  considerably,  and  also  moves  somewhat.  An  ad¬ 
ditional  weak  line  is  seen  right  between  the  main  peaks;  the  origin  of 
this  one  is  not  known,  but  some  possible  explanations  will  be  discussed 
in  Section  4.  4  . 

Because  of  the  line  overlap,  it  is  difficult  to  measure  first 
moments,  and,  therefore,  we  have  again  used  the  peak  positions  of  the 
lines.  The  peak  positions  of  the  lower,  strong  line  shown  in  Fig.  13a, 
are  corrected  for  the  demagnetizing  shift  according  to  (1-49)  and 
(1-57)  .  Since  this  salt  is  almost  cubic,  the  dipolar  parts  of  the  lattice 
sum  cancel  and  do  not  contribute  to  the  first  moment.  However,  as 
we  saw  in  Section  1. 3  ,  isotropic  exchange  will  shift  the  first  moment, 
of  this  line,  and  from  (1-36)  we  find 


h  <6v  >  *  “  6A 

(b 


"  -  b^) 


(4-4) 


I 


MAGNETIC  FIELD  IN  GAUSS 


TEMPERATURE 

Fig.  13  Measured  Peak  Positions  of  Resonance  Lines  in  NiSiF^  •  6H2O 
Compared  to  Calculated  First  Moments  • 
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Here  we  have  assumed  the  exchange  to  have  appreciable  strength 

only  to  the  six  nearest  neighbors  of  each  spin.  The  best  fit  to  the  ex- 

perimental  values  in  Fig.  13a  is  seen  to  be  for  the  choice  of  A 
-18 

=  -  3.9  '  10  erg  .  This  is  a  ferromagnetic  interaction.  Antiferro¬ 
magnetic  exchange,  A  >  0  ,  would  have  shifted  the  line  the  other  way. 

Similarly,  the  demagnetization-corrected  peak  position  of  the 
upper,  weak  line  is  plotted  in  Fig.  13b  .  Again,  dipolar  interaction 
does  not  change  the  first  moment,  and  the  exchange  shifts  it  by 


h  <Av>  = 


^  2k*  21^ 

6A  (b  -  b  ^  ) 

^1  *2  *3 

(b  Ub  '‘+b  ■*) 


(4-5) 


The  choice  of  A  =  -  3.9 


10"^®  erg 


fits  the  measured  values 


reasonably  well. 

The  different  shifts  of  the  two  lines  can  be  esqplained  as  follows: 


In  the  lower  line,  the  exchange  term  S^.  will  cause  a  shift  propor- 

*1  *3 

tional  to  (b  -  b  ).  The  spin-flip  term  S  .  causes  a  shift  in 

*1  ^2.'^ 

the  opposite  direction  proportional  to  (b  -  b  ;  .  Since  the  propor- 

*2  *3 

tionality  constants  are  the  same,  the  net  shift  goes  as  (b  -  b  )  , 
which  is  zero  for  T  — >  oo  and  T  — ^  0  .  For  the  upper  line,  how¬ 


ever,  spin-flips  with  near  neighbors  become  less  probable  as  the  two 

*2  *3 

upper  levels  are  depopulated,  and  the  shift  proportional  to  (b  -b  ) 

*1  *3 

is  insufficient  to  cancel  the  shift  (b  ~  b  )  caused  by  the  S^|  S^ 
terms. 

Qualitatively,  we  may  think  of  a  spin  resonating  in  the  lower 
transition  as  having  a  fictitious  spin  ^  .  This  spin  will  see  only  like 
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neighbors  at  temperatures  so  low  that  the  uppermost  level  is  ef  • 
fectively  depopulated,  and  isotropic  exchange  cannot  then  shift  the  line. 
A  spin  resonating  in  the  upper  transition,  however,  can  bethought  of  ds 
seeing  mainly  unlike  neighbors  at  low  temperatures,  and  exchange  will 
then  change  the  first  moment. 

Anisotropic  exchange  would  have  given  a  very  different  first 
moment  shift.  If  we  were  to  assume  anisotropic  exchange  of  the  form 


.  (1)  _ 
u 


(4-t) 


where  we  write 


tan  ^ 


(4-7) 
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We  see  from  the  graphs  of  McMillan  and  Opechowmki  that  our 
measurements  will  only  fit  'Q-  =  45°  +  5°  which  corresponds  to 

A. .  B  A. .  to  within  4  20  4 . 

We  also  looked  for  possible  temperature  shift  in  line  position 
when  made  an  angle  of  cos^  0  =1/3  with  the  hexagonal  axis 
and  the  two  lines  overlapped.  For  this  direction,  the  demagnetisation 
correction  should  vanish  in  ttie  first  order.  We  saw  no  Uneshift 
larger  than  the  uncertainty  in  measurements,  and  this  may  be  taken  as 
additional  indication  of  the  absence  of  anisotropic  exchange. 

Anisotropic  exchange  having  the  same  angular  dependence  as 
dipolar  interaction  could  still,  in  principle,  be  present.  This  would 
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not  affect  first  moments  ,  since  the  crystal  is  almost  cubic  in 
form. 


4.3.  Second  moment  chances. 


The  second  moments  of  the  resonance  lines  in  NiSiF^  ■  6H2.O 

75 

have  been  calculated  at  infinite  temperature  by  Ollom  and  by 
Ishiguro,  Usui  and  Kambe.  McMillan  and  Opechowski^  used  this 
salt  as  an  example  for  their  formulas  for  the  temperature  dependence 
of  the  moments,  and  from  (A-10),  they  find,  for  the  transition  a^  -  a^ 
when  Hq  is  along  the  crystalline  axis 

h^<Av^>  =  6A  [Z  Mb  +  2b  )  +  5  (b  ^  -  b  ] 


-4  4  1 

+  0. 626  P-A  ■■  <  1  +  4Z  (2b  +  b  ) 

d° 


2  2a.  a,  a.  ^^2  ^2 

Z"^  (9b  -  6b  ‘  ^  +  b  ^  -  12b  ^  ^ 


2a.  ^2^4 

+  4b  ^  +  4b  ^  ^  ) 


(4-8) 


*1  ®2  *3 

Here  we  have  used  g^=g^  =  giZ=(b  fb  fb  ),  and  the 
result  of  the  lattice  summation 


^  ^  (3  cos^  0.j  -  =  0. 626/d^  , . 


(4-9) 


where  the  lattice  constant  is  d  =  6.  21  A  . 
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The  second  central  moment  can  he  found  by  combining  (4-8) 
and  (4-4) 

^  Central  -  <iv*  >  -  (h  <  4  V  > 

=  6A^  [Z‘^  (b*^  ♦  2b*^)  -  Z"^  (b*^  -  b*^)^  ] 


4^  4^  I  &  St 

•f  0.626  fiX  I  1+4Z  ■^Zb  ^  +  b  ^  ) 


-2  ^1  *1*2  ^2  *1*3 

-  Z  ^  (9b  -  6b  ^  ^  +b  ^  -  12b  ^  ^ 


Za^  *2*3 

♦  4b  ^  +  4b  ^  ^  ) 


(4-10) 


The  square  root  of  (4-10)  is  plotted  in  Fig.  14  as  a  function  of 

~  -18 

temperature  for  V  s  29.80  kMc/s  and  A  -  -  3.9  *  10  erg  .  This 
is  compared  to  measured  values  of  the  square  root  of  the  second  cen¬ 
tral  moment  taken  on  the  non-overlappixig  low  side  of  the  lower  line. 
The  agreement  is  seen  to  be  good.  Since  about  90  ^  of  the  second 
moment  comes  from  the  exchange  term,  this  is  a  good  check  on  the 
magnitude  of  A  found  from  the  first  moment  shifts,aiid  the  presence  of 
any  anisotropic  exchange  must  be  ruled  out. 

It  is  seen  that  even  if  the  line  clearly  narrows  at  very  low  tem- 


I 


peratures,  the  second  moment  remains  somewhat  larger  than  predicted 
by  the  theory.  This  is  probably  caused  by  imperfect  elliptical  shape 
of  the  samples,  and  possibly  by  nonlinearity  from  too  nnuch  sample 
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saturation  or  sample  heating  at  0.4°  K  .  A  perfect  spheroidal  sample 
shape  will  not  influence  the  second  central  moment  measured  in  this 
orientation. 

A  similar  formula  could  be  given  for  the  second  moment  of  the 
other  line,  but  this  weak  line  disappears  below  1°  K  ,  and  the  formula 
cannot  very  well  be  compared  to  experimental  results. 

At  directions  where  the  two  lines  overlap  we  find  the  second 
moment  from  (1-25)  with  S  =  1 

h^<Av^>  =  f(l+3f)  2  •  <4-11) 

j 

where 

'  ■  ViTT-Wnr- 

b  b 

— >  j  for  T  —4  00  (4-12) 

In  using  formula  (1-25)  for  this  case,  we  have  neglected  small 

2  2 

terms  in  the  matrix  elements  between  the  levels  of  order  D  / (gPHQ)  . 

When  the  field  is  parallel  to  one  of  the  rhombohedral  axes  (di¬ 
rection  in  Fig.  11  )  the  lattice  sum  over  92  nearest  neighbors  is 
found  to  be 

2  (3cos^0.j  -  1)^ =  224-  10^^  cm"^  .  (4-13) 

j 
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The  square  root  of  the  second  moment  for  this  direction  of  the 


field  is 


8  P  { 


1 

(3cos^  0.,  -  1)^-’^ 

- r - 

r.. 


-] 


•4  380  gauss  for  T  —4  oo 


(4-14) 


At  the  other  resonance  crossover  point  u^en  Hq  is  in 
direction  of  Fig.  11  the  lattice  sum  is 


2 


(3  cos^  0.j  - 


57. 10^^  cm"^ 


(4-15) 


and  the  square  root  of  the  second  moment  is 


1 
,Z 


-|1 

(  <  AH-  >  )-  =  j^f  (1  -  3f)  •  5. 60  •  10^  gauss 


193  gauss  for  T  — 4  oo 


(4-16) 


Formulas  (4- 14)  and  (4- 16)  as  functions  of  temperature  are 
shown  in  Fig.  15  together  with  some  measured  values  of  the  square 
roots  of  the  second  moments.  The  agreement  is  fair.  But  most  of  our 
measured  values  are  too  small,  and  this  shows  how  easy  it  is  to  under¬ 
estimate  the  intensity  in  the  wings  of  an  exchange -narrowed  line.  For 


■Q  1 1  b,  however,  the  measured  width  is  much  larger  than  predicted  at 
77^.  This  is  probably  caused  by  beginning  spin- lattice  relaxation  broad¬ 


ening  .  At  room  temperature  we  find  all  line  s  to  be  significantly  broadened. 
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Below  1°K  the  sample  magnetization  makes  the  moment  calcula¬ 
tion  uncertain  since  the  field  is  not  along  the  spheroid  axis.  We  are 
also  likely  to  have  experimental  difficulties  here  because  of  too  much 
signal,  and  therefore,  no  second  moment  measurements  were  made  at 
the  lowest  temperature. 

4.4.  Discussion  of  the  nickel  fluosilicate  results. 

Our  measurements  of  resonance  line  moments  are  in  good 

agreement  with  the  theory  if  we  assume  an  isotropic,  ferromagnetic 

~  -18 

interaction  of  strength  A.j  =  -(3.  9  4  0.4)  *  10  erg  to  the  six  nearest 
neighbors  of  each  ion.  We  will  now  see  how  this  compares  with  the  re¬ 
sults  from  other  investigations. 

The  specific  heat  of  NiSiF^  '  6H2O  has  been  measured  by 
76 

Benzie  and  Cooke,  and  the  susceptibility  has  been  found  from  the  op¬ 
tical  Faraday  effect  by  Bequerel  and  Opechowski.  These  results 

75  77 

have  been  reanalyzed  by  Ollom  and  Van  Vleck  '  using  the  crystal¬ 
line  field  splitting  constant  D  given  by  paramagnetic  resonance  ex¬ 
periments.  They  found  good  agreement  by  assuming  some  exchange 
that  would  give  a  Curie  temperature  of  0  =  0.1^.  From  (3-24)  this 
corresponds  approximately  to 

^  “  -3.5.10"^®  erg 

when  we  take  the  number  of  nearest  neighbors  to  be  z  =  6  .  This  is  in 
close  agreement  with  our  result  for  A  . 
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But  earlier  determinations  of  exchange  strength  from  line  mo> 

ments  served  only  to  confuse  the  situation  because  of  inadequate 

measurements.  At  room  temperature  we  find  the  lines  to  be  about 

1800  gauss  wide  for  Hq  along  the  crystalline  axis,  but  at  77^K  they 

are  only  about  800  gauss  wide.  The  resonance  is  obviously  broadened 

by  fast  spin-lattice  relaxation  at  room  temperature,  and  one  cannot 

71 

use  the  lineshapes  published  by  Holden,  Kittel  and  Yager  to  deter- 

66  75 

mine  the  exchange  as  Ishiguro  et  al.  and  Ollom  have  done. 

Trying  to  reconcile  the  linewidths  with  the  magnetization  and 

77  75 

specific  heat  data,  Ollom  and  Van  Vleck  !  introduced  the  possibili¬ 
ty  of  anisotropic  exchange.  This  would  supposedly  broaden  the  lines 
without  giving  much  contribution  to  susceptibility  and  specific  heat.  In 
fairness,  however,  it  should  be  said  that  they  were  not  very  happy  with 
this  explanation,  and  in  the  end  they  rather  suspected  the  linewidth 
data. 

Other  support  for  the  idea  of  anisotropic  exchange  seemed  to 

67 

come  from  the  pressure  experiments  at  room  temperature  by  Walsh. 
He  expected  to  see  the  resonance  line  exchange -narrowed  when  the 
hydrostatic  pressure  caused  the  crystalline  field  splitting  D  to 
vanish.  He  attributed  the  failure  to  observe  this  to  anisotropic  ex¬ 
change,  but  we  have  demonstrated  that  lifetime  broadening  at  room 
temperature  is  the  real  cause  and  that  the  exchange  is  almost  com¬ 
pletely  isotropic. 
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In  Fig.  12,  showing  the  shape  of  the  resonances,  there  appears 

a  weak  line  between  the  two  main  peaks.  This  line  was  also  observed 
78 

by  Hempstead.  It  could  possibly  be  caused  by  two-photon  transitions, 
but  we  found  its  relative  intensity  not  to  vary  when  we  changed  the 
microwave  power  level. 

It  is  more  likely  that  it  is  caused  by  some  variation  in  the  ex¬ 
change  and  crystalline  splitting  parameters  so  that  for  a  few  percent 
of  the  spins  |D|  <  |A|.  Such  variations  can  be  caused  by  imper¬ 
fections  in  the  crystal  lattice,  by  strain  around  cracks  and  impurities, 
etc.  The  resonance  line  for  these  ions  will  appear  exchange -narrowed, 
and  it  will  fall  just  between  the  two  main  lines. 

From  the  relatively  narrow  lines  observed  in  Zn-diluted 
72 

NiSiF^  •  6H2O  ,  it  does  not  appear  possible  that  the  variation  in  O 
is  sufficient  to  fulfill  the  condition  |D|  <  |A|  for  the  required  fraction 
of  the  spins.  It  is  more  likely  that  some  pairs  of  spins  are  coupled  by 
an  exchange  five  times  stronger  than  the  average. 

We  investigated  whether  crushing  of  a  crystal  into  a  fine  pow  - 
der  would  strongly  affect  D  or  A  by  comparing  the  linewidth  of  the 
powder  to  the  linewidth  calculated  from  measurements  on  a  single 
crystal  in  various  field  directions.  The  agreement  was  good,  and  no 
change  of  the  constants  by  this  treatment  was  detected. 


CHAPTER  5 


MEASUREMENTS  ON  COPPER  POTASSIUM  SULFATE 
AND  COPPER  SULFATE 

5.1.  Description  of  copper  potassium  sulfate. 

Hydrated  copper  potassium  sulfate,  K2Cu(SO^)2  *  ^^20  ,  is  a 
member  of  the  long  series  of  isomorphic  double  sulfates  known  as 
Tutton  salts.  The  other  salts  in  the  series  can  be  obtained  by  sub¬ 
stituting  for  Cu^”*^  ;  Mg,  Zn,  Cd,  Ni,  Co,  Fe  or  Mn  ;  for  :  Rb,  Cs, 
(NH^)  or  T1  ,  and  the  sulfate  can  be  replaced  by  selenate.  The  salts 
are  monoclinic  with  two  ions  per  unit  cell. 

79 

From  goniometer  and  density  measurements,  Tutton  has 
given  the  topic  axial  ratios  of  K2Cu(SO^)2  ’  8®^  the^di  - 

mansions  of  the  unit  cell  by  multiplying  these  numbers  by  (Z/N)"^  = 
1.489  •  lO”®  cm  where  N  =  6.  025  ♦  10^^  is  Avogadros  number.  We  find 

a  =  9.  040  A  b  =  12.069  A  c  -6. 141  A. 

o  o  o 

The  monoclinic  angle  is 

p  *  104°  28*  . 

80 

Hofmann  has  x-ray  analyzed  a  number  of  the  salts,  particu¬ 
larly  the  ammonia  compounds.  For  (NH^)2Mg(SO^)2  ' 

5 

found  the  space  group  to  be  Czh  so  that  we  can  regard  the  position 
of  the  second  molecule  in  the  unit  cell  as  being  derived  from  the  first 
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,  <2  •  0)  followedby  a  reflection 
in  the  ac -plane.  The  ions  are  surrounded  by  distorted  octahe¬ 

drons  of  water  molecules.  four  of  which  lie  nearly  in  a  square,  each  at 
o  .  ,  o 

a  distance  1.  9  A  from  the  Mg  and  the  other  two  at  2. 15  A  .  The 

crystalline  electric  field  should,  therefore,  have  approximately  tetrag¬ 
onal  symmetry  about  the  line  joining  the  more  distant  waters.  How- 

81 

ever,  this  structure  has  been  criticized  by  Wycoff,  and  anyway,  we 
expect  internal  dimensions  to  be  somewhat  different  in  K2Cu(SO^)2  ’ 

m^o. 

2 

The  field  of  tetragonal  symmetry  will  split  the  D- state  of  the 

free  Cu^^  ion  and  leave  a  Kramers'  doublet  some  20000  cm~^  below 

82 

the  other  orbital  levels.  The  tetragonal  symmetry  will  be  reflected 

via  the  spin-orbit  coupling  in  an  axial  g -tensor  for  each  of  the  two  ions 

in  the  unit  cell.  However,  if  rhombic  distortions  are  also  present, 

there  will  be  g -tensors  with  three  principal  axes,  g  ,  g  and  g  . 

X  y  z 

The  electron  spin  resonance  behavior  of  the  copper  Tutton  salts 

83  84 

has  been  discussed  by  Bleaney  and  co-workers.  ’  Each  magnetic 
complex  can  be  described  by  the  spin- Hamiltonian  (S  *  'j)  i 

Xo  +  TTP^iro"?  (5-1) 

63  65 

The  copper  isotopes  Cu  (natural  abiindance  69  ^  )  and  Cu 

(abundance  31  ^  )  have  both  nuclear  spin  1  »  .  The  tensor  B  gives 

the  electron-nucleus  interaction.  The  tensor  gives  the  electro¬ 
static  interaction  with  the  nuclear  electric  quadrupole  moment.  The 


one  by  a  translation  from  (0,  0,  0) 


to 


no 

last  term  in  (5>1)  is  the  direct  interaction  between  the  nuclear  mag¬ 
netic  moment  and  the  external  field. 

Susceptibility  measurements  on  the  Tutton  salts  are  compli¬ 
cated  by  the  presence  of  two  ions  per  unit  cell  with  different  magnetic 
axes..  If  tetragonal  symmetry  for  each  ion  is  assumed,  g-values  can 
be  deduced  from  the  susceptibilities,  but  this  assumption  may  not  be 
warranted.  In  resonance  measurements,  however,  one  can  find  unam¬ 
biguously  the  magnetic  axis  of  each  complex  and  the  principal  values  of 
the  g-tensors,  but  this  is  tedious  work  if  the  symmetry  is  rhombic. 

In  K2Cu(SO^)2  ‘  ^^2^  the  directions  for  maximum  g  =  g^  is  in 

o  84 

the  crystalline  ab-plane  making  angles  of  about  45  with  the  a-axis. 

In  the  concentrated  salt,  g^  is  reported  to  be  2.  36,  but  we  have  meas- 

ured  g^  =  2.41.  In  the  plane  perpendicular  to  this  direction,  g  varies 

from  2.05  to  2.12  which  signifies  a  considerable  rhombic  component  of 

the  crystalline  field.  In  Zn-diluted  salts,  it  is  found  that  g^  3  2.44  and 

o  83 

that  the  magnetic  axis  makes  angles  of  a  =(42^2)  with  the  a-axis. 

The  spin-spin  interaction  broawlens  the  resonance  in  concentrated 
salts,  and  only  in  magnetically  diluted  salts  are  thehyperfine  lines  seen 
separated.  Normally,  the  allowed  hyperfine  transitions  are  those  in  which 
the  nuclear  orientation  doe  s  not  change  ( 6m  =  0) ,  but  if  the  quadrupole  inter  - 
action  is  comparable  to  the  magnetic  interaction,  transitions  corresponding 
to  Am  -  4  2  are  also  possible.  In  addition.  Am  =  4  1  transitions  are 
allowed  if  the  external  magnetic  field  is  not  along  a  rhombic  axis. 
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The  hyperfine  and  quadrupole  interaction  constants  are  only 

83 

different  for  the  two  isotopes,  and  their  weighted  means  are 
B,  =  (103  +  5)  lO"**  cm'^  , 

Z  " 

,  .4  .1 

-  (34Hh5)M0  cm  (measured  in  ab -plane)  , 

P  =  4p  =  (11  +  1)  10“*  cm~^  , 


P  =  P  is  small 
X  y 


With  the  field  along  the  magnetic  z-axis  of  one  type  of  ion,  the 
hyperfine  lines  corresponding  to  dm  >  0  of  this  ion  are  shifted  in 
energy  by 


m  *  +  3/2 

h  AV  a  +  3/2  B, 
z 

AH  = 

-  138  gauss 

m  =  +  1/2 

h  Av  a  +  1/2  B 

z 

AH  • 

-  46  gauss 

m  a  -  1  /2 

h  Av  a  -  1/2  B 

z 

Ah  a 

46  gauss 

m  a  -  3/2 

h  AV  =  -  3/2  B^ 
z 

AH  = 

138  gauss  .  (5-2) 

Here  we  have  used  g^  ~  2. 41  . 

For  other  orientations,  the  splittings  are  smaller  and  more 
lines  may  appear. 

Some  weak  exchange  has  been  observed  in  K2Cu(SO^)2  *  > 

This  will  be  discussed  in  connection  with  the  results  of  our  first  mo¬ 
ment  shifts  in  Section  5.  3  . 

5.  2.  First  moment  changes. 

We  have  measured  an  approximately  spherical  sample  of 
K2Cu(SO^)2  ’  6H2O  at  30.00  kMc/s,  so  mounted  in  the  waveguide 
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spectrometer  that  the  external  field  could  be  rotated  in  the  crystalline 
ab'plane.  The  field  direction  could  then  be  made  parallel  to  the  tetrag¬ 
onal  axis  of  one  of  the  magnetic  complexes,  and  we  observed  two 
resonances  with  g  =  2.41  (called  line  I  below)  and  with  g  =  2. 11 
(line  II)  .  When  the  field  was  rotated  45°  away  from  this  position,  the 
two  lines  did  overlap  as  expected. 

Line  II  is  relatively  narrow  and  appears  symmetric  at  all  tem¬ 
peratures.  Line  I  is  much  broader  because  of  the  larger  hyperfine 
splitting,  and  it  becomes  very  unsymmetric  at  low  temperatures  as 
shown  in  Fig.  16.  No  separation  of  the  hyperfine  components  is  seen  , 
and  the  line  narrows  at  low  temperature  to  a  half  width  of  about  200 
gauss  as  compared  to  the  predicted  hyperfine  width  of  280  gauss.  In 
Fig.  17a  we  show  the  first  moment  shifts  of  the  two  lines  when  the 
field  is  in  the  direction  of  maximum  line  separation,  and  in  Fig.  17b  we 
show  the  change  in  line  position  when  the  field  is  along  the  crystalline 
b-axis  and  the  two  lines  overlap. 

When  another  sample  of  this  salt,  which  had  considerably 
better  spherical  shape,  was  measured,  we  found  the  resonance  at  the 
line  overlap  position  to  broaden  below  1°K  and  split  up  into  a  number 
of  separate  peaks.  The  positions  of  the  peaks  did  not  correspond  to 
the  positions  of  the  hyperfine  lines  calculated  from  (5-1).  We,  there¬ 
fore,  believe  that  these  peaks  must  be  the  magnetostatic  modes  briefly 
e3q>lained  in  Section  1.  8  ,  and  which  will  be  discussed  again  in  Section 
5.5  . 


I 


Fig.  16  Shape  of  Resonance  Line  I  in  K-Cu(SO.),  ‘  6H.O  as  Function 


Fig.  17  First  Moment  Shift  of  Resonance  in  K2Cu(SO^)2  ’ 

Fig.  17a.  Field  Angle  45^  With  Respect  to  Crystal  Axes  a 
Fig.  17b  Field  Parallel  to  b-Axis. 
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An  exact  calculation  of  the  line  shifts  from  the  first  moment 
theory  is  hardly  possible  because  of  the  rhombic  character  of  the  two 
g-tensors.  We  will  use  formulas  (1-24)  and  (1-43)  with  the  measured 
g -values  and  use  at  first  the  usual  criterion  that  unseparated  lines 
mean  like  spins.  Different  g -values  in  directions  perpendicular  to 
T^q  will  cause  small  changes  in  transition  probabilities  and  possibly  in 
lineshapes  ,  but  any  such  effect  will  be  ignored. 

For  the  field  along  the  b-axis,  all  spins  are  like  and  the  lattice 
sum  is  found  to  be 


I 


(3  cos^  0 


l)/r 


ij 


3.7  •  10^^  cm‘^  . 


(5-3) 


From  (1-24),  the  first  moment  shift  should  then  be  for  g  *  2.  23, 
and  gpH  =  1  cm”^ 


T  = 

0 

0.  5 

1.0 

2.0 

4.0 

8 

o 

<  AH>  = 

57 

51 

36 

19 

10 

0  gauss 

This  calculated  shift  is  plotted  in  Fig.  17b  ,  and  we  see  that  the 
agreement  with  the  experimental  values  is  not  too  bad. 

¥rhen  we  observe  two  lines,  we  have  to  calculate  the  lattice  sums 
separately  for  like  and  unlike  neighbors.  With  IIq  in  the  ab-plane 
making  an  angle  of  45^  with  the  sums  are 


£  (3  co.^  -  l)  /  ry’l 

j 


=  -6.6’  10^^  cm  ^ 


like 


(5-4) 
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[Y  (3  cos^  Q.  -  1)  /r.,  ]  =  5.  2  •  10^^  cm'^ 

^  unlike 

The  dipolar  shift  is  from  (1-43)  , 


(5-5) 


<  AH> 


dip 


1 


like 


4  Tanh 


'  ZTt" 


gj_P  [ 


(5-6) 


For  line  1  we  find  the  shift  by  using  =  2.41,  = 

g Hq  =  1  cm'^  ,  and  for  line  II  we  have  to  use  gy^  =  2*  11  »  g  . 

.1 

and  g«P  =:  1  cm  .  Some  numerical  values  are  ; 


2. 11  and 
=  2.41 


T 

0 

0.5 

1.0 

2.0 

4.0 

o 

8 

<^”>Idip  = 

51 

47 

36 

20 

11 

0  gauss 

^  ”  ^dip  » 

54 

41 

28 

14 

7 

0  gauss  . 

These  calculated  dipolar  shifts  are  indicated  by  the  broken  lines 
in  Fig.  17a  .  The  difference  between  <  6H'>  and  the  experimental 
values  must  be  caused  by  exchange  between  unlike  ions. 

The  first  moment  change,  because  of  exchange  to  z  unlike 
nearest  neighbors,  is  from  (1-43). 


< 


6H>. 


ex 


j  Tanh ( 


2  kT  ’ 


*^ik> 

«i  P 


(5-7) 
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The  dipolar  plus  exchange  shifts  give  the  best  fit  to  the  ex  - 

~  -18 

perimental  values  for  line  II  if  we  choose  A.^i  =-1.4*  10  erg  , 
which  is  a  ferromagnetic  exchange.  Then  we  have  for  z  =  4,  nearest, 
unlike  neighbors 


T 

.  0  0. 5 

1.0 

2.0 

4.0 

°K 

<  Ah 

=lex  = 

-70 

-38 

-20 

0  gauss 

<  AH 

=ll.x  = 

-96 

-56 

-29 

0  gauss 

<*H>«  + 

are 

shown  in 

full  lines  in  Fig.  17a, 

and  it  is  seen  that  the  agreement  with  the  experimental  shifts  of  line  I 
is  very  poor.  Even  if  we  assumed  the  dipolar  sums  to  be  in  error,  it 
does  not  seem  possible  to  explain  the  first  moment  shifts  of  both  lines 
with  the  same  isotropic  exchange  between  unlike  spins.  Neither  will  it 
help  to  postulate  anisotropic  exchange,  since  any  reasonable  type  of 
this  will  shift  the  lines  the  same  way  and  by  about  the  same  amount. 

We  will  now  discuss  qualitatively  a  possible  explanation  for  the 
difference  in  line  shift.  Isotropic  exchange  between  the  spins  reso¬ 
nating  in  line  I  will  not  change  the  first  moment  of  the  line  at  low 
temperature  if  we  do  not  restrict  spin-flips.  This  we  would  have  to 
do  if  the  hyperfine  interaction  were  so  strong  that  we  could  see  the 
line  split  into  21  4  1  =  4  separate  components.  We  do  not  observe  a 
split  line  1  in  K2Cu(SO^)2  ’  •  This  must  mean  that  spin-flips 

are  allowed  at  least  between  neighboring  hyperfine  components  and 
that  a  re-arrangement  within  the  spin  system  can  easily  take  up 
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the  energy  difference  (Am  =  1)  .  But  it  does  not  necessarily  mean 

3 

that  flips  between  spins  of  nuclear  quantum  number  m  s  and  m  = 

3 

-  ^  are  permitted  with  the  same  probability,  since  this  involves  an 

energy  difference  of  (Am  =  3)B  . 

z 

The  idea  here  is  essentially  the  same  as  the  cross  relaxation 

8  5 

discussed  by  Bloembergen  and  co-workers.  Spin-flips  that  do  not 
conserve  energy  are  not  completely  prohibited,  but  they  occur  with  a 
probability  that  decreases  exponentially  as  the  energy  difference  in¬ 
creases.  Here  our  interest  is  in  a  case  where  the  probability  for 
spin-flips  of  Am  =  3  ,  and  perhaps  Am  =  2  ,  is  sufficiently  reduced 
so  that  these  processes  do  not  influence  the  resonance  lineshape.-  But 
considered  as  a  cross  relaxation,  these  flips  may  still  occur  at  a  very 
fast  rate  compared  to  the  spin-lattice  relaxation. 

If  exchange  spin-flips  between  the  hyperfine  components  are 
to  be  only  partly  allowed,  we  must  require  the  order  of  magnitude  ex¬ 
change  field  to  be  somewhat  less  than  280  gauss  from  each  of  the 
two  like  neighbors  6. 14  ^  away.  This  limit  is  considerably  higher 
than  the  dipolar  interaction  between  the  same  spins.  From 

A  S  .  S  . 

we  find  |A|  <  12  •  10"^®  erg. 

If  we,  for  example,  assume  that  spin-flips  of  Am  =  2  and 
Am  B  3  must  be  excluded,  then  on  tiie  average  3/8  of  the  spins  of 


I 
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type  1  would  appear  \inlike  to  other  spins  of  type  I  .  The  resulting 
first  moment  shift  <  A  H  >  is  given  by  (5-7)  when  we  use  z  = 

3 

2  '  To  get  calculated  shifts  in  agreement  with  the  experimental 

~  -18 

data,  we  choose  A  =  >  10  *  10**  °  erg  and  find 
T  *  0  0  .  5  1.0  2.0  4.0  00°  K 

<AH  >  jgjj  ~  -168  -150  -105  -57  -28  0  gauss 

The  sum  ,  <AH  >jjj£p  +  <  +  <  AH  >  is  shown 

as  the  dash-dot  line  in  Fig.  17a  . 

~  -18 

The  choice  of  A  «  -  10  *  10  erg  is  in  good  agreement  with 
the  condition  derived  from  (5-8).  But  it  should  be  stressed  that  a 
different  restriction  on  the  spin-flips  will  give  a  different  result. 

The  hyperfine  broadening  of  line  II  is  much  smaller,  and  all 
spins  of  this  kind  can  probably  be  treated  as  sufficiently  similar  so 
that  the  exchange  between  them  will  not  shift  this  line. 

Because  of  the  hyperfine  contribution  to  linewidth,  no  attempt 
was  made  to  measure  the  second  moment  of  the  lines  or  compare  it 
to  calculated  values. 

5.  3.  Discussion  of  exchange  in  K2Cu(SO^)2  ’  AH2O  . 

Measurements  of  first  moment  shifts  of  resonance  lines  will 
only  give  information  about  exchange  between  non-identical  spins.  This 
may  be  of  value  in  complicated  crystals  to  separate  the  various  ex¬ 
change  interactions.  We  attempted  unsuccessfully  to  explain  the 
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observed  line  shift  in  K2Cu(SO^)2  ‘  basis  of  some  iso¬ 

tropic  exchange  between  the  two  ions  in  the  imit  cell.  However,  we 
could  get  agreement  if  we  assumed  that  the  hyperfine  interaction  for 
spins  in  line  1  was  strong  enough  to  reduce  the  probability  for  the 
spin-flips  within  the  line,  and  in  this  way  we  could  guess  the  order  of 
magnitude  exchange  interaction  between  ions  in  the  same  lattice  posi¬ 
tions.  Before  we  go  on  to  compare  these  results  to  other  data  on  the 
strength  of  the  exchange,  we  should  note  that  this  suggestion  of  partly 
hyperfine -hindered  spin-flips  will  also  help  to  eaqplain  the  observed 
asymmetry  of  line  1.  Some  of  the  spins  are  free  to  spin-flip  with 
their  like  neighbors;  their  resonance  frequency  is  not  shifted,  and  the 
line  should  be  exchange -narrowed.  This  would  account  for  the  rather 
sharp  peak  seen  In  Fig.  16  for  T  =  0.45^  .  Other  spins  are  restricted 
in  their  spin-flips;  their  resonance  frequencies  are  then  shifted  towards 
lower  fields  by  the  ferromagnetic  interaction.  The  maximum  shift 
would  be 

“like  • 

AH  ■  -  =  -450gaues  ,  (5-9) 

and  this  is  about  the  length  of  the  line  tail  we  see  in  Fig.  16  . 

The  Curie  temperature  of  K2Cu(SO^)2  *  ^H20  has  been  derived 

from  susceptibility  and  specific  heat  measurements  on  powdered 

Q-  87 

samples.  OeKlerk  found  0  a  0.054  a;  data  of  Kramers  et  al. 

88 

seemed  to  be  in  agreement  with  this,  but  Garrett  and  Benaie  and 
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Cooke  have  found  0  3  0.034^  and  0  =  0.  035°K,  respectively.  The 

interaction  is  ferromagnetic.  Dipolar  interactions  contribute  to  this 

59 

Curie  temperature,  and  Daniels  has  given  a  rather  complicated 
formula  for  0^.^  in  Tutton  salts,  taking  accotint  of  the  two  different 
ions  in  the  unit  cell.  However,  the  Curie  temperature  from  dipolar 
interactions  will  only  be  of  the  order  0.  OOl^K  ,  and  in  view  of  the  vari¬ 
ations  in  the  measured  0  ,  it  is  of  little  consequence  to  neglect  the  di¬ 
polar  contribution  to  0  . 

We  can  relate  the  Curie  temperature  to  the  strength  of  the  ex¬ 
change  by  the  approximate  formula  (3-241.  Since,  in  this  case,  weprob- 
ably  have  an  exchange  to  =  2  like  nearest  neighbors  and  to 
z^  =  4  unlike  nearest  neighbors,  we  should  write 

(A,  z,  +  A,  z  )  S  (S  +  1) 

g  - - L-i - -  (5_10) 

3  k 

In  the  previous  section  we  found  that  the  first  moment  changes  could 
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best  be  explained  by  assuming  A2  =  '*1>4  *  10  erg.  Then  we  can 
solve  (5-10)  for  A^  and  find 


^2  *2 
*1 


-  7  •  10“^®  erg  for  0  »  0.  035°K 
-12  •  10’^®  erg  for  0  *  0.  054®K 


(5-11) 


This  is  surprisingly  close  to  the  very  speculative  value  A^  * 
-18 

-  10  *  10  erg  derived  on  the  assumption  of  some  restriction  on 
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spin-flips  between  hyperfine  lines.  However,  no  great  faith  should  be 
put  on  this  particular  value  of  ;  all  that  can  safely  be  said  is  that 
it  appears  likely  that  the  exchange  between  like  nearest  neighbors  in 
K2Cu(SO^)2  *  is  somewhat  stronger  than  the  exchange  between 

unlike  neighbors. 

This  also  appears  reasonable  from  the  crystalline  structure. 

80 

From  the  published  atomic  positions  (which  have  been  criticized),  we 

find  the  most  likely  exchange  path  between  like  neighbors  as 

o  o  o 

..  ^.2A  3. 6  A  1.9  A 

- °H,0 - °H,0 -  C**  ■ 


Unlike  neighbors  would  have  to  be  coupled  via  the  longer  path 


Cu 


++ 


o 

1.9  A 


o 

4.  4  A 


o 

1.9  A 


»2® 


»20 


Cu 


++ 


and  one  would  guess  that  this  interaction  is  weaker. 

5.  4.  Description  of  copper  sulfate. 

The  crystal  structure  of  CuSO^  '  has  been  investigated 

89 

by  Beevers  and  Lipson.  .  The  unit  cell  is  triclinic  and  has  the  fol¬ 
lowing  dimensions  : 

o  o  o 

a  =  6.12A  b  =  10.7A  c  =  5. 97  A 

o  o  o 

a  =  82°16'  p  =  107°26'  Y  =  102®40' 

There  are  two  inequivalent  copper  ions  in  the  unit  cell  on  the 
positions  (0,  0,  0)  and  (’j  ’  ?  ’  them  is  surrounded  by  an 
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approximate  octahedron  of  oxygen  atoms,  belonging  to  four  water  mole¬ 
cules  and  two  sulfate  ions.  We  expect  the  crystalline  field  to  be  mainly 
cubic,  but  with  a  tetragonal  component  along  the  line  joining  the  two 
sulfate  oxygens.  The  angle  between  the  tetragonal  axis  of  the  two  cop¬ 
per  ions  is  about  82°  . 

As  in  K2Cu(SO^)2-  6^0,  the  field  of  tetragonal  symmetry  will 
leave  a  Kramers'  doublet  with  an  axial  g-tensor  as  the  lowest  state  of 
the  Cu^^  ion. 

From  susceptibility  measurements,  one  can  deduce  the  values 

g^  3  2.47  and  g^  =  2.  06  ,  and  the  directions  of  the  magnetic  axes 

are  in  good  agreement  with  what  is  expected  from  the  crystalline 
90 

structure.  These  g-values  have  not  been  checked  directly  by  reso¬ 
nance  measurements,  since  the  presence  of  considerable  exchange 
interaction  between  the  unlike  ions  tends  to  broaden  the  lines  and  make 

them  overlap.  Even  with  a  frequency  as  high  as  35  kMc/ s,  Bagguley 
91 

and  Griffiths  are  just  barely  able  to  separate  the  lines  associated 
with  the  two  non -equivalent  crystal  sites.  At  room  temperature,  the 
apparent  g-factors  are  g^  =2.37  and  g^  =  2. 12  .  For  directions  of 
the  field  where  the  two  resonances  overlap,  the  line  is  strongly  ex- 
chang»>narrowed. 

No  hyperfine  lines  have  been  previously  reported  in  CuSO^ 

*  5H2O  since  this  salt  cannot  be  magnetically  diluted  by  Zn,  Mg  or 
Cd  .  However,  we  expect  the  hyperfine  splitting  to  be  of  the  same 
order  as  in  K2Cu(SO^)2  *  6H2O  ^ich  has  similar  g-yalues. 
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From  susceptibility  measurements,  it  has  been  observed  that  the 

exchange  leads  to  an  anisotropic  Curie  temperature  of  about  -0.6^K.^^ 

This  is  supported  by  specific  heat  measurements  by  Geballe  and 
9  3 

Giauque^  who  found  a  specific  heat  maximum  around  1^  .  But  the 
entropy  associated  with  the  integrated  specific  heat  peak  turned  out  to 
be  only  R  ^  2,  half  of  the  expected  value  for  an  ensemble  of  identic 
cal  spins  of  S  =  •  Geballe  and  Giaugue  then  suggested  that  the  ions  in 

the  unit  cell  are  different  and  only  ions  of  one  kind  couple  together  to 
give  a  Curie  temperature  0  =  -0.6°  . 

94 

This  idea  has  been  developed  further  by  Miedema  et  al.  in  a 
recent  paper  that  arrived  only  after  we  had  started  measurements  on 
CxiSO^  ■  5H2O  .  They  measured  susceptibility  and  specific  heat  of  this 
salt  down  to  0.  03°K  and  found  two  Curie  temperatures  and  specific  heat 
anomalies.  The  relatively  strong  antiferromagnetic  coupling  giving 
rise  to  0  3  -0.  6°  they  associate  with  the  ions  on  crystalline  sites 

,  0  )  (eystem  I).  These  spins  are  linked  together  in  a  direction 

parallel  to  the  c-axis  via  two  water  oxygens  as  follows  : 

o  00 


Cu 


1.95  A 


2.6  A 


1.95  A 


'H2O 


'H2O 


Cu . 


A  much  weaker  ferromagnetic  coupling  giving  another  Curie 
temperature  of  about  0. 02°K  they  associate  with  the  spins  on  sites 
(0,  0,  0)  (iystem  11),  which  are  coupled  via  a  minimum  of  three  oxygens 
between  equal  Cu^^  neighbors. 
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The  two  systems  are  coupled  together  with  an  exchange  strength 
that  is  very  roughly  estimated  from  the  resonance  lineshape  to  be 
about  0.15  cm  ^  at  room  temperature.  We  find  this  exchange  to 

increase  with  decreasing  temperature,  and  it  is  perhaps  50^  stronger 
at  helium  temperatures. 

The  coupling  between  the  spin  systems  has  a  consequence  for 

94 

susceptibility  measurements  not  pointed  out  by  Miedema  et  al.  At  a 
temperature  so  low  that  system  I  already  is  in  the  antiferromagnetic 
state,  this  coupling  between  unlike  spins  will  tend  to  align  system  11 
antiferromagnetically.  The  observed  weak  ferromagnetic  coupling  in 
system  11  must  then  be  the  net  effect  of  a  ferromagnetic  coupling  be¬ 
tween  like  spins  almost  cancelled  by  a  coupling  to  the  antiferromag¬ 
netically  ordered  system  of  tmlike  spins. 

5.  5.  Measurements  on  copper  sulfate  and  discussion  of  observed 
magnetostatic  modes. 

We  have  measured  an  almost  perfect  sphere  of  CuSO^-  5H2O 
in  the  waveguide  spectrometer.  However,  it  was  very  hard  to  x-ray 
orient  the  sample  since  we  had  no  clue  from  external  faces  to  explain 
the  lAue  patterns  from  this  triclinic  crystal.  In  the  end,  we,  there¬ 
fore,  looked  at  the  resonance  of  the  sample  in  the  cavity  spectrometer 
at  room  temperature  ^ere  it  was  rather  easy  to  move  it  around  until 
we  saw  the  two  resonances  fairly  well  separated.  Then  the  field  had 
to  be  almost  parallel  to  one  of  the  tetragonal  axes,  and  this  was  the 
orientation  we  used  for  the  low  temperature  work. 
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But  at  77*^K  the  two  resonances  had  merged  into  one  broad  line. 
This  might  be  caused  by  : 

a)  The  principal  values  of  the  g-tensors  change  in  this  temperature 
range . 

b)  The  exchange  betweenunlike  spins  changes. 

With  kind  help  from  Dr.  Foner  at  M.  I.T.  we  attempted  to  find 
out  what  was  happening  by  measuring  the  crystal  at  105  kMc/s  at 
300°K  and  77°K  .  At  this  frequency  the  resonances  were  far  apart  at 
room  temperature,  but  they  broadened  and  moved  somewhat  towards 
each  other  as  the  temperature  was  lowered.  Although  die  evidence  is 
not  quite  conclusive,  the  exchange  appeared  to  increase  when  the 
sample  was  cooled.  Qualitatively,  this  is  what  we  expect  to  happen 
when  the  lattice  contracts  and  the  atoms  move  closer  together. 

If  we  assume  that  the  effect  of  the  exchange  in  this  case  can  be 
approximated  by  the  Markoffian  random  modulation  model  developed  by 

•7 

Anderson,  we  see  from  his  Fig.  2  that  w  /w  must  change  from  about 

0.  5  at  room  temperature  to  at  least  0.7  at  77°K  .  This  increase  of 

40  ^  in  the  exchange  is  rather  more  than  one  might  expect  from  a 

slight  thermal  contraction  of  the  lattice. 

A  similar  and  even  larger  increase  in  the  exchange  is  observed 
97 

by  Abe  in  K2CuCl^  *  and  we  have  confirmed  this  result.  At 

room  temperature  one  has  to  measure  at  a  frequency  larger  than  50 

kMc/s  to  see  two  separate  resonances  from  the  two  differently  oriented 

98 

magnetic  complexes  in  the  unit  cell  of  this  salt.  For  a  frequency  of 
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30  k^c/s,  the  composite  line  is  about  1200  gauss  wide  at  300^  ,  but 
at  77°K  we  find  the  width  narrowed  to  about  200  gauss.  This  is  the 
order  of  magnitude  exchange -narrowing  we  expect  from  the  interaction 
giving  a  reported  Curie  temperature  of  0.  8^  K  . 

Fig,  18  shows  the  change  of  the  resonance  in  CuSO^*  5H2O 
at  low  temperatures  when  the  field  is  in  a  direction  where  the  two 
resonances  are  separated  at  room  temperature.  The  first  moment  of 
the  line  clearly  changes  as  we  reduce  the  temperature,  but  we  will  not 
attempt  to  explain  this  quantitatively  from  spin- spin  interaction,  since 
the  exchange  couplings  are  complex  and  the  g -values  and  sample 
orientation  are  uncertain.  However,  the  line  behavior  below  1^  is 
striking .  The  line  broadens  and  splits  up  into  a  number  of  separate 
peaks.  A  similar  behavior  of  the  resonance  is  seen  for  other  orienta¬ 
tions  of  the  magnetic  field.  We  have  also  observed  peaks  of  the  same 
kind  in  a  good  sphere  of  K2Cu(SO^)2  *  ^^2^  0.  5^  K  when  Hq 

was  in  a  direction  where  the  two  resonances  overlap. 

We  first  thought  that  these  were  hyperfine  lines  although  their 
positions  did  not  agree  very  well  with  the  spin -Hamiltonian  (5-1)  for 
K2Cu(SO^)2  ’  ^^2^  ‘  random  modulation  model  of  ex¬ 

change -narrowing  by  Anderson  and  Weiss, ^  one  may  be  misled  to  be¬ 
lieve  that  the  spin-flips  should  be  frozen  out  and  the  modulation  cease 
at  T  s  0  .  This  is  not  so.  If  we  are  to  observe  the  resonance,  some 
spins  must  be  excited  by  the  microwave  field  and  these  spins  are  then 
free  to  flip  with  their  neighbors  in  the  gro\md  state  provided  that  the 


126 


exchange  between  them  is  stronger  than  the  possible  difference  in  hy- 
perfine  or  Zeeman  energy.  Hyperfine  lines  masked  by  spin-spin 
interaction  at  high  temperature  should,  therefore,  not  appear  separated 
when  T  — ^  0  ,  and  this  cannot  explain  the  observed  lineshape. 

It  is  more  likely  that  these  peaks  are  the  magnetostatic  modes 
that,  in  principle,  may  appear  in  the  resonance  of  any  magnetized  body 
regardless  of  whether  the  spin  alignment  is  caused  by  strong  exchange 
or  by  low  temperature  and  strong  field.  In  K^u.Cl^  •  2H2O  ,  such 

19 

modes  have  been  observed  by  Abe  et  al.  to  branch  out  from  the  main 
resonance  below  2.  2°  K  when  measured  at  25kMc/s.  The  micro - 
wave  field  in  our  experiment  is  probably  sufficiently  nonuniform  to 
make  excitation  of  magnetostatic  modes  possible.  If  so,  why  should 
they  not  then  appear  around  Z°K  as  in  K^CuCl^  *  2H2O?  The  satura¬ 
tion  magnetizatipn  of  CuSO^*  ^^2^  about  20  ^  larger  than  in 
K2CU  Cl^  *  ^^2^  ’  also  did  the  work  in  a  stronger  field  than 

Abe  et  al. 

This  can  be  explained  by  the  difference  between  the  ferromag¬ 
netic  coupling  in  K2CU  CI^  ’  ^^2^  external  field  in 

aligning  the  spins  and  the  antiferromagnetic  coupling  in  spin  system  I 
in  Cu  SO^  ’  5H2O  that  opposes  alignment.  The  nnagnetization  p(T) 
relative  to  the  saturation  magnetization  of  K2CU  Cl^  *  2H2O  in  a 
magnetic  field  Hq  can  be  found  by  solving  graphically  the  equation 

p  =  Tanh  (  (<1  g  p  Hq  pA)/kT  )  .  (5-13) 
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The  strength  of  the  exchange  A  is  given  by  the  condition  that 
for  Hq  — ^  0  the  equation  (5-13)  has  to  give  spontaneous  magnetiza¬ 
tion  for  the  Curie  temperature  T  =  0.  8^  K  ,  so  we  find  A/k  =  0.8  . 
Then  we  get  for  g  ^  Hq  =  25  kMc/s 


T 

= 

2.2 

2.0 

1.8 

1.  5 

1. 

0 

0 

p  =  M(T)/M^ 

= 

0.40 

0.45 

0.  51 

0.62 

0. 

86 

P(A  =  0) 

= 

0.265 

0.29 

0.32 

0.  38 

0. 

54  . 

Because  of  the  strong  exchange,  the  magnetization  is  consider¬ 
ably  stronger  at  a  given  temperature  than  xdiat  it  would  be  in  the  ex¬ 
ternal  field  alone,  which  is  given  as  p^^  ^  q^. 

In  CuSO^  *  5H2O  ,  spin  system  I  can  give  only  a  small  contri¬ 
bution  to  the  magnetization  since  the  antiferromagnetic  coupling  in  this 
system  giving  a  Curie  temperature  of  -0.  60°  is  comparable  to  die 
force  of  the  external  field.  Then  only  below  1°  K  is  system  U  80-90 
aligned  for  gpH  =  30  kMc/s  ,  and  only  there  can  this  system  alonjS 
give  a  magnetization  comparable  to  that  of  K2CU  Cl^*  6H2O  around 
2°  K  . 

The  maximum  separation  of  modes  in  CuSO^  •  5H2O  is  about 

200  gauss.  This  is  a  reasonable  splitting  since  hrom  Walker's 
20 

theory,  the  maximum  separation  should  be  6H  =  2x  M  .  At  T  m  0.  5  , 
the  magnetization  in  CuSO^  *  5H2O  is  M  w-j  *  30  gauss  for 
gP  Hq  =  1  cm''^  so  that  2v  M  m  200  gauss. 
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In  K2Cu(SO^)2  ‘  ,  the  saturation  noagnetization  3 

S  gP  N|  is  only  about  32  gauss  compared  to  about  47  gauss  in 
K^Cu  Cl^  '  2H2O  .  Since  the  exchange  in  this  salt  is  weak,  it  also  has 
to  be  cooled  below  l^K  before  the  magnetization  becomes  comparable 
to  that  of  K^Cu  Cl^  •  2H^O  around  2*^K  .  We  tried  to  measure  the 
position  of  the  modes  in  K2Cu(SO^)2*  ^  function  of  the  tem¬ 

perature,  but  we  could  not  clearly  establish  that  the  splitting  was  pro¬ 
portional  to  the  magnetization  as  required  for  Walker  modes  because 
the  modes  were  well  separated  only  below  0.  5^  K. 

The  conclusion  that  our  observed  peaks  are  magnetostatic 
modes  is  also  supported  by  the  fact  that  we  see  them  only  in  the  most 
spherical  samples.  In  samples  of  more  irregular  shape,  the  modes 
will  be  poorly  defined,  broad  and  difficult  to  excite.  The  use  of  poor 
spheres  may  be  the  reason  why  we  did  not  see  them  in  Nd(C2 
SO^)^  *  9H2O  and  NiSiF^  *  6H2O  where  the  saturation  magnetization 
is  28  and  86  gauss,  respectively.  We  then  have  the  paradox  that  a 
sample  of  less  than  perfect  shape  may  actually  give  the  narrowest 
resonance  line  when  measured  below  1^  K. 

Whether  magnetostatic  modes  can  be  seen  as  separate  lines 
obviously  also  depends  upon  how  wide  the  imiform  resonance  is. 
Excited  modes  of  nearly  the  same  energy  will  give  only  one  broad  line. 
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CHAPTER  6 


RESONANCE  LINE  BROADENING  FROM  VARIATION 
IN  CRYSTALLINE  FIELD 

6.1.  Introduction. 

The  resonance  linewidth  experiments  described  in  this  chapter 

were  undertaken  to  investigate  the  practical  possibility  of  a  frequency- 

stable,  zero-field  maser  oscillator  suggested  by  Professor  Bloember- 
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gen.  His  idea  is  as  follows: 

A  paramagnetic  material  for  a  maser  oscillator  used  as  a 
frequency  standard  should  have  very  narrow  resonance  lines.  Inter¬ 
action  between  like  spins  can  be  reduced  to  give  negligible  line 
broadening  by  diluting  the  crystal  with  non-magnetic  material.  But 
the  nuclear  moments  of  protons  and  other  magnetic  nuclei  will  still 
give  a  local  field  variation  of  a  few  gauss.  This  will  normally  result 
in  a  minimum  linewidth  of  the  order  of  20  Me  /  s ,  and  when  such  a  salt 
is  used  in  a  maser  oscillator,  the  resulting  frequency  will  not  be 
particularly  stable. 

But  consider  an  ion  with  integral  spin  in  a  crystalline  field  of 
rhombic  or  lower  symmetry.  The  2S  •h  1  spin  levels  are  then  com¬ 
pletely  split  by  the  crystalline  field,  and  this  makes  maser  operation 
possible  without  an  external  field.  If  a  weak  field  Hq  is  applied, 
the  magnetic  moment  operator  has  no  diagonal  matrix  elements  in  the 
representation  where  the  spin-Hamiltonian  is  diagonalized.  The 
energy  levels  will  then  only  depend  on  Hq  in  the  second  order,  and 
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the  local  field  fluctuations  from  unlike  spin  neighbors  can  only  cause 
a  negligible  broadening  of  the  resonance.  Interaction  between  like 
ions,  however,  will  broaden  the  line  because  of  the  spin-flips,  but  this 
effect  can  be  reduced  by  diluting  the  crystal. 

An  example  of  a  possible  zero-field  maser  system  is  the  Ni^^ 
ion  of  spin  S  =  1  in  the  rhombic  field  of  hydrated  double  sulfates.  The 
structure  of  the  monoclinic  Tutton  salts  were  described  in  Section  5. 1. 
There  are  two  ions  in  the  unit  cell  and  their  magnetic  axes  are  related 
by  reflection  in  the  symmetry  plane  ac  of  the  crystal.  This  is  shown 
schematically  in  Fig.  19a,:  where  the  magnetic  axes  for  each  complex 
are  labelled  x,  y,  and  z  ;  the  macroscopic  susceptibility  axes  are 
,K2  ,  andK^  ,  and  the  monoclinic  crystal  axes  are  a,  b,  and  c. 

The  magnetic  properties  of  each  complex  can  be  described  by 
a  spin- Hamiltonian^ 

Kq  =  gP  i?o  S  +  D(S^^  -  I )  +  E(S^^  -  Sy2)  .  (6-1) 

When  the  external  field  is  along  the  z-axis  of  one  spin,  its 
energy  levels  are 

1 

\  J  =  yD  +(E^  +g^P^  ^0^^^  y  D  +  E  +  g^p^  ^0^^ 


X3  =  jD-(E2+gVHQ^)^j,-::^0  jo  -  E  -  gVH^^/ZE  .  (6-2) 
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The  rhombic  field  mixes  the  spin  states  so  that  all  transitions 

X. .  __  X  .  are  allowed. 

1  J 

Griffiths  and  Owen  ^  have  investigated  a  number  of  the  Ni  - 
Tutton  salts.  The  crystalline  field  parameters  D  and  E  and  the 
direction  of  the  magnetic  axes  varies  somewhat  from  salt  to  salt,  and 
these  parameters  are  a!so  fimctions  of  temperature.  For  Ni(NH^)2  ' 
(SO^)^  •  .  they  find  at  90®  K 

g  =  2.25  D  =  -l.  99cm"^  E  =  -0.48cm"^  . 

The  energy  levels  (6-2)  with  these  parameters  are  shown  in 
Fig.  19b  as  functions  of  the  magnetic  field  in  the  z-direction.  The  z- 
axes  of  the  two  complexes  almost  coincide  in  this  salt;  the  angle 
Zj  0  z^  shown  in  Fig.  19a  is  7®  . 

A  zero-field  maser  using  Ni(NH^)2(SO^)2  *  would 

oscillate  at  about  29  kMc/s  (X^ —  X^)  and  it  would  require  a  pump 
frequency  of  about  75  kMc/s  (X^-_  X2)'  The  theoretical  linewidth  of 
this  salt  when  highly  diluted  should  be  comparable  to  the  4  kc/s  width 
of  the  NH^  -transition  in  the  beam  maser.  However,  this  would 

7 

require  the  crystalline  field  to  be  constant  to  1  part  in  10  .  Strains 
and  crystal  distortions  around  defects  and  impurities  will  probably 
make  such  narrow  lines  impossible.  But  even  for  a  width  of  1  Me/ s, 
Bloembergen  estimated  that  the  frequency  stability  of  such  a  maser 
oscillator  might  compare  favorably  to  the  stability  of  about  1  :  10^^ 
of  the  NH^  beam  maser.  This  is  because  the  solid  state  maser 
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oscillates  at  a  higher  power  level,  gives  off  negligible  noise  and  is 
little  influenced  by  temperature  and  pressure  variations. 

Before  it  could  be  decided  whether  a  zero-field  maser  would 
be  worth  building,  we  needed  more  information  about  linewidths  in 
diluted  crystals.  In  the  next  section  some  measurements  on  nickel 
Tutton  salts  will  be  described. 

6.2.  Measurements . 

For  the  design  of  a  zero-field  maser,  the  resonance  linewidth 
in  exactly  zero  field  is  of  interest.  This  can,  in  principle,  be 
measured  with  a  frequency- swept  microwave  source.  However,  the 
microwave  path  in  a  low  temperature  spectrometer  has  to  be  very  long,  a 
fewhundred  wavelengths  at  a  frequency  of  30  kMc/s.  The  small  reflec¬ 
tions  at  every  waveguide  joint  and  bend  makes  the  response  of  the  system 
vary  widely  with  frequency  if  not  continuously  retuned  .  In  practice ,  only 
lines  narrower  than  a  few  Me /s  can  be  measured  directly  this  way. 

We  attempted  to  get  around  the  problem  of  reflections  in  the 
waveguide  by  changing  the  frequency  step-wise,  tuning  up  the 
transmission  line  in  which  the  sample  was  sitting  for  each  frequency 
and  measuring  the  difference  in  transmitted  power  between  periods 
of  no  field  at  the  sample  smd  of  a  field  of  a  few  thoussmd  gauss 
there.  Since  the  measurements  were  made  at  liquid  helium 
temperature,  we  could  make  the  field  coil  of  superconducting  niobium 
wire.  A  coil  of  7650  turns  of  0.005  "  wire  wound  directly 
around  the  stainless  steel  waveguide  was  calculated  to  give 
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a  field  of  about  2000  gauss  for  a  current  of  1  A,  This  was  the  maxi¬ 
mum  field  the  niobium  could  stand  in  the  superconducting  state.  If 
this  field  is  applied  along  the  magnetic  z-axis,  we  find  from  (6-2)  that 
it  would  shift  the  position  of  the  resonance  (\.  ^  —  \  ^)  by  2600  Mc/s  . 
We  used  a  circuit  built  by  Dr.  Jeener  to  switch  the  field  on  and  off 
once  a  second  and  to  integrate  the  difference  between  transmitted 
microwave  power  during  the  field  "on"  and  "off  periods. 

However,  this  spectrometer  design  was  not  successful.  The 
resonant  behavior  of  the  long  waveguide  and  the  dispersion  effect  from 
the  sample  resulted  in  a  great  number  of  apparent  absorptions.  We, 
therefore,  had  to  give  up  the  attempt  to  measure  line  widths  in  exactly 
zero  field  and  began  measurements  at  constant  frequencies  and  small, 
swept  fields.  Then  the  lineshape  recorded  in  field  units  had  to  be 
graphically  converted  to  the  more  meaningful  frequency  units,  since  the 
resonance  shift  was  not  a  linear  function  of  the  field.  The  linewidths 
in  zero  field  could  only  be  found  by  extrapolating  the  results  in  finite 
fields. 

We  used  a  waveguide  reflection  spectrometer  with  simple  crys¬ 
tal  detection  similar  in  principle  to  the  one  described  in  Section  2. 1  . 
The  fixed  cavity  version  could  not  be  used  since  we  had  to  measure  at 
a  number  of  frequencies  close  to  the  zero-field  resonance.  For  the 
same  reason,  we  chose  aVarian  X-12  klystron  followed  by  a  crystal 
harmonic  generator  as  our  microwave  source.  This  combination  would 
give  a  useful  output  over  the  frequency  ranges  25  -  33  kMc/s  and 
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38  -  43  kMc/s,  which  is  a  much  wider  band  than  could  have  been 

covered  directly  by  any  single  tube.  The  magnet  field  available  was 

limited  by  power  supply  arcing  to  about  6000  gauss.  This  field  limit 

restricted  the  number  of  resonances  that  could  be  observed. 

We  measured  at  4.  2^K  some  nickel  ammonia  double  sulfates 

with  various  percentages  of  dilution  by  Zn,  Mg,  and  Cd  .  The  crystals 

were  generally  cut  rectangularly  with  tapered  ends  .and  pressed  in  place 

against  the  waveguide  wall  by  a  polyfoam  plug .  With  the  crystal 
— k 

b-axis  vertical,  the  magnet  was  rotated  until  the  field  was  along  the 

susceptibility  axis  .  This  direction  is  close  to  the  z-axes  of  both 

ions  and  the  two  resonances  overlap.  Around  30  kMc/s,  we  could 

2  2  2  2  ^ 

observe  the  resonance  (Xj  —  X^)  =  2(E  +  g  P  Hq  )  .  Working 

with  the  third  harmonic  from  the  multiplier  around  40  kMc/s,  we  could 
usually  also  see  the  transition  (X  ^  —  X.  j)  *  D  +  (E^  +  g^P^  ^0^^^  * 

In  the  concentrated  Ni(NH^)2  ‘  (SO^)2  *  6H2O  ,  the  linewidth 
was  very  large,  2300  4  500  Mc/s  .  This  is  more  than  we  e3q>ect 
from  dipolar  interaction  alone  and  signifies  some  exchange  coupling. 

In  the  diluted  salts,  the  linewidth  was  less,  but  dilution  below  about 
1  ^  Ni^^  concentration  did  not  seem  to  reduce  the  width  any  more.  The 
residual  linewidth  found  in  salts  of  nominal  Ni^^  concentration  of 
0.1  ^  together  with  measured  values  of  the  crystalline  field  parameters 
is  listed  in  Table  (6-1).  Approximate  ionic  radii^^^  of  Ni^^  and  di¬ 
lutants  are  included  for  purpose  of  comparison  . 
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Dilutant 

Ionic 

Radius 

|DI 

|E| 

>•> 

8 

Mg'^'*‘ 

0.75A 

1.895+0.003  cm"^ 

0.466+0.001  cm  ^ 

2.2440.01 

70+20  Mc/s 

Zn”'”'* 

0.83 

1.98  +0.01 

0.48240.001 

2.23+0.01 

130+30 

Cd'*’**’ 

0.99 

>  2or  <  1 

0.38540.001 

2.16+0.01 

230+50 

0.74 

0.48 

2.25 

2300+500 

TABLE  (6-1) 

The  apparent  g  listed  here  is  considerably  smaller  in  the 
Cd'diluted  salt  than  the  usual  value  2.  25  for  Ni^^  .  This  probably 

reflects  a  change  in  the  orientation  of  the  magnetic  axis,  and  the 

||{ 

small  g  can  be  explained  if  the  angle  K.^  Oz^  in  Fig.  19a  has  in¬ 
creased  from  the  value  3.5°  in  the  undiluted  salt.  No  attempt  was 
made  to  measure  the  orientation  of  magnetic  axes  in  the  diluted  salts. 

The  residual  linewidth  was  not  field  dependent  down  to  the  min¬ 
imum  field  of  about  500  gauss  that  we  could  work  with.  The  widths 
must,  therefore,  be  caused  by  variations  of  the  crystalline  field  param¬ 
eters  D  and  E  ,  rather  than  by  a  variation  in  g  or  the  orientation  of 
the  magnetic  axes.  The  transitions  and  (X2WX2)  de¬ 

pend  on  D  and  E  in  different  ways.  No  difference  in  the  widths  of 
these  lines  was  found,  and  the  absolute  variations  in  O  and  E  must, 
therefore,  be  approximately  the  same.  The  residual  linewidth 
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remained  essentially  the  same  when  the  magnet  was  rotated  so  that  the 
resonance  frequency  was  a  very  strong  function  of  field  angle.  From 
this,  we  find  that  the  possible  variation  in  orientation  of  the  magnetic 
axes  is  less  than  1/4*^  . 

6.3.  Discussion. 

The  resonance  line  moments  in  undiluted  Ni  Tutton  salt  have 
been  calculated  by  Ollom^^  and  Stevens^^^  for  the  high  temperature 
limit.  For  Ni(NH^)2  (SeO^)  *  6H2O  in  zero  field,  the  limiting  value 
of  Stevens'  expression  for  the  second  moment  of  the  transition 
—  \  2)  caused  by  dipolar  interaction  is 

A  E  ^  >  - >  [0.  252  11  +  0.  270  +  0.  044  •  4]  lO""^  cm"^ 

H  ”^0 

=  0.29  (kMc/s)^  .  (6-3) 

We  expect  the  second  moment  for  Ni(NH^)2(SO^)2  •  6H2O  to 

be  approximately  the  same.  Our  experimental  value  in  near  zero 

field  for  the  transition  (  — •  X  is  <AE^=  1.4  (kMc/s)^  ,  however, 

and  this  means  that  most  of  the  observed  linewidth  must  be  caused  by 

exchange.  Stevens  assumes  a  ferromagnetic  exchange  of  A  =  -  5  * 

-18 

10  erg  to  each  of  the  two  nearest  like  neighbors  and  finds  then 

agreement  with  the  second  moment  measurements  of  Griffiths  and 

^  101 
Owen. 

The  line  moments  as  function  of  temperature  can  be  found  in 
the  way  outlined  in  Chapter  1,  but  the  calculation  is  complicated  by  the 
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two  different  ions  in  the  unit  cell  and  the  rhombic  crystal  field  that 
mixes  the  spin  states.  Accurate  comparison  between  calculated  mo¬ 
ments  and  experimental  data  obtained  in  low  fields  is  difficult,  because 
of  the  nonlinear  relation  between  field  and  resonance  frequency  and 
because  of  the  field  dependence  of  the  transition  probability. 

Kittel  and  Abrahams  have  considered  the  line  shape  in 
randomly  diluted  crystals.  They  find  the  second  moment  from  di¬ 
polar  interaction  to  be  proportional  to  the  spin  concentration  f 

<  Av^>,  ,  =  f<Av^>  (6-4) 

dll  cone 

The  fourth  moment  decreases  relatively  less  with  dilution  so 
that  the  lineshape  will  change  from  approximately  Gaussian  for  f  =  1 
to  the  Lorentzian  with  cut-off  wings  for  f  <  0.  01  . 

The  exchange  contribution  to  the  second  moment  in  diluted 
salts  will  be  in  the  form  of  satellite  lines  for  the  small  fraction  of 
spins  that  happen  to  have  another  spin  as  nearest  neighbor.  The  reso¬ 
nance  of  all  the  other  spins  will  not  be  influenced  by  the  exchange. 

For  Ni-Tutton  salt  diluted  to  f  •  0.  001  ,  the  theoretical 
second  moment  in  zero  field  is  from  (6-3)  and  (6-4) 

1 

[  <  A  17  Mc/s  . 

Most  of  the  contribution  to  the  second  moment  ctxnes  from  the 
wings  of  the  line,  and  the  theoretical  half -width  is  considerably  less  than 
17  Mc/s.  We  observed  a  minimxun  of  about  100  Mc/s  for  the  linewidth 
and  the  second  moment.  The  resonance  lines  in  diluted  salts  must,  there¬ 
fore,  be  broadened  by  a  variation  in  the  crystalline  field. 
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We  see  in  Table  (6~1)  that  the' minimum  linewidths  observed 

80 

vary  with  the  dilutant.  It  is  noted  that  increased  sise  of  the  divalent 

metallic  ion  will  reduce  the  stability  of  the  Tutton  salt.  We  found  that 

the  cadmium  salt  disintegrated  after  a  few  weeks  at  room  temperature. 

The  trend  to  reduced  stability  is  reflected  in  greater  variation  in 

crystal  field  as  seen  from  the  table. 

Other  experiments  might  have  been  undertaken  to  find  the  effect 

on  linewidth  if  or  Tl^  is  substituted  for  ,  (SeO^)  for 

(SO^)""  or  D2O  for  H^O  .  But  it  is  very  unlikely  that  lines  any> 

where  narrow  enough  for  the  ultrastable  maser  oscillator  can  be  found 

in  hydrated  nickel  salts.  This  view  is  supported  by  the  minimum  line- 

72 

width  of  30  gauss  m  95  Mc/s  found  by  Hoskins,  Pastor  and  Trigger 
in  1000  :  1  Zn-diluted  NiSiF^  *  6H2O  .  Similar  linewidths  are  found 
in  diluted  nickel  double  nitrates. 

In  the  crudest  calculation  of  the  spin-Hamiltonian  parameters  from 

first  principles ,  the  water  dipoles  surrounding  the  nickel  ion  is  replaced 
107 

by  point  charges.  This  will  give  an  electric  field  of  known  symmetry  at 

the  position  of  the  magnetic  ion,  and  the  stark  splitting  of  its  orbital  levels 
can  be  calculated  as  a  function  of  point  charge  strength.  Since  the  polari¬ 
sation  of  the  water  molecules  in  the  crystal  is  unknown,  the  static  field 
must  be  found  by  fitting  the  orbital  levels  to  spectroscopic  data.  The  spin- 
orbit  coupling  is  then  applied  as  a  perturbation.  In  tiiis  way ,  one  can,  in 
principle,  relate  the  constants  D  and  £  and  their  variations  to  certain 
rhombic  distortions  of  the  point  charge  model .  This  will  not  be  attempted 
here. 
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Qualitatively,  the  variation  in  D  and  E  and,  therefore,  in  the 
low-symmetry  part  of  the  crystalline  field,  is  about  1  part  in  300.  In 
the  point-charge  model  these  field  components  will  vary  with  at  least 
the  third  power  of  the  interatomic  distance.  Thus,  a  variation  of  this 
distance  of  only  1  part  in  1000  will  result  in  the  observed  spread  in  the 
values  of  D  and  E  .  A  similar  displacement  of  the  point  charge  in 
the  perpendicular  direction  will  only  give  a  variation  in  the  orientation 
of  the  magnetic  axes  of  about  0.  001  radian  ^  (  l/20)^.  Such  an  angular 
variation  would  have  a  negligible  effect  on  the  width  of  the  resonance 
line,  and  so  the  observed  broadening  is  caused  by  a  variation  in  the 
crystalline  field  parameters  D  and  E  . 


CONCLUSION 
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In  the  first  five  chapters  of  this  report,  a  nvimber  of  para¬ 
magnetic  resonance  lineshapes  observed  under  conditions  where 
h  V  >  kT  have  been  analyzed  and  explained  in  terms  of  spin? spin 
interaction.  It  is  shown  that  this  technique  may  be  an  important  sup¬ 
plement  to  other  methods  of  investigating  interaction  in  paramagnetic 
salts.  The  advantages  of  the  line  moment  method  at  low  temperature 
are: 

a)  The  exchange  interaction  can  be  found  in  sign  and  mag¬ 
nitude.  Possible  anisotropy  of  the  exchange  can  be  detected. 

b)  Exchange  interaction  between  unlike  neighbors  only  can  be 
measured.  If  the  total  exchange  is  known  from  other  data,  the  ex¬ 
change  between  like  neighbors  can  be  separated  out. 

c)  The  method  works  best  at  rather  weak  interactions  where 
specific  heat  and  susceptibility  experiments  are  difficult  to  perform 
with  sufficient  accuracy. 

It  should  be  stressed,  however,  that  the  lineshapes  can  only  be 
completely  analyzed  in  relatively  simple  systems,  preferably  the  ones 
with  only  one  magnetic  ion  per  xmit  cell. 

The  disadvantages  of  the  method  are: 

a)  A  complex  He  cryostat  is  needed  and  it  is  difficult 
to  mount  the  sample  inside  the  cavity  in  good  heat  contact  with  the 

3 

He  reservoir. 
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b)  Spherical  samples  are  needed  to  avoid  demagnetizing  effects. 
However,  magnetostatic  modes  are  easily  excited  in  spheres,  and  this 
may  obscure  the  true  line  shape. 

The  experiments  described  in  Chapter  6  were  unsuccessful  in 
the  sense  that  we  did  not  find  lines  narrow  enough  for  the  ultra-stable 
zero-field  maser  proposed  by  Professor  Bloembergen.  It  is  difficult 
to  relate  the  observed  small  variation  in  crystalline  field  to  definite 
distortions  of  the  crystal  lattice. 
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APPENDIX  I 


THE  MOMENT  FORMULAS  OF  McMILLAN  AND  OPECHOWSKI  FOR 
INTERACTION  BETWEEN  LIKE  NEIGHBORS®’^® 

3 

Following  Pryce  and  Stevens,  McMillan  and  Opechowski  intro¬ 
duce  the  projection  operators  P^  and  Pp  that  describe  the  mani - 
folds  of  unperturbed  energy  states  E^  and  E^  for  the  spin  system. 
Such  operators  are  defined  by 


where  |p,k)  satisfies  |r»^)  =  p,k).  We  can  then  write 

5? 


(A-1) 


(A- 2) 


and  when  the  perturbation  is  small, 

-  \  t 

P  a 


A  ^ 

M^  »  ^  P-  M  P 
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«»P 
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(A-3) 


(A-4) 
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a.p 


where  ^  means  summation  over  all  values  of  a  and  p  such  that 

®»P  e  e 

Ep  -  E^  =  hv  where  v  is  the  unperturbed  frequency. 

E3q>ansion  of  the  density  matrix  gives 


P  XP 


=  b 


E  P 
a  a 


[1-  -a 


P  P 


kT 


+  . . .  ]  . 


(A- 5) 


A-2 


In  the  first  approximation,  only  the  first  term  in  the  bracket  is 
kept,  and  this  approximation  is  identical  to  our  Equation  (1-14).  The 
commutator  expressions  (1-11)  and  (1-12)  for  the  moments  around  v 
then  become 


E 

h  <Av  X  ^  ^ 

a,p 

-  P  P  MP-  M)  (A-6) 

a  Q  p  '  '  ' 

E 

h^<Av^>j  =  ^  b  ®  Trace  (PpJC^^^  PpX^^^PpMP^.M 

a,p 

-  2P  P  MP-  P-M  +  P  X^^^  P  X^^^  P  M  P-  M),  (A-7) 
a  aP  pa  a  ap 

where 
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Trace  P  MP^M  . 
a  p 


(A-8) 


As  in  Section  1. 2  ,  these  operators  have  to  be  evaluated  with 
respect  to  the  eigenfunctions  ja)^  ]b)2  •  • .  •  where  |c)^ 

satisfies 


X.^°^  |c).  =ajc)i 

3 

How  this  can  be  done  has  been  indicated  by  Pryce  and  Stevens,  but  the 

18 

detailed  calculation  takes  about  20  pages  in  McMillan's  M. Sc. Thesis 
and  it  will  not  be  repeated  here. 
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McMillan  and  Opechowaki  apply  (A-6)  and  (A- 7)  to  the  case 
when  each  unperturbed  spin  has  R  nondegenerate  energy  values  , 
a^  >  ■  •  •  ap^  .  Q  pairs  of  these  levels  have  a  separation  hv  ,  and  the 
lower  levels  of  each  of  these  pairs  form  a  set  G  . 

Their  general  formulas  for  the  first  two  moments  in  the  first 
approximation  are  then 
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A24(efg)  =  (e  lm|e  )  (f  |tn|f)  ^  ^ I  |ef  ) 

ijk 

[(e'g|3c!J^^|e'g)  +(f  g|3c!^Vg)  "  («g  l«‘jk^^^|eg) 

-(fg  |K  Vg)l 

A25(efg)  =  l(e|m|e  )|^  ^  [(e  f\y^^}  |e  f)  -  (ef  jy. |cf)] 

ijk 

[(e'gbCji^^^^|c'g)  -  (eg  |yj^^leg)]  , 


where  the  notation  used  is  : 

I  I  ,  ♦ 

|e  )  denotes  the  single  spin  state  of  energy  a^i  *  a^  +  hv  » 

(e  |m|e)  denotes  the  matrix  elennent  between  states  |e)  and  |e  )  of 
the  microwave  field, 

(ef  |y..|ab)  denotes  the  interaction  matrix  element  between  states 
|e).lf)j  and  |a).  |b)j  , 

^  means  summation  over  all  R  states  |e)  of  a  single  spin, 
e 

g  means  summation  over  the  states  in  the  set  G  , 

21  means  summation  over;  all  a  and  b  so  that  •n  +  *  *0  +  *f » 

ab 


A-6 


S 

ab 


means  summation  over;  all  a  and  b  so  that  a_  +  a.  =  a  -fa, 

a  D  e  I 

and  a  is  in  the  set  G  . 


means  summation  over  all  a  and  b  so  that  a_  +  a.  =  a^i  +  a. 

a  b  e  1 

♦ 

=  a  +  a,  +  hv 
e  I 


The  summations 


can  be  simplified  by  remembering 


(A- 11) 


It  is  not  immediately  clear  how  these  formulas  can  be  modified 
to  include  also  the  interaction  between  nonequal  neighbors,  since  this 
requires  a  change  in  the  Boltzmann  factors  as  well  as  the  matrix 
elements. 
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APPENDIX  II 


McMILLAN'S  FORMULA  FOR  THE  FIRST  MOMENT 


IN  THE  SECOND  APPROXIMATION 


McMillan  also  calculated  the  first  moment  to  a  higher  approxi¬ 
mation  by  keeping  the  term  P^3C^^^  P^/kT  in  (A-5).  His  result  is 


h  <  A  V  >2  =  X/Y 


(A- 12) 


where 
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The  aummations  in  A^g  and.A^^  appear  to  be  troublesome, 

E, 

A..  A.  I  ,  -which  will  give  con- 

,  ijkl 

tributions  o£  order  .  The  same  is  the  case  £or  the  term  A  gg  . 

I£  these  terms  remained  in  the  end  result,  it  would  mean  that  the 


correction  from  the  second  approximation  depended  on  N  ,  that  is,  the 
size  of  the  crystal.  However,  these  terms  will  cancel  if  we  expand 


h  <  A  V,  >2 


a  +  b/kT 

7TOTT 


a 

c 


+  ^  - 


ad 

c 


) 


=  h  <  A  V  >j  +  ^  (b  -  ^)  .  (A- 13) 

Even  if  b  and  d  contain  terms  in  that  are  very  large 
compared  to  the  terms  in  N  ,  we  can  still  choose  T  so  large  that 
b/kT  «1  and  d/kT  «  1  ,  and  the  expansion  is  valid. 
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